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Abstract 

In this thesis, I use the strong coupling expansion to investigate the multiflavor lattice 
Schwinger models in the hamiltonian formalism using staggered fermions. In particular, 
I am interested in analysing the mapping of these gauge theories onto quantum spin-1/2 
antiferromagnetic Heisenberg chains. Exploting this mapping, the chiral symmetry breaking 
patterns are studied and the spectra are computed. The extrapolation to zero lattice spacing 
of the results compare favorably qualitatively and quantitatively with the weak coupling 
studies of the gauge models in the continuum. 

I studied the one- flavor lattice Schwinger model, in the strong coupling limit, showing 
that it is equivalent to an Ising model with long-range Coulomb interaction. Even though the 
continuous chiral symmetry is explicitly broken by the staggered fermions, a discrete chiral 
symmetry remains and appears in the lattice theory as a translation by one lattice site which 
forbids fermion mass and must be broken in order for the lattice model to exhibit the features 
of the spectrum of the continuum theory. The effects of the anomaly in the continuum 
appear on the lattice through spontaneous symmetry breaking: indeed I show that the 
ground state spontaneously breaks translations by one lattice site. The dicrete symmetries 
are carefully defined on the lattice. The masses of the low lying bosonic excitations as well 
as the chiral condensate of the model are computed up to the fourth order in the strong 
coupling expansion. Very good agreement between lattice calculations and continuum values 
is found. 

The two-flavor lattice Schwinger model is analysed using the same computational scheme 
adopted in the one-flavor case. The problem of finding the low lying excitations is showed 
to be equivalent to solving the spin-1/2 antiferromagnetic Heisenberg chain. In fact the 
Heisenberg Hamiltonian is the effective Hamiltonian of the gauge model at the second order 
in the strong coupling expansion. The ground state of the spin model represents the vacuum 
of the gauge theory in the strong coupling limit. Two kinds of excitations can be created 
from the ground state; one kind involves only spin flipping and has lower energy since 
no electric flux is created, the other involves fermion transport besides spin flipping and 
thus has higher energy The massless excitations of the gauge model are the spinons of 
the antiferromagnetic chain. The excitation masses can be expressed in terms of spin-spin 
correlators evaluated on the ground state. Either the massless or the massive excitations 
have identical P- and G-parity to the low lying excitations of the continuum two flavor gauge 
model. Good agreement between lattice and continuum results is found. The ground state 
is translationally invariant: both the isoscalar or the isovector chiral condensates are zero 
to every order in the strong coupling expansion as it happens in the continuum theory as 
a consequence of the Coleman theorem. In addition the vacuum expectation value of the 
umklapp operator is non zero as it happens for the corresponding operator in the continuum 
and this is the only relic of the chiral anomaly on the lattice. 

The A/"-flavor lattice Schwinger models are analysed generalizing the N = 2 case and 
at the second order in the strong coupling expansion arc effectively described by spin-1/2 
SU(J\f) antiferromagnetic Heisenberg chains. The generalization is quite straightforward, 
even if the theory is much different depending on if Af is even or odd. 
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This thesis analyses the strong coupling limit of the multiflavor lattice Schwinger models and 
their relationship with generalized quantum spin-1/2 SU(J\f) antifcrromagnctic Hciscnberg 
chains. My original contributions are summaryzed in chapters 4, 5 and 6. In chapter 4 the 
issues of chiral symmetry breaking and the mass spectrum are studied for the one-flavor 
lattice Schwinger model. Part of this work appears in the pubblication 



F. Berruto, G. Grignani, G. W. Semenoff and P. Sodano, Phys. Rev. D57, 5070 
(1998). 



In chapter 5 the two-flavor lattice Schwinger model and its mapping onto the quantum 
spin-1/2 antiferromagnetic Heisenberg chain are analysed in detail. At variance with the 
one-flavor theory, the chiral symmetry breaking is now explicit and the spectrum of the 
gauge model is derived starting from the ground state of the Heisenberg model. Many of 
the results concerning the analysis of this model can be found in 



• F. Berruto, G. Grignani, G. W. Semenoff and P. Sodano, Phys. Rev. D59, 034504 
(1999). 



F. Berruto, G. Grignani, G. W. Semenoff and P. Sodano, hep-th/9901142 



Chapter 6 is devoted to a generalization to the TV-flavor lattice Schwinger models and 
their mapping onto quantum spin-1/2 SU(Af) antiferromagnetic Heisenberg chains. 
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Chapter 1 



Introduction 



Quantum field theory and condensed matter physics had an amazing development in this 
century and each one borrowed from the other new ideas. For instance the concept of 
spontaneous symmetry breaking was originally introduced to study ferromagnetism and the 
Higgs phenomenon appears studing the Meissner effect in the theory of superconductivity £Q . 
Mathematical analogies between the two disciplines were pointed out by Wilson [2] who 
joined the concepts of renormalization, typical of particle physics, with the concepts of 
universality and scaling, originally introduced to describe second order phase transitions. 
Later Wilson |2] introduced a lattice formulation of gauge theories. Once an euclidean field 
theory is defined on the lattice it becomes a statistical mechanics problem. In general a 
quantum theory in d dimensions is equivalent to a classical statistical mechanical system 
in d + 1 dimensions. The phase diagram of a statistical system determines its physical 
properties. Landau introduced the general description of all phase transitions as changes of 
symmetry 0]. Fluctuations grow in a substance as its critical point is approached and they 
interact very strongly with each other so that the system is strongly correlated. The order 
of a phase transition is crucial: only if the phase transition is continuous, it is possible to 
obtain from a lattice field theory a relativistic continuous theory. The lattice approach to 
gauge theories provides a powerful tool to study non-perturbative properties. 

One of the most important properties of non-Abelian gauge theories, such as Quantum 
Chromodynamics (QCD), is asymptotic freedom [5]. At short distances or high momenta 
quarks are weakly interacting and it is possible to perform a perturbative expansion. At large 
distances or small momenta the coupling constant increases indefinitely confining quarks. 
Confinement is an observed property of the strong interactions and it is an unproven, but 
widely believed feature of most non-Abelian gauge theories in four and lower space-time 
dimensions. The theoretical mechanism behind confinement is difficult to analyze since it 
escapes the weak coupling analysis: it is intrinsecally a strong coupling and non-perturbative 
phenomenon. 

Wilson 3, proposed on the lattice a mechanism for quark confinement. Gauge theories 
have been analysed in the strong coupling limit on a (3+I)-dimensional space-time lattice. In 
the weak coupling limit the gauge field behaves like a normal free massless field and quarks 
are unbound. In the strong coupling limit the gauge field is massive and the quarks are 
bound. There should be a confinement-deconfinement phase transition at some intermediate 
value of the coupling constant. 

In the continuum, closely related to confinement there is dynamical chiral symmetry 
breaking, which creates the 7r mesons and keep their masses light. Chiral symmetry 
is only an approximate symmetry of particle physics since the up and down quarks are 



light but not massless. Otherwise the pion, which is the Goldstone boson arising from the 
symmetry breaking would be strictly massless rather than being just light. The theory 
behind chiral symmetry breaking is only partially understood because this phenomenon 
occurs at strong coupling where the standard tools of perturbation theory are not reliable. 
There are conjectures |Z] that if one could adjust the charge of the electron e in ordinary 
QED, there would be a critical value of e at which a phase transition would take place and 
the electron would become very massive. A transition to heavy electrons would break the 
approximate chiral symmetry of QED which is there because electrons are very light. 

The strong coupling limit [8] of lattice gauge theories though far from the scaling regime 
is often used to study the qualitative properties of a gauge model. Two important fea- 
tures of the spectrum of non-abelian gauge theories appear there. The strong coupling limit 
exhibits confinement |3|, which is related in a rather natural way to the gauge symmetry 
and compactness of the non-abelian gauge group. Moreover, it is straightforward to show 
that strongly coupled gauge theories exhibit dynamical chiral symmetry breaking. Wilson 
fcrmions [§] explicitly break chiral symmetry but for what concerns staggered fcrmions 10 , 
even if the continuous chiral symmetry is broken explicitly, a discrete axial symmetry sur- 
vives the lattice regularization and appears in the lattice theory as a translation by one 
spacing. Several quantitative investigations of gauge theories using strong coupling tech- 
niques have been performed and there have been attempts to compute the mass spectrum 
of realistic models such as QCD |11| . The strong coupling expansion is an expansion in the 
inverse of the gauge coupling constant. Due to asymptotic freedom, the continuum limit 
is found where the coupling constant is small, so that one should not expect a priori that 
the strong coupling expansion gives accurate quantitative information about the continuum 
gauge theory. However, the expansion does have a finite radius of convergence J12) : this indi- 
cates that its properties are shared by the model for a finite range of the coupling parameter 
and certain informations can be obtained by analytic continuation to regions outside the 
convergence radius. 

It has been recognised for some time that the strong coupling limit of lattice gauge the- 
ories with dynamical fermions is related to certain quantum spin systems. The relationship 
between gauge and spin systems is particularly evident in the hamiltonian formalism |l.'i| and 
already appeared in some of the earliest analysis of chiral symmetry breaking in the strong 
coupling limit 14 . In particular mesons emerge as spin waves when spontaneous chiral sym- 
metry breaking takes place. It was noticed that there are several formal similarities between 
some condensed matter systems with lattice fermions — and in particular certain antiferro- 
magnetic spin systems — and lattice gauge theories in their strong coupling limit 15 . The 
quantum spin-1/2 Heisenberg antiferromagnet was recognised to be equivalent to the strong 
coupling limit of a E/(l) lattice gauge theory J|)] and can also be written as a kind of SU(2) 
gauge theory |17j . Furthermore, staggered fermions resemble ordinary lattice fermions used 
in tight binding models in condensed matter physics at half-filling, i.e. when there is one 
electron per site, and placed in a background U(l) magnetic field it (mod 2ir) — 1/2 of a 
flux quantum — through every plaqucttc of the lattice. In (2 + l)-dimensions the parallel 
between lattice fermions in condensed matter and gauge theories was already recognized in 
the first work [SI] on the Azbel-Wannier-Hofstaeder problem and then has been discussed 
in the context of the so called flux phases of the Hubbard model. It is actually true for all 
d > 2 + 1 22 , i.e. the staggered lattice fermion approximation of the relativistic Dirac field 
with the minimal number of flavors is identical to a condensed matter fermion problem with 
simple nearest neighbor hopping in a background magnetic field which has 1/2 flux quantum 
per plaquette. In the condensed matter context, the magnetic flux can be produced by a 
condensate, as in the flux phases of the Heisenberg and Hubbard models ^H]- A 1/2 flux 
quantum per plaquette for ordinary lattice spacing is yet an experimentally inaccessible flux 
density. Nonetheless it could be achieved in analog experiments where macroscopic arrays 
of Josephson junctions, for example, take the place of atoms at lattice sites and their ground 
state in two dimensions is expected to be a flux phase 23 . 



The analogy between Dirac and tight-binding condensed matter fermions was exploited 
to hnd an exact mapping of the strong coupling limit of a large class of lattice gauge theories 
with dynamical staggered fermions onto certain quantum antiferromagnets, both the con- 
ventional Heisenberg antiferromagnets with spin S and generalized antiferromagnets with 
spins taking values in Lie algebras other than SU(2) [21]. The precise structure of the 
resulting antiferromagnet depends on the number of colours in the gauge group — since 
S = Af c /2 — and also on the number of fermion flavors. There are four different cases, (a) 
An U(Af c ) gauge theory with an even number Af of lattice flavors of staggered fermions is an 
SU(Af) antiferromagnet in the strong coupling limit. The number of colours Af c determines 
the representation of the antiferromagnet as being the one represented by the Young tableau 
with Af c columns and Af/2 rows, (b) An SU(Af c ) gauge theory with Af c > 3 is an U(Af) 
antiferromagnet. The representation of the U(l) subgroup of U(Af) is not constrained as it 
was in the (a) case where the local U(l) charge had to vanish. The representation of the 
SU(Af) subgroup of U(Af) can now take on any representation on a given site, subject only 
to the constraint that the average fermion density is Af c Af/2 and the number of fermions on 
each site is given by the number of boxes in the corresponding Young tableau, (c) The SU(2) 
lattice gauge theories are special. In the strong coupling limit additional terms compared 
to SU(Af c > 3) appear and one obtains for Af lattice flavors antiferromagnets associated 
with the simplectic group Sp(2Af) and for Af = 1 a spin-1/2 Heisenberg antiferromagnet. 
(d) A gauge theory with an odd number Afoi lattice flavors of staggered fermions in the 
strong coupling limit is an SU(Af) antiferromagnet where the size of representations on 
neighboring sites differ. One has to divide the sites in two sublattices such that the nearest 
neighbors of all the sites of one sublattice are in the other sublattice; if this classification in 
two sublattices is possible, the lattice is called bipartite. On the bipartite lattice one has the 
representation of SU(Af) with Young tableau with (Af + l)/2 rows and Af c columns on one 
sublattice and with (Af — l)/2 rows and Af c columns on the other sublattice. There exist 
analogs in spin systems of strong interactions behavior, particularly for what concerns con- 
finement and chiral symmetry. Laughlin 25 recently provided an appealing interpretation 
of confinement in quantum antiferromagnets. The physical idea is that the phase diagram of 
antiferromagnets consists of competing ordered phases regulated by a nearby quantum crit- 
ical point. Exactly at the critical point the low lying elementary excitations of the magnet 
are gauge fields and particles with fractional quantum numbers analogous to the spinon and 
holon excitations found in spin chains. Away from the critical point, even if very close, the 
"fundamental constituents" bind at low energy scales to make the familiar collective modes 
of ordered states which one renormalizes. The existence of "fundamental constituents" in 
conventional materials and models may be indirectly inferred by high-energy spectroscopy 
and inconsistences in sum rules exactly the way the existence of quarks is inferred in par- 
ticle physics. There exist physically identifiable objects behaving like U(l) quarks out of 
which the elementary excitations are built. In Laughlin's scenario quark-like objects and 
gauge fields would be liberated at the critical point and would become the true elementary 
excitations. 

In quantum antiferromagnets, the appearence of Neel order is the analog of chiral sym- 
metry breaking. Exploiting the mapping existing between gauge and spin models, rigorous 
results about antiferromagnets can be used to prove spontaneous broken chiral symmetry in 
certain strong coupling gauge theories. Whenever the ground state of the antiferromagnet 
is not translationally invariant, the discrete chiral symmetry is spontaneously broken. The 
possibility of studing in detail the patterns of chiral symmetry breaking on the lattice is 
very important. In particular, an interesting question is how the axial anomaly |26| appears 
in the lattice regularization. Anomalies occur when all of the symmetries of a classical 
field theory cannot be realized simultaneously in the corresponding quantum field theory. 
The classical example is the Adler-Bell-Jackiw anomaly in a vector like gauge theory where 
gauge invariance and global axial symmetry are incompatible. This is a manifestation of 
the fact that it is impossible to find a regularization of ultraviolet divergences preserving 
both axial and gauge invariance. If gauge invariance is preserved, the observable manifesta- 
tion of the anomaly is the failure of the conservation laws for axial vector currents and the 



resulting abscence of their consequences in the spectrum of the quantized theory. Nielsen 
and Ninomiya |27j proved that fermion theories on a lattice have an equal number of species 
of left- and right-handed Weyl particles in the continuum limit, thus preventing an axial 
anomaly. Lattice field theories are manifestly gauge invariant by construction and therefore, 
to produce continuum theories, they must find a way to violate axial current conservation. 

In the continuum an interesting mechanism of confinement in Quantum Electrodynamics 
(QED) with massless fermions in (l + l)-dimcnsions was discovered by Schwinger [2S, which 
demonstrated that the "photon" acquires a mass and the spectrum of the model does not 
exhibit any free electron asymptotycal state. The "photon" mass is proportional to the 
electromagnetic coupling constant if charges are totally screened by the vacuum polarization. 
After this investigation QED in (1 + l)-dimensions with massless fermions has been called 
the Schwinger model. Moreover Schwinger speculated that in four dimensions the photon 
mass would be zero for any electromagnetic charge e less than e c ; for e greater than e c the 
photon mass would be non-zero and vary with e. At the critical coupling constant e c a 
change of phase from zero to non-zero photon mass should take place. 

In this thesis it is my purpose to investigate the one-flavor and the multiflavor lattice 
Schwinger models in the hamiltonian formalism using staggered fermions. Exploiting these 
toy models I will clarify the mechanisms of chiral symmetry breaking on the lattice and the 
vestiges of the chiral anomaly of the continuum theory. Moreover the spectra of these models 
will be computed in the strong coupling limit. Extrapolating the spectra to the zero lattice 
spacing limit by means of Pade approximants |29| , I find that the lattice compare favorably 
— both qualitatively and quantitatively — with the weak coupling continuum studies of the 
multiflavor Schwinger models. The success of the procedure is not completely surprising 
since the continuum theory of free massive mesons can be described in terms of quarks and 
antiquarks pairs bound by gauge strings and it is precisely this idea which is embodied in 
lattice methods. The A/"-flavor lattice Schwinger models in the strong coupling limit are 
effectively described by generalized quantum SU(J\f) spin-1/2 Heisenberg antiferromagnetic 
Hamiltonians. Spinons are the fundamental spin-1/2 excitations of quantum antiferromag- 
nets. In the same way as the Schwinger models spectra do not exhibit free quarks (electrons) 
but just massive bosons, the "photons" or better mesons, also the Heisenberg models exci- 
tations contain an even number of spinons and so carry integer spin. One original result of 
my work has been to show explicitly how spinons appear in lattice gauge theories. 

In chapter 2 I shall review the lattice approach to gauge theories. In particular boson, 
fermion and gauge fields are defined on the lattice. The "fermion doubling" problem is 
analysed and Wilson and staggered fermions are introduced. The hamiltonian formulation 
of gauge theories and the strong coupling limit are described. Furthermore the realization of 
the chiral anomaly on the lattice is discussed and the Nielsen and Ninomiya no-go theorem 
is demonstrated. 

Chapter 3 is devoted to a detailed analysis of quantum antiferromagnetic spin chains 
with special emphasis to the topics useful to understand the connection between Schwinger 
and Heisenberg models. First I demonstrate that spin wave theory fails in one dimension 
due to strong infrared divergences. Then I introduce the Haldane-Shastry model where it is 
particularly easy to see the spinons and analyse the spectrum. The Bcthe ansatz solution of 
the spin-1/2 antiferromagnetic Heisenberg Hamiltonian is then reviewed in detail. I study 
the complete spectrum of finite size Heisenberg antiferromagnetic chains of 4 and 6 sites 
and I verify that the thermodynamic limit solution provided in j.'iOj is already very well 
reproduced in the finite size chains. Moreover, I write down explicitly the ground state of 4, 
6 and 8 sites chains. Spin-spin correlators relevant for evaluating the mass spectrum of the 
two-flavor lattice Schwinger model are calculated. In the last section, generalized Heisenberg 
antiferromagnets with symmetry group SU(J\f) are introduced. 

In chapter 4 the one-flavor lattice Schwinger model is studied. The discrete chiral sym- 



metry on the lattice is broken explicitly and this mechanism reproduces the effects that in 
the continuum are generated by the axial anomaly. The discrete symmetries of the model 
are analysed in detail and the vacuum of the model is found to be the ground state of the 
antifcrromagnetic Ising spin chain with long range Coulomb interaction. The spectrum of 
the continuum theory is very well reproduced when the lattice results are extrapolated using 
the Pade approximants. 

In chapter 5 I analyse the two-flavor lattice Schwinger model. In the strong coupling 
limit the gauge model is effectively described by the quantum spin-1/2 antifcrromagnetic 
Heisenberg model. The vacuum of the gauge theory in the strong coupling limit is the ground 
state of the spin model: a quantum disordered state very different from the classical Neel 
state. It is not possible to write down explicitly this ground state in the thermodynamic 
limit, but its energy can be computed using the Bcthe ansatz as explained in chapter 3. The 
massless excitations of the gauge model are the spinons of the Heisenberg antiferromagnet, 
while the massive excitations are created by acting on the ground state with pertinent 
operators that have the right quantum numbers of the continuum theory and generate 
charge transport. The excitation masses of the gauge model are completely expressed in 
terms of vacuum expectation values (VEV's) of spin-spin correlation functions. The chiral 
symmetry breaking pattern is explicit in this case, due to the staggered fermions coupled 
by the gauge field. The isoscalar and isovector chiral condensates, that in the continuum 
are zero due to the Coleman theorem, are zero also on the lattice. The vacuum expectation 
value of the umklapp operator is non-zero due to the explicit breaking of [7^(1); this can be 
viewed as a manifestation of the anomaly on the lattice. 

In chapter 6 the analysis of chapter 5 is generalized to the Schwinger models with SU(Af) 
flavor group. In the strong coupling limit the gauge models are equivalent to generalized 
SU(AT) quantum Heisenberg antifcrromagncts. Even if most of the analysis is a straightfor- 
ward generalization of the M =1 case, some surprising difference arises between the M even 
and odd models. Some concluding remarks and a short summary of results is provided. 

The appendix aims at elucidating some features of Pade approximants relevant for a 
better understanding of chapters 4 and 5. 



Chapter 2 



Lattice gauge theories 



It is my purpose to review in this chapter the lattice approach to gauge theories. I do not 
pretend to be exhaustive, but I want to illustrate the topics used in my approach to the 
lattice Schwinger model. The reader interested in knowing more about lattice gauge theories 
could usefully look up Refs.[5llST]. 

Quantum field theory is so far the most appropriate scheme for describing the strong, 
electromagnetic and weak interactions between elementary particles. It has been known for 
long time that electromagnetic interactions are described by a quantum gauge field theory, 
but the fundamental role played by the principle of gauge invariance in the construction of 
a theory for the strong and weak interactions has been recognized only much later. 

QCD, the theory of strong interactions, is a gauge theory based on the unbroken non- 
abelian SU(3) group. The group SU(3) has eight generators and so there are eight "mass- 
less" gluons carrying a colour charge which mediate the strong interactions between the 
fundamental constituents of matter, the quarks. QCD is an asymptotically free theory 0. 
Asymptotic freedom allowed to carry out quantitative perturbative calculations of observ- 
ables which are sensitive to the short distance structure of QCD. Quarks have never been 
seen free in nature and only colour neutral baryons are observed: quarks are confined. 
A demonstration that QCD accounts for quarks confinement can only come from a non- 
perturbative treatment of the theory; confinement is indeed a consequence of the dynamics 
at large distances where perturbation theory breaks down. 

Wilson [3J showed [J that in the strong coupling limit QCD confines quarks, however 
this is not a justified approximation when studying the continuum limit. There are nu- 
merical simulations which strongly suggest that QCD accounts for quark confinement |81| . 
Moreover, there are other fascinating questions that one would like to answer: does QCD 
account for the observed spectrum? Are there other particles allowed by QCD but not yet 
experimentally observed? Does QCD account for chiral symmetry spontaneous breaking? 
In order to answer the above questions it is necessary a non perturbative treatment of QCD 
which at the moment is only provided by the lattice formulation. 

The lattice formulation of QCD by Wilson opened the way to the study of non-perturbative 
phenomena using numerical methods. His purpose was to study the long distance properties 
of QCD in a format where the short distance properties which lead to ultraviolet divergences 
are regulated by the lattice cutoff. The space-time discretization provides a natural cut-off 

1 Independently, gauge theories were discussed on a lattice by Wegner 1321 . who elevated the global up <-> 
down symmetry of the ordinary Ising to a local symmetry and in an unpublished work, which deals mostly 
with Abelian gauge theories, by Polyakov in 1974. 



scheme as wave-lengths shorter than twice the lattice spacing, a, have no meaning and the 
momentum domain is restricted to a region bounded by n/a. The lattice provides then a 
non-perturbative regularization of ultraviolet divergences. Once a lattice field theory has 
been formulated, the original field theory problem becomes one of statistical mechanics. 
With a finite lattice there are a finite number of variables. It is then possible to study 
various physical interesting quantities: energy spectrum, correlation functions, critical ex- 
ponents etc., for example by computer simulations based on the Monte Carlo method 31 . 
By now lattice gauge theories have become a branch of particle physics in its own right, and 
their intimate connection to statistical mechanics raised the interest of elementary particle 
phisicists as well as of condensed matter physicists. 

Since the beginning of lattice gauge theory, chiral symmetries have been perplexing. The 
issues revolve around anomalies and fermion doubling. For vector-like theories such as QCD 
the problems are largely resolved. The Wilson approach j2j breaks chiral symmetry explicitly 
by adding a symmetry breaking term to give all the doublers a mass which becomes infinite 
with the cut-off scale. The Kogut and Susskind fcrmions JU|, also called staggered fcrmions, 
eliminate the doubling by reducing the Brillouin zone, i.e., by doubling the effective lattice 
spacing. The situation is much more clouded for the full Standard Model. Here chiral 
symmetry plays a fundamental role, with neutrinos maximally violating parity. Nielsen and 
Ninomiya demonstrated a no-go theorem [27| whose most important consequence is that the 
weak interaction cannot be put on the lattice. For a general class of fermion theories on a 
lattice in the hamiltonian formalism an equal number of species of left- and right-handed 
Weyl particles, i.e. neutrinos, necessary appears in the continuum limit. They also proved 
that this disease is not peculiar to the lattice theory, but some trouble appears universally 
in all regularization schemes, for instance dimensional regularization. In any lattice theory 
of chirally invariant fermions with locality there is an equal number of production and 
annihilation of Weyl fermions. Thus there is no net production, so that the axial charges 
are conserved. An analogy or a simulation exists between the Weyl fermion theory and 
gapless semiconductors, where two energy bands have point-like degeneracies. In section 
(|2.4|l I shall discuss thoroughly this no-go theorem. At the moment no reliable lattice 
computation for the Standard Model does exist. 

In section 1)2.1(1 I discuss the differences arising between boson and fermion fields on the 
lattice and I analyse the fermion doubling, an obstacle to put fermions on a lattice. Gauge 
fields on the lattice are discussed. 

Section 1)2.2(1 is devoted to introduce staggered fermions, a particular type of lattice 
fermions which partially avoid the fermion doubling problem. Moreover, hamiltonian for- 
mulation of lattice gauge theories is introduced. 

Section 1(2.3(1 explains the strong coupling approach to lattice gauge theories in the hamil- 
tonian formalism. 

In section ()2.4(l the Nielsen and Ninomiya no-go theorem is reviewed and the chiral 
symmetry and chiral anomaly on the lattice are analysed and the implications and analogs 
of these topics in condensed matter theory are discussed. 



2.1 Bosons and fermions on the lattice 



Any field (scalar, fermionic or gauge) can be defined on a lattice, but the result of putting 
a bosonic or a fermionic field on a lattice is distinctly different. While the procedure for a 
bosonic field is straightforward, placing the Dirac equation on a space-time lattice presents 
some surprisingly difficult problem. 



In this section I shall first illustrate the way to put the scalar field on the lattice. Then I 
shall discuss the so called doubling problem |H] of lattice fermions. A prescription to cure it 
is discussed, the Wilson fermions 9 1. while section l|2.2|) is devoted to discuss the Kogut and 
Susskind fermions or staggered fermions |10| . which have been used in my approach to the 
lattice Schwinger model. The link variable nature of gauge fields on the lattice is explained. 



2.1.1 Bose fields on the lattice: scalar fields 



Let us study how a scalar field 4>{x) can be defined on a lattice. Consider the continuum 
field theory in Euclidean space 



S(<f>) = / d d 



where 



lid^f + V^) 



V{4>) = imV + ^ 4 



(2.1) 



(2.2) 



A matter field is attributed to the lattice sites and it is natural to approximate a continuous 
field by its values at the lattice sites 



(a;) — ► 4> x 



(2.3) 



Clearly, in order for the lattice field 4> x to be a good approximation of a continuous field 
configuration </>(cc), the lattice spacing should be much smaller than the characteristic size 
of this configuration. The derivative can be replaced by 



d u 



1 



J x-\-\i Yx ) 



(2.4) 



where p, is a d-vector of length a in the direction of /i. The d dimensional integration is 
replaced by a sum 

'd d x^a d J2 ( 2 -5) 



/■ 



so that the scalar action on the lattice becomes 

„d-2 d 



S{4>) = £ 



n 






>.,-,, -<:,)-+a d [?Y<l£ + ^4>i 



(2.6) 



It is instructive to go to momentum space to see the spectrum of the free field theory. For 
this I take the Fourier transform 



d d k 



(27T) 



af-'m . 



(2.7) 



Since it is meaningless to consider wavelengths less than twice the lattice spacing, the above 
integral is taken over only one "Brillouin zone" of the reciprocal lattice 



7T 7T 

- < L < - for each [i 
a a 



(2. 



where k ll = k- jX. Substituting Eq. (|2.7|) into Eq. (|2.6() the free field action can be written as 



s ^ )= \Swf 



E4 . o / ak u 



<K-k)<Kk) 



(2.9) 



Each mode contributes to the action in the momentum space a quantity 



(2.10) 



S(k) - m 2 



-7r/a 



7r/a 



Figure 2.1: The dispersion relation S'(fe) for a free scalar field. The solid line, fc 2 , is for the 
continuum theory, the dotted line for the latticized system. 

rather than the standard m 2 + k 2 . Nevertheless these two expressions have the same con- 
tinuum limit as they coincide at the minimum value of k = (see fig. I|2.1|l 'l. 



The lattice action Eq. (|2.6() is quantized by using the Feynman path- integral formalism, 
in which the expectation value of a product of fields is given by 



< o\cf> xt 4> X3 . . . 4> xl \o >= — / J^[rf0x]^!^ 2 



-S(<j>) 



(2.11) 



where 



Z 



Hid '4> x 



-s(fl 



(2.12) 



The meaning of the integrals should be clear as one recalls that the usual functional integrals 
are actually defined on a discretized space-time lattice and an appropriate continuum limit 
is taken at the end. If one rescales the field as 



the lattice action scales as 



b' x = V\(f> x , 



5(0) = ~SV) 



with 



W) = E 



~ij- Z^Wx+it - (p x ) +a \—4> x + -q> x 






(2.13) 

(2.14) 

(2.15) 



i.e. the coupling constant A has become an overall factor in the action. In this way Eas. (|2.11|l 
and (|2.12|) may be written as 



< o|<M 



■■ 4 ! io>=|/n[^x]c^ 

X 



X2 . . . ifj xi e 



(2.16) 



Note that Ea. l|2.16|) has the same structure of the partition function in statistical mechanics, 
once the identification 

\^ Ps ir (2 ' 17) 

is made. The strong coupling expansion (in powers of A -1 ) corresponds to the high temper- 
ature expansion in statistical mechanics. 



2.1.2 Fermi fields on the lattice 



I shall now introduce the procedure to put Fermi fields on the lattice. A new phenomenon, 
the "fermion doubling" 8 , 31 appears discretizing the Dirac equation, i.e. additional fermionic 
species are generated. Let us illustrate the fermion doubling problem starting from an exam- 
ple on a spatial lattice and continuum time in 1+1 dimensions. The massless Dirac equation 
in 1+1 dimensions in the continuum reads 



with -0 



4>i 

4 ; 2 



iip = -iad x ip = -ij5d x ip 
and a = 75 = o\ . Choosing for the plane waves 

^±(z,i)=e- i(fcx - Bt) X± 



with 

the dispersion relation is 



75X± = ±X± 



(2.18) 

(2.19) 
(2.20) 



E{k)=±k , -oo</c<oo (2.21) 

and the excitations are left- and right-handed particles and antiparticles — see fig. (|2.2(l . 




Figure 2.2: Spectrum of the free, continuum Dirac fermion in 1 + 1 dimensions. 



Let us consider the naively latticized Dirac equation. I place the spinor tp = I , | on 
each site of a spatial lattice and change the partial derivative d x with a finite difference 



iijj(x) 



2« 



7S[V»(s + l)-ip(x-l)] 



(2.22) 
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The dispersion relation fig. (|2. 3(1 is now 



E{k) = ± 



sin(fco) 



(2.23) 




Figure 2.3: Spectrum of the naively latticized Dirac fermion in 1 + 1 dimensions. 



This is a relativistic dispersion relation for i:a< 1 

E(k) ~ ±k + 0(k 3 a 2 ) 
and for ka = n — k' a, with k' a <C 1, 

E(k) ~=Ffc' + 0(fc' 3 a 2 ) 



(2.24) 



(2.25) 



There are two two-component Dirac particles in the continuum limit, and the total chiral 
charge of the fermions is zero. Consider, for example, the excitations associated with the 
field ip+ on the lattice. There are right-movers (k ~ 0) and left-movers (fc ~ 7r), and since 
chirality (helicity for particles) is just velocity in 1 + 1 dimensions this lattice field describes 
a pair of fermions with net chirality zero. 

The same phenomenon happens in 3+1 dimensions, even if it is more difficult to visualize 
it. The Dirac Hamiltonian for massive fermions now reads 



H = 



■ d 

-iai— 1- rap 

oxi 



(2.26) 



with 



a; =757* , P = 7° • (2.27) 

The field operator describing the multifermion system will be denoted by ipa(x) where 
a = 0, 1,2, 3 is the Dirac index and the three-component vector x denotes a space point. 
The Hamiltonian of the free fermion field is obtained from the single particle Hamiltonian 
Eq.||H§ as 

H = [ d 3 xiPi{x)H a0 Mx) (2-28) 
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or rewriting Ea. (|2.28J) with a symmetrized derivative, one gets 

H = J#x4>\(x)lmp+±(^--^)a l \ M*) ■ ( 2 - 29 ) 

Let us now introduce a regular cubic lattice in space with lattice spacing a and elementary 
cubes of volume a 3 . On a finite lattice with total volume V = (aL) 3 the lattice points have 
integer indices x = {xi ; I = 1, 2, 3} satisfying < xi < L — 1. The continuous field will 
be replaced by a discrete variable associated with elementary cubes, and the volume integral 
becomes a finite sum 



iW^ii , d 3 x^J2 a3 ■ ( 2 - 3 °) 

J X 

The discrete field variables are assumed to satisfy the anticommutation relations 

{VW, i>Hy\ = {ll>tg, ^Iff} = , [iplg, V>w} = ^a/3<% ( 2 -31) 

and the continuous derivatives are discretized as 

d 1 

— i/> a (g) -* -{ip a3+ j;~ ii) (2.32) 

where I denotes a unit vector on the lattice in direction I. Applying these rules in the 
Hamiltonian Eq. (|2.29|l one obtains the naive discretized form 

H = J2 a ) ^\s m Msx + Y a Yl^ls+M^ 5 ^ ~ ^d^h^ss+i \ ■ ( 2 - 33 ) 
2 I 1=1 ) 

The free field Hamiltonian describes a collection of harmonic fcrmionic oscillators in mo- 
mentum space. Assuming periodic boundary conditions in all three orthogonal directions of 
the L 3 cube, the allowed momentum components are discrete points in the Brillouin zone, 
namely 

p=ipi = ^m ; I = 1,2,31 , 0<m<L-l . (2.34) 

Due to periodicity, instead of < pi < 2tt, the Brillouin zone can also be equivalently defined 
in the interval — n < pi < tt. Fourier transformed field components read 

X 

VW = 7^pE e ^ Xi ^P ' ( 2 - 36 ) 

v 

In momentum space the Hamiltonian lj2.33|l reads 

H = , , 3 ^2 ^\p \ m P+ ~ a i s[n Pi \ ^sp ■ (2.37) 

The matrix in curly brackets in Ea. l|2.37(l has eigenvalues 



E f =± 



\ 



1 3 
— ^sin 2 P/ . (2.38) 



;=i 



Ea. (|2.38|l gives the possible energies of the free naive lattice fermions. In the continuum 
limit a — » 0, Eq. (|2.38|l approaches the usual dispersion relation for free fermions, since the 
coordinate vector is ax and correspondingly, the momentum vector is k = -p. For fixed 
momentum k and a —* one has 

Ef y = m 2 + k 2 + 0{a 2 ) . (2.39) 
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This is as expected, but still the naive lattice Hamiltonian Ea. (|2.33|l is not entirely satisfac- 
tory. Due to the fact that sin(p + w) = — ship one has 

El = E} + ~ (2.40) 

where if is one of the eight vectors 

tt = { (0,0,0),(7r,0,0),(0,7r,0),(0,0,7r) , 

(7r ) 7r J 0),(7r,0,7r),(0,7r J 7r),(7r,7r,7r) } . (2.41) 

These are the eight corners of the part of the Brillouin zone satisfying < pi < ir (I = 1,2, 3). 
The consequence of Eq. (|2.4()|> is that in this naive lattice formulation there are eight fermion 
states per field component. This phenomena is called fermion doubling although a more 
appropriate name would be fermion octupling. The basic reason for fermion doubling is 
that the Dirac equation is of first order. In the case of free fermions one could, perhaps, 
tolerate this unwanted proliferation of degrees of freedom, but in an interacting theory the 
extra fermions influence the physical content in a non-trivial way, because the additional 
states can be pair-produced through the interactions of the fermion field. For example, even 
if the external particles in a process are the states at the zero (0, 0, 0) corner of the Brillouin 
zone, the states of the other zones appear in virtual loops. 

In order to cure the disease, Wilson introduced an additional second order term in the 
lattice Hamiltonian 

H = J2a 3 {ipl s (m+—)l3 7 s'ipsx 

X 

3 
I 



Ya J2^l s+ i( a M^ ~ ^\M)i^&3+^ 



1=1 

3 



YY<&l S+ ^M8x+iPUPM ss J} ■ (2.42) 



2__ 
i=i 



The Wilson parameter r is assumed to be in the interval < r < 1. In momentum space, 
instead of Eq. Q2.37J I one has now 



y j p i=i 



pi) 



a 
and the correspondingly energy eigenvalues are 



+ -aismpi} 7 fip S p (2.43) 



Eff = ± 



\ 



3 

r 



]T(l-cos W )] 2 + E4 sin2 ^ ■ (2-44) 



a * — ' ' — ' a 

i=i i=i 



In the continuum limit a — ► the mass m is replaced by m + — n n , where n, is the number 
of momentum components equal to n. Therefore, the states with n^ ^ become infinitely 
heavy. The only physical fermion state with finite energy is at the zero corner of the Brillouin 
zone. In section Ij2.2|l I shall provide a detailed description of a different method to cure 
the doubling disease, i.e. the staggered fermions which are the closest possible analog to the 
lattice fermions in condensed matter physics. 



2.1.3 Gauge fields on the lattice 



Let us now discuss how to put gauge fields on the lattice and briefly review how one arrives at 
the gauge invariant action in the continuum with the simple example of QED. The starting 
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point is the action of the free Dirac field 

S F = ( &i$[x){i'fd li -rn)<il>{x) . (2.45) 

The action is invariant under the transformation 

ip(x) — ► e lA %j)(x) (2.46) 

${x) — > 4>{x)e- lA (2.47) 

with A independent of x (i.e. a global transformation). 

The next step is to require to the action to be also invariant under local U(l) trans- 
formations with ip (ip) transforming independently at different space-time points. This is 
accomplished by introducing a four- vector potential A^x) and replacing the ordinary four 
derivative d^ by the covariant derivative T)^ defined by 

V^ = d^ + ieA^ . (2.48) 

The resulting new action 

S F = I d i x^(x){i-/ l2 V l " - m)%l){x) (2.49) 

is then invariant under the local transformations 

i/)(x) — ► e iA(x) ip(x) (2.50) 

^{x) — ► ^(x)e-' A ^ (2.51) 

A^x) — * e iA ^(A fl (x)--d fi )e- iA ^ . (2.52) 



The crucial property which ensures the gauge invariance of Ea. H2.49(l is that, while A^ 
transforms inhomogeneously, the transformation law for the covariant derivative Eq. (|2.48|l 
is homogeneous 

X> M — ► e lA{x) V^ lA{x) . (2.53) 

By introducing a four-vector field A^ gauge invariance of the action (|2.49(1 is ensured; one 
now must introduce a kinetic term which allows A^ to propagate 

S G = -\ j d A xF^F» v (2.54) 

F^u is the gauge invariant field strenght tensor 

F„ v = d„A v - d v A^ . (2.55) 

The gauge invariant action describing the dynamics of the fermionic and gauge field is 

Sqed = S f + S g . (2.56) 



Let us start our approach to lattice QED by considering the lattice action of a free 
Dirac field. To parallel as closely as possible the steps in the continuum formulation, I shall 
consider Wilson fermions where every lattice site is occupied by all Dirac components ip a . 
The results presented in the following are valid also for staggered fermions. The Hamiltonian 
Ea. (|2.42|) is invariant under the global transformations 

i> ax — ► e iA ^ ax (2.57) 

^ ax — ? Q ,e- jA (2.58) 

where e is an element of the U(l) group. The next step is in requiring the theory to 
be invariant under local 17(1) transformations, with the group element e " depending on 
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the lattice site. Because of the non diagonal structure of the second and third term in the 
Hamiltonian Ea. H2.42(l (whose origin is the derivative in the continuum formulation) one is 
forced to introduce new degrees of freedom. Since the group elements e x do not act on the 
Dirac indices, it is sufficient for the following argument to focus our attention on a typical 
bilinear term 

«H + f • ( 2 - 59 ) 

In the continuum it is well known how such bilinear terms should be modified in order to 
arrive at a gauge invariant expression. Since the bilinear ip(x)ij)(y) transforms under gauge 
transformations 

^{x)ip(y) — ^(x)e- iA ^ ) e iA ^ip(y) (2.60) 

one must include a factor depending on the gauge potential which compensates the Eq. (12.6011 
gauge variation. This factor, called Schwinger line integral is given by 

U(x,y)=e ie S* dz » A,i W (2.61) 

and the line integral is carried out along a path connecting the point x with the point y. 
U(x, y) is an element of the U(l) group. Under a gauge transformation Ea. (|2.52[l . Ea. l|2.61|l 
transforms as 

U(x,y) -^e lK{s) U(x,y)e- lK{S) . (2.62) 

From these considerations, one concludes that the following bilinear expression in the fermion 
fields ip and ip is gauge invariant 

^(x)U{x iy :)ij{y)=^{x)e le ^ dz ^ {s H{y) . (2.63) 

Let us suppose now that y = S+e. From Eq. (|2.63J) one concludes that the bilinears should 
be defined as 

4>{x)ij)(x + t) ► i)(x)U(x,x + e)i)(x + e) (2.64) 

^{x + t)ip{x) — ► lp(x + t)U\x 1 x + t)ip{x) (2.65) 

with 

U{x,x + e) = e le? - A ^ (2.66) 

and e-A = Y^^. 

To arrive at a gauge-invariant expression for the fermionic action on the lattice, one 
should make the following substitution in Ea. H2.42(l 

il> s+ j{loli + ir)ips — ► ip s+ [iloli + ir )U 3+ i 3 Tpx (2.67) 

i>Alolfi - ir)tl> s+ r — ► i^dloli - W^+jl^+f (2.68) 

where 

and U# -, ris an element of the U(l) gauge group and can therefore be written 



where <fr l s is restricted to the compact domain [0, 2tt]. The right hand side of Eas. (|2.67l2.68|l 
are now invariant under the following set of local transformations 

(2.71) 
(2.72) 
(2.73) 
(2.74) 
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4>x - 


-» e iA ^ s 


iPs ~ 


-> ^e"^ 


x : x-\-l 


- e iAg U g ^- iA *+r 


x-\-l,x 


-» e iA *+rU^ r ~e- iA * 



a 



x,x+l 

> — 



u 



x+l,x 



X + l 



O- 

x 



x + l 



Figure 2.4: Link variables. 



At variance with the matter fields discussed before, the group element U- -,,- lives on the 
links connecting two neighbouring lattice sites; hence I shall refer to them as link vari- 
ables. The link variables have a direction and one can introduce a graphical representation 
- fig.JHJ). 

The Hamiltonian Ea. 1)2.42(1 with the Wilson fermions can be written in the following 
gauge invariant form 

3r\ 



H = VV{ (m+—)i/) aS ip aS 



•h Y} ^as+^oli + ir) al3 U s+r ^ ps 



2a 



1=1 



i> a dlOll ~ i^apUg^f/)^} 



(2.75) 



By requiring that U s -, ,-and U~,?-. transform according to Eas. (|2.73l2.74|) . one naturally 
implements [/(l)-gauge invariance in the Hamiltonian and, a posteriori, one can state that 
the link variables are elements of the 27(1) gauge group due to the requirement that their 
gauge transforms must also be elements of U(l). 

Taking the continuum limit a — > of Eq. (|2.75fl one should get the usual continuum 
Hamiltonian Ea. ((2.29(l . This requirement is fulfilled by establishing a relation between the 
link variables and the vector potential A^x). The vector potential A^(x) at the lattice site 
x is real valued and carries a Lorentz index. The same happens for $- which parametrizes 
the link variable Ug^+a- The difference is that $- takes only values in the interval [0, 2tt], 
while the values taken by the vector potential A^ix) in the continuum theory extends over 
the entire real line. 

Since A^ carries the dimension of an inverse length, while $ M is dimensionless one may 
try the ansatz 



$ M (x) = caA tl (x) 



(2.76) 



where a is the lattice spacing and c is a constant to be determined. In the continuum limit 
a — > the range of A^ will be infinite. 

By scaling m,f/j and ip as m — » am, ip — ► a= ip and ip — > a zip and replacing U^^+p f° r 
small values of a with 

Ux,s+p — 1 + icaA^x) . (2.77) 

Ea. ((2.75(1 reduces to the continuum Eg. 1(2.49(1 if one sets c — e. From here on I use the 
notation 

U„(x) = U^+it = e ieaA -W . (2.78) 

With this identification it is now an easy matter to verify that for a->0 U^(x) transforms 
as follows under gauge transformations 



U^x) -» e <A *J7 /t (2)e- <A »+" = e ieaA ^ s) 
where A~[(x) comes from the discretized version of Eo. ((2.52(l . 



(2.79) 
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To complete the construction of lattice QED one has to construct the lattice version of 
the kinetic term Ea. lj2.54jl for gauge fields which should be strictly gauge invariant and a 
functional of link variables only. Such gauge-invariant functionals are easily constructed by 
taking the product of link variables around closed loops on the euclidean space-time lattice. 
Moreover, because of the local structure of the integrand in Ea. (|2.54|l one should focus his 
attention on the smallest possible loops on the euclidean space time lattice. Hence one is 
led to consider the product of link variables around an elementary plaquette as shown in 
fig. l|2.5|) . One then defines 



x + v 



Vtix) 



X + fl + v 




/\ U v {x + fl) 



Figure 2.5: A plaquette. 



U^(x) = U„(g)U v {g + p)UUx + V)Ut{x) 



(2.80) 



where I have path ordered the link variables. The path ordering is irrelevant in the Abelian 
case, but becomes important in QCD. Inserting Eq. (|2.79|l into Ea. H2.80|) one gets 



U^(x) 



a i&d 1 F tl , v (x) 



(2.81) 



where F^ (x) is a discretized version of the continuum field strength tensor 



F^(S) = ~[(A„(£+p)-A v (S))-{A li (x + i7)-A fi {x))} ■ ( 2 - 82 ) 

From Ea. H2.81|) . in the limit of small lattice spacing a, one has 

i£ £ l 1 - \( U ^) + U U S ))] - I E a 4 F^(x)F^(x) (2.83) 



x H<v 



where the sum appearing on the left hand side extends over the contributions coming from 
all distinct plaquettes on the lattice. The lattice action [3] for the gauge potential in a 
compact form reads 

(2.84) 



^oo = ^z{i-^+4)} 



where Up, the plaquette variable, stands for the product of link variables around the bound- 
ary of the plaquette P, taken in a counterclockwise direction. 
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2.2 Staggered fermions 

In what follows I shall introduce the staggered fermion formalism JO] that I used to study 
the lattice Schwinger models. In this formalism one eliminates the unwanted fermion modes 
by reducing the Brillouin zone, i.e. by doubling the effective lattice spacing. In the case of 
a (i-dimensional lattice, one subdivides it into elementary <i-dimensional hypercubes of unit 
length. At each site within a given hypercube one places a different degree of freedom and 
repeats this structure periodically throughout the lattice. Since there are 2 d sites within a 
hypercube, but only 2? components of a Dirac field, one needs 2^ different Dirac fields to 
reduce the Brillouin zone by a factor of i 

In the continuum, the massless Dirac spinors are invariant under the chiral rotations 

^(i)^e ils V(i) , ^(x) — > Jj{x)e~ l ^ e (2.85) 

The continuum axial symmetry l|2.85|) is broken explicitly by the staggered fermions, but a 
discrete axial symmetry remains. It is a chiral rotation of 9 — tt/2 and corresponds to the 
continuum transformation 

i/j(x) — > 75^) i ^(x) — ► ~ip(x)j5 (2-86) 

and appears on the lattice as a translation by one site. 

Let us start to see how staggered fermions work in (1 + I)-dimensions. The main idea is 
to reduce the number of degrees of freedom by using a single component Fermi field ip x on 
each site of the lattice 

ty*,lM = , {1>tM = S xv (2.87) 

so that the lattice Dirac Hamiltonian reads 

^ = -£E^+i^-^+i) (2-88) 

and the equations of motion are 

ipx = -i[H D ,ip x ] = tH^z+i _ Tpx-i] ■ (2.89) 

la 

If one decomposes the lattice into even and odd sublattices (characterized by x even and 
odd), one can identify a single two component Dirac spinor by associating the upper com- 
ponent with even sites and the lower component with odd sites (or vice versa). In this way 
the equations of motion become 

ipi{x) = ^-[i>2{x + l)-fa{x-l)\ (2.90) 

la 

Mx) ^ ^-[Mx + I) - Mx - I)} ■ (2.91) 

la 

Comparing Eq. (J2.90JI and Eq. (|2.91Jl with the the continuum Dirac equation, one immediately 
sees that the staggered fermion formalism avoids the species doubling problem by doubling 
the lattice spacing and so halving the size of the Brillouin zone. 

Let us now discuss the symmetries of the lattice Dirac Hamiltonian Ea. H2.88|l . It is 
invariant under translations of the spatial lattice by an even number of sites. Translation 
by two lattice spacings is the ordinary continuum translation, whose generator is 



P 



dxip^d x ip = -i dx{ip\d x ipi + ^\d x i>i) (2.92) 
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and it does not mix upper and lower spinor components. Consequently, one expects the 
lattice generator not to mix the two sublattices 

p = -i ^(4+2V>z +V4V>*+ 2 ) • (2.93) 



The Hamiltonian Ea. H2.88|) is also invariant under translations by an odd number of 
sites. The lattice generator of this symmetry is 

Qs = J>4V^+i + 4+i^) (2.94) 

X 

which in the continuum would be Q5 = J dx^j^. The matrix 75 applied to the bispinor 
interchanges upper and lower component, this on the lattice corresponds to translation by 
one lattice site. 

I illustrate now staggered fermions on a (rf - l)-dimensional hypercubic lattice and con- 
tinuum time ( hamiltonian formalism) which are obtained by spin-diagonalization |3$j of 
the naively latticized Dirac Hamiltonian 



/ = - ~ E (^ ( f + 3)°i*> ^ ( f ) " ^ (*)<*' ^ ( f + 3)) (2-95) 



H 

2 

x,j 



where ip a (x) are the fermion field operators, a — 1, 2, . . . 2^ d ' 2 > is the spin index, x = 
(xi,x<z . . .X(d-i)) = ^2j=i x jj (the Xj are integers) refers to lattice sites, j are unit vec- 
tors and a J are the 2^- d / 2 ^ x 2^- d / 2 ^ Hermitean Dirac matrices - here [d/2] is the largest integer 
less than or equal to d/2. Hamiltonian (|2.95|) describes a fermion hopping problem in a 
Tj(2[ d / 2 1) background gauge field given by the unitarjQ matrices a- 7 . The crucial observation 
which allows spin diagonalization is that this background field has only U(l) curvature, i.e. 
if one considers the product around any plaquette 

a J a fe (a i ) t (a* : ) t = -1 ■ (2.96) 

Hence the only information carried by the a-matrices which is invariant under a space 
dependent change of phase of the fermion fields is that a product of a's around any plaquette 
is —1. This allows diagonalization of the a's using a gauge transformation such as 

iP(x) -> (a 1 )* 1 (a 2 )* 2 . . . (a^-Vy^-ViPix) (2.97) 

resulting in the Hamiltonian 

Hf = -^(-l) 5 ^ 1 " (^(x + jH(x)-^(x)^j(x+j)) (2.98) 

which describes 2' ld / 2 \ species of lattice fermions with background U(l) magnetic flux n 
through every plaquette of the lattice. This flux is contained in the background U(l) field 

U°[x,j] = (— iy^p =1 Xp . Each species of fermion must have the same spectrum as the 
original one given by the Dirac Hamiltonian Ij2.95|) . This allows reduction of the fermion 
multiplicity by a factor of 2'- d ' 2 ' (by dropping the fermion spin index a). The result resembles 
a condensed matter hopping problem with a single species of fermion where there is a 
background magnetic field 7r per plaquette. 

2 Dirac's a-matrices are Hermitean, 

a't = a* , 

and, due to their anticommutator algebra, 

{0%^} =25 ij , 
they are also unitary. 
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Discrete Chiral symmetry interchanges the even (Y) p _ j a; p =even) and odd Q2p=i x v — °dd) 
sublattices. The substitutions 

ip(S) -> (-l)^t=hi x ^(x + j) (2.99) 

for j = 1, . . . , (d — 1) leave the Hamiltonian in Ea. (|2.98() invariant. A candidate for Dirac 
mass operator, which changes sign under the transformations in Eq. i|2.99|) . is the staggered 
charge density operator 

fi = m'Y^{-l)^i=^ x '-ijj\x)xjj{x) . (2.100) 



To obtain the continuum limit and the number of fcrmion species, one first divides the 
lattice into 2 d ~ 1 sublattices according to whether the components of their coordinates are 
even or odd. For example, when (d — 1) = 3, I label 8 fermion species as 

■(/'(even, even, even) = ip\ , ip(even, odd, odd) = ip2 
0>(odd,even, odd) = ipz , ^(odd, odd, even) = -04 
■(/'(even, even, odd) = -05 , ip (even, odd, even) = ipe 
^(odd, even, even) = -07 , ^(odd, odd, odd) = tps • (2.101) 

Then, if one adds the mass operator 1)2.100(1 . in momentum space the Hamiltonian H2.98J1 
has the form 



with 

and 8x8 matrices 



H f = / d?k ^(k) {A 1 sinfej + Bm) ip(k) (2.102) 

Jn B 

^(k) = (4(k),...,4(k)) (2.103) 



with a 1 the Pauli matrices; I used the Fourier transform 

W) = J aB ^-^e-^m (2.105) 

where £Ib = {k = (fei, fe, ks), —tt/2 < ki < it/ 2} is the Brillouin zone of the (even, even, even) 
sublattice. The fermion spectrum is 



w(Jfe) 



. ^sin 2 fc 4 +rn 2 , (2.106) 



and only the region ki ~ is relevant to the continuum limit. I have normalized ?/>(fc) so 
that 

{i>{x),^{y)} = 8{x-y) , {^(k),^(l)} = S 3 (k- I) . (2.107) 



a 2 
-o< 



If one defines f3 = ( n 2 ) ancl the unitary matrix 



and 



— Ki;2 !-?' 



ip = Mip' (2.109) 
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with 



the Hamiltonian is 



Hi = / d 3 k 



# = ( 'i i . ( 2 - n °) 



</>. 



/Or 



V'" n M a 1 sinfc + /3m / V </>& / V 



the low momentum limit, sin /e^ ~ fcj, with fermion density 1/2 per site so that the Fermi 
level is at the intersection point of the positive and negative energy bands, one obtains 2 
continuum Dirac fermions. Furthermore, the staggered charge operator gives a Dirac mass 
of differing sign for the two species. 

If one considers AfAfc lattice species in d dimensions, i.e. fermions ip 1 ^ with a = 1, . . . , TV 
and A — 1, . . . ,Mc a lattice flavor and color index, respectively, this yields AfAf c 2^ d ~ 1 ' ) /2^ d ^ 2 ^ 
continuum species of Dirac fermions where the lattice fermion density is AfAfc/2 per site. 
For the Dirac mass operator one may choose 

» = ^(-l)^ Z{x *^{x)rn a p^ B {x) , (2.112) 

X 

with m a p = diag(mi, . . . , to/vOj and this gives the fermion spectrum 

oj(k) 



«/-l 



. >J sin 2 fci + m Q (2.113) 



where — ir/2 < ki < ir/2 and a = 1, . . . ,Af. 

Here, I count the number of flavors of fermions obtained in the continuum limit by noting 
that Ea. (|2. 98(1 describes a 2 d ~ ^component fermion. In d dimensions the Dirac matrices are 
[d/2] dimensional, therefore the continuum limit of Ea. lj2.98fl describes 2 d ~ 1 /2'- d ' 2 ' species 
of Dirac fermions. 

d no. of flavors no. of spinor components 

2 1 2 

3 2 2 ( 2 - 114 ) 

4 2 4 



2.3 Strong coupling lattice gauge theory 



This section is devoted to a detailed analysis of the strong coupling limit of hamiltonian lat- 
tice gauge theories. I shall give a general introduction to the subject by considering a gauge 
theory with colour group U{Nc) or SU(Afc) and S'?7(A/')-fiavor groups. The Hamiltonian 
formalism is particularly suitable for illustrating the relationship between several gauge and 
spin systems and was already exploited in some of the earliest studies of chiral symmetry 
breaking in the strong coupling limit |15) . The strong coupling limit, though far from the 
scaling regime is often used to study qualitative properties related to confinement and chi- 
ral symmetry breaking. In chapter 4, 5 and 6 I shall study the one-flavor and multiflavor 
lattice Schwinger models. The one-flavor model will be mapped in the strong coupling limit 
onto an antiferromagnetic Ising model with long range Coulomb interaction. The multiflavor 
Schwinger models are effectively described by spin-1/2 antiferromagnetic S'C/(A/') Heisenberg 
models. 
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2.3.1 Hamiltonian formulation of lattice gauge theory 

Let us consider an hamiltonian lattice gauge theory \l'6\ where space is a (d— l)-dimensional 
hypercubic lattice with oriented plaquettes [a:, i, j] with corners x, x + i, x + i + j, x + j and 
orientation i x j. The gauge field U[x, i] associated with the link [x, i] is a group element 
in the fundamental representation [Afc] of SU(Afc) and, if the color group is U(Afc), also 
carries a representation of U(l). It has the property U[x,—i] = U'[x — i,i]. The electric 
field operator E a [x,i] associated with link [x,i] has the Lie algebra 

[E a [x,i],E b [y,j]] = if abc E c [x,i]5([x,i]-[y,j]) (2.115) 

and 

E[x,-i] = -tf[x -i,i]E[x -i,i]U[x -%,i] (2.116) 

where 

E[x,i] = E a [x,i]T a (2.117) 

with T a = (T a ) t , a = . . . ,7V^ - 1 the generators of the Lie algebra of U(Nc) obeying 

[T a ,T b ] =if abc T c , (2.118) 

T° = 1 the Afc x ■N'c un it matrix representing U(l) and T a , a ^ 0, in the representation 
[Afc] of SU(Afc)- It generates the left-action of the Lie algebra on U[x, i), i.e., 

[E a [x,i],U[yJ]] =-T a U[x,i]5{[x,i]-[y,j]) . (2.119) 



The Hamiltonian for lattice gauge theory with staggered fermions and discrete chiral 
symmetry is 

H= J2 ^E a [x,t} 2 + I J^ ±Tr(U[x,tW[x + i,flU*[x+j,i\U*[x,3\) + h.. 

[x,i],a \[x,i,j] 

J2tix,i]*l>A(x + i)UAB[x,i}tl>%(x)+h.c.\ (2.120) 

[x,i] J 

where Tr(-) is the N'c x Afc matrix trace and where 

*M = 4H EJ " iXp ' h*,-i]= t * [x -iA {2A21 ^ 

and the generators of static gauge transformations are 

(d-i) 
G a {x)= Y, E a [x,i}+^(x)n B ^( X ) (2.122) 

i=-(d-l) 

with o=l,... ,Af c — 1. If the gauge group is U(Nc) rather than SU(Nc) there is also the 
U(l) generator 

(d-i) 
Q°(x)= J2 {E [x,i}-E°[x-i,i})+-[^(xU%(x)] (2.123) 

with the fermionic U(l) charge operator ordered so that, like E, it changes sign under charge 
conjugation. (Note that there is no ordering ambiguity for the Q a , a ^ 0.) These generators 
obey the Lie algebra 



[G a (x),g b (y)}=if abc g c (x)6(x-y) . (2.124) 
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The Hamiltonian is gauge invariant, 

[G a {x),H]=Q . (2.125) 

The dynamical problem of lattice gauge theory is to find the eigenstates of the Hamilto- 
nian operator Eq. (|2.120|l which are also gauge invariant, i.e. which obey the physical state 
condition 

G a (x)\^ phys > - . (2.126) 



2.3.2 Strong coupling 

In order to perform a strong-coupling expansion |31j . let us rewrite Ea. (|2.12()|l 

H = Ho + Hi + H 2 

where H = £ ^-E 2 , H x = ^(t^U^ + h-c.) and H 2 = £ ^1 {TrUUUU + h.c). Each term 
is gauge invariant, 

[g a (x),H t ] = o. 

Therefore, if one finds a gauge invariant eigenstate of Hq, perturbations in Hi and H 2 remain 
gauge invariant. Hq is the sum of group manifold Laplacians for each link. If |0 > is a singlet 
of the algebra l|2.115[) . i-e. 

E a [x,i}\0>=0, 

then 

H Q \0>=0 

and 

Q 2 

H o U[x,i]\0 >= y —C 2 (N c )U[x,i\\0 > (2.127) 

In Ea. (|2.127Jl the Casimir operator is 

c 2 (M c )i = J2 TaTa 

a 

and a runs either from (for U{N C )) or 1 (for SU{Mc)) to Mq - 1. 

Consider the empty vacuum which is a singlet of 1)2.115(1 and which has no fermions, 

V>2(x)|0 >= , Vx,A,a. 

It is necessary to find the lowest energy (Eq= eigenvalue of Ho) eigenstates 1^ > of Hq 
which are gauge invariant, 

^ Q (a;)|*>=0 Va,x, 

with the constraint that the fermion states are half-filled, i.e. half of the fermion states are 
filled. 

I shall first review the simplest case of SU(Afc > 3) gauge theory with 1 lattice flavor 
(TV = 1) and 2 d /2[ d / 2 l continuum flavors JT]. The density is A/"c/2 fermions per site. The 
lowest energy eigenstates of Hq with Eq = are the states which are singlets of the electric 
field algebra (|2.115|1 and, since they are gauge invariant they must also be color singlets, 
i.e. singlets of the algebra l(2.124|l . Then, they must also be singlets of the algebra of the 
fermion currents ip^(x)T a ip(x). One can form a "baryon" (=color singlet) at a site by cither 
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leaving it unoccupied or putting J\fc fermions with antisymmetrized singlet wave-function 
by applying the creation operator 

S\x) = e Al ... AMc ^ Ai {x)...^ AMc {x) (2.128) 

{tAi,...Atf IS the usual antisymmetric tensor). 

Fermi statistics allows at most one singlet per site. Otherwise one can distribute the N /2 
singlets arbitrarily (N is the total number of lattice sites). Thus, there are Nl/((N/2)l) 2 
degenerate ground states with a typical state being 

\{px}x >= n s ^ x )\° > 

occupied x 

I label them by the eigenvalues p x of the local Fermion number operators 
p{x) = ^(x),^(x)} = ^(x)^{x)-Nc/2. 

All matrix elements in the vector space of degenerate vacua of the first order Hamiltonian 
vanish, 

<{p*}x|#iIK}*>=o 

so second order perturbation theory must be considered and the matrix with elements 
< {px}x\H 2 \{p' x } x > - < {p x } x \H 17 - —H 1 \{p' x } x > . 

must be diagonalized. The first term is zero and diagonalizing the second term is equivalent 
to diagonalizing the four-fermion Hamiltonian 

H sB = -K^^ A {x + l)i>B{x)^ B {x)i) A {x + I) = I<Y,P(x + i)p{x) (2.129) 

[x,i] [x,i] 

(up to an additive constant), with 

t 2 
K = 2Kr - — > , (2.130) 

in the space of pure fermion states where p(x) has eigenvalues 

Px = ±Nc/2 
and 

X 

This is the antiferromagnetic Ising model. 

Now I consider arbitrary Af. Then the density is AfcN/2 fermions per site and the 
"baryon" (= color singlet) creation operator is 

Sl... aMc (x)^e Al ... AM ^^)--^Z C c \x) . (2.131) 

Since this operator is symmetric in the lattice flavor indices ai, . . . , aw c , it carries an irre- 
ducible representation of SU{N) with Young Tableau with He columns and 1 row. Fermi 
statistics allows at most M singlets on a given site. Thus, the allowed representations of the 
flavor SU(Af) algebra at one site are the empty singlet and those with the Young tableaux 
given in fig. I|2.6(l and which are distinguished by the fermion numbers 

p x =N C {2v-M)/2, v = 0,1,2,..., N (2.132) 

24 



Mc 



N 



Figure 2.6: Young tableaux representations allowed at a given site. 

respectively. A ground state of Hg is obtained by creating AW/2 color singlets. The states 
are labeled by the local fermion densities p x and the vector in the corresponding SU(N) 
representation at each site and are degenerate. 

Again this degeneracy must be resolved by diagonalizing perturbations. The first-order 
perturbation to the vacuum energy vanishes and diagonalizing the second order perturba- 
tions is equivalent to diagonalizing the effective Hamiltonian 

H cS = -KJ2^a(x + iW B {x)^ {x)^ A {x + i) = KY,J afi {x + i)J 0a (x) (2.133) 



(up to an additive constant) where the operators 
1 



J^(x) = - [r A Hx), i> A (x)\ = <(*)<(*) - 8^ J -f = (J**(aO)t 
are generators of the Lie algebra of U(A/"), 



(2.134) 



a,3 (x), J Q,/3 '(y)] = S(x - y) ^ a 'j afi '{x) - S 13 ' a J a ' ^ {x)^j . (2.135) 



J" 



Thus, the effective Hamiltonian is that of a U(AT) quantum antiferromagnet with represen- 
tations given in fig. i|2.7[l and with the constraint ^2 x p x = 0, and I have shown that it 
is equivalent to the strong coupling limit of SU(Afc) lattice gauge theory with M lattice 
flavors of staggered fermions. 



The SU(2) color group is peculiar and this can be traced back to the fact that for 2 
colors the singlet creation and annihilation operators involve only 2 fermion field operators, 
and in second order perturbation theory terms are generated describing the hopping of these 
singlets |24|. Due to this Eg. I|2. 133(1 is not the right effective Hamiltonian but 



H t , 



-K 



E (^a(^b(x + i)KH* + *)V£(aO - £ac£bd^ a Hx + i)^(x)^(x + i)^ D (x] 



KJ2 (j ap (x)J Pa {x + I) + Sf{x)S f i a (x + l) 

[x,i] 



(up to a constant), with 



-,a/3 



S+{x) = CABlpA {XWB ( X ) > 



Q t^\„/,/3t/ 



(2.137) 
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S^(x)=e AB n(^ B (x) = (Si a (x))^ , (2.138) 

and the J a @(x) as above. One can check that the operators J al3 {x) and S±(x) obey the 
relations of the Lie algebra of the symplectic group Sp(2A^), i.e. Eqs. I|2.135|) and 

~Sf (x) , S c l p ' (y)] =6(x- y) (> Q ' J q/5 ' (a;) + S a0 ' J pa ' (x) + 5 aa ' J w ' (x) + 5^' J aa ' (»)) , 

(2.139) 
'j a ^x),S^ (3 '(y)]=S(x-y)^ a 'sf'(x)+6^'Sl a '(x)) , (2.140) 

J a0 {x),S°L p \y^ = -8{x - y) (8 aa ' S?/ {x) + 5 afy S ( l a ' [x]) , (2.141) 

hence H2.136|l is a Sp(2iVi) antiferromagnet. Especially for J\f = 1 it is the spin 1/2 Heisen- 
berg antiferromagnet as originally found in |24|. 



If, instead of SU(Afc) one had the gauge group XJ(Nc), gauge invariance would require 
that one imposes the extra constraint 

0°(aOI* >=0 , Vs. 

This can be fulfilled without electric fields, i.e. with p x — Vie, only when J\f is even. Then 
the allowed SU{Nc) representation at each site is the one with the Young Tableau given in 
fig. (|2.7(l . and, with the constraint 

E JQQ w = ° . 



J a/3 (x) generate the SU(Af) algebra in this representation. Thus, U{Nc) gauge theory is 
equivalent to an SU(Af) quantum antiferromagnet in the representation fig. (|2.7J) . This 
actually is true for all J\fc > 1 {■N'c — 2 is not special here as no hopping of color singlets is 
allowed for U(Nc) color groups [24jV 

-* A/c *- 
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Figure 2.7: Young tableau with Mc columns and Af/2 rows. 



As a concrete example, let us consider the case Af = 2 which gives 4 continuum flavors 
in d = 4. The SU(Nc < 3) gauge theory is equivalent to a U(2) antiferromagnet. Using the 
identity 

2<5 / 3 7 (5 Q , e = a a p ■ <7 7e + 5 a p8 1<L (2.142) 

(with a the Pauli matrices) one can change basis in the U(2) Lie algebra to find the effective 
strong coupling Hamiltonian 

HcS = tY1 (V (aO^Oc) -^(x + i)oip{x + i)+ p(x)p(x + ?)) (2.143) 

[x,i] 



2(> 



where 

(V 1 ^ = ^a a ^ A ) (2.144) 

and the representation at a given site has either p x — and the spin S = A/c/2 representation 
of 5(7(2) or ^ = ±iV c and the spin 5 = 0. 

If the gauge group is U(Nc) then the density is constrained, p x — and each site is 
occupied by the S = .A/c/2 representation of SU(2) and (|2. 143(1 is the Hamiltonian of the 
quantum Heisenberg antiferromagnet. 

The most subtle case is that of a U(Mc) lattice gauge theory with odd N ' ■ There I first 
solve the problem of minimizing the contribution of U(l) electric fields to the energy. This 
is done by solving the auxilliary problem of minimizing the energy function 



with the constraint 



where 



E 



2 5 2 

[a-Y'J 

E°[x,i}-E°[x-l,i\\ + p(x) -0 (2.145) 



P {x) = ^{xm{x)-McM/2 



If Afc'N is odd, the charge density p{x) has only half-odd-integer eigenvalues p x . Since 
they have no zero eigenvalue, the constraint Ea, 12.145|) cannot be satisfied unless the U(l) 
electric field is non-zero. This is also true of SU{Nc) color singlet states, since they always 
contain an integral multiple of Ac particles. Then the operator ip^ip has eigenvalue which is 
an integer multiple of J\fc and when J\f is odd, the charge density p{x) again has no nonzero 
eigenvalues. The configuration which minimizes the U(l) energy is the most symmetric one, 

e°\x ±i]=± — f-:ns£i«> 

with the accompanying eigenvalues of the charge density 

If, instead of U(l) electric fields, one had color non-singlets, the energy would be higher. 
In chapter 4 the one-flavor lattice Schwinger model is studied as a concrete example with 
7V C = 1 and TV = 1. 



2.4 Chiral anomaly on the lattice 



We shall see in section (|2.4.1t that under some mild assumption about the lattice action, 
there are always an equal number of left- and right-handed Weyl particles in the continuum 
limit of a lattice fermionic theory. In the continuum the axial or chiral charge Q5 is defined 
as the integral of the chiral density of the current j$(x) = ip(x)^ fi j5tp(x) 

Q 5 = [d 3 xf 5 (x)= fd 3 x^(x,t) l5 4,(x,t) (2.146) 



According to Ea. p. 146(1 Q5 is the number of right-handed fermions minus the number of 
left-handed fermions minus the number of right-handed antifcrmions plus the number of 
left-handed antifcrmions. In the case of massless fermions Q5 is conserved at the classical 
level, not only in the case of free theories, but also in the case of interacting theories like 
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QED or a- models |S]. At the quantum level the chiral symmetry is broken by the anomaly. 
For example in massless QED the divergence of the axial current is given by |26| 

«(*) = ^e^ pa F^(x)F?°(x) (2.147) 

There is no way to avoid the physical consequences of the anomaly; one might also define a 
conserved axial vector current, but it is not gauge invariant and the physical consequences 
are unchanged. The axial anomaly appears in perturbation theory as a consequence of the 
linear divergences in the triangle graph with an internal fermion loop, and either three axial 
vector current insertions or one axial vector and two vector current insertions. The one-loop 
quantum corrections to a chirally symmetric classical action are not chirally symmetric. It 
is impossible regularize a quantum field theory preserving both axial and gauge invariance. 
Since lattice gauge theories are manifestly gauge invariant, it is an interesting question how 
axial anomaly appears there. In general the fermion doublers do contribute to the anomaly. 
Since the doublers appear with opposite chiralities, one might conclude that the anomaly 
on the lattice is zero. However the question is more involved [31| . 

I am interested in analyzing how the anomaly manifests itself in the lattice Schwinger 
models with different flavor groups. I already illustrated that staggered fermions explicitly 
break the Ua{^) symmetry but possess a discrete chiral symmetry corresponding to trans- 
lations by one lattice spacing. As we shall see in chapters 4, 5 and 6 in the continuum 
the Ua{^) symmetry is broken by the anomaly either in the one-flavor or in the multiflavor 
Schwinger models. In chapter 4 we shall see that on the lattice in the one-flavor case the 
discrete chiral symmetry is broken spontaneously. At variance, in the two-flavor model we 
shall see in chapter 5 that the ground state is translationally invariant, i.e. does not break 
the discrete chiral symmetry. The only relic of the continuum axial anomaly is the non-zero 
VEV of the umklapp operator. 



2.4.1 The Nielsen-Ninomiya no-go theorem 

Nielsen and Ninomiya |27j demonstrated that fermion theories on a lattice have an equal 
number of species of left- and right-handed Weyl particles, i.e. neutrinos, in the continuum 
limit. The consequence of this behavior is a no-go theorem for putting theories of weak 
interactions on a lattice. Moreover it is not possible in strong interaction models to solve 
the species doubling problem on a lattice in a chirally invariant way. The appearence of 
equally many right- and left-handed species of Weyl particles is an unavoidable consequence 
of a lattice theory under some mild assumptions. 

The hypothesis of this theorem are locality of the theory and exact conservation of 
discrete valued quantum numbers and charges assumed to have a density defined from a 
finite region. 

If one is interested in the strong interactions and wants to study a Dirac particle rather 
than a Weyl particle, there is no unsurmountable problem because a Dirac particle can be 
thought of as a composite of two Weyl particles, a right-handed one and a left-handed one. 
However there is no way to construct chiral invariant lattice models for QCD. 

Let us consider the general class of lattice fermion theories for which the bilinear part of 
the action for the N-component complex fermion field ip(x) is of the form 



S F = 



: JdtJ2 ^(xWx) - JdtJ2 v^)# (£ - y)^(y) ( 2 - 148 ) 

x x,y 
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/ 4>i 

with H an Hamiltonian where interactions of quartic or higher degree in tp are neglected. In 
fact, one can define the number of species of Weyl fermions depending only on the bilinear 
part of the action, while the interactions do not change the dispersion relation but just cause 
scattering processes. The action Sf should fulfill the following three conditions 



• Locality of the interaction, i.e. H(x) — > fast enough as \x\ — > oo that its Fourier 
transform H (p) is a smooth function. Thus the eigenvalues u>i(p) (i — 1, . . . , N) are 
smooth except in the case that two eigenvalues tOi(p) and u>i+i(p) coincide 

• Translational invariance on the lattice 

• Hermiticity of the N x N matrix Hamiltonian H. 

The assumptions made for the charges Q (lepton number, for instance) are 

• Exact conservation of Q, even at the scales where the lattice cut-off is relevant 

• Q is locally defined Q — ^ s p{x). The charge density p(x) is a function of the field 
variables ip(y) related to y within a bounded distance from x. 



• Q is quantized 

• Q is bilinear in the fcrmion field ip(x). 



The weak interactions and the Standard Model exhibit right- and left-handed particles that 
do not have the same hypercharge, so that parity and charge conjugation are broken in weak 
interaction processes. This implies that if nature were indeed built on a lattice, it should be 
possible to realize Weyl fermions with different quantum numbers for left- and right-handed 
ones. 

One has thus to understand where is the problem: one must give up the idea that nature 
is based on a fundamental lattice cut-off or some assumption of the no-go theorem, or weak 
interactions phenomenology, i.e. the usual understanding of parity violation. 

Nielsen and Ninomiya gave two proofs of their no-go theorem, a very technical one based 
on algebraic topology and homotopy and a more intuitive topological proof. The proof is 
mainly supported by the fact that the momentum space of lattice theory is periodic, i.e. 
there is the Brillouin zone 

7T 7T , 

~-<Pi<- (2.149 

a a 

where the end points w/a and —ir/a are identified. From a topological point of view pi runs 
on a circle S , thus the momentum space makes up a hypertorus S 1 (g> S 1 <g> S . 

Let review the intuitive proof of the no-go theorem for a general class of 1 + 1 dimensional 
lattice fermion theories, whose action is the 1 + 1 dimensional version of the 3 + 1 dimensional 
case Ea. (12 .148(1 . One thus takes i/j to be a complex N-component field in order to keep 
generality. A generic 1 + 1 dimensional Hamiltonian has non degenerate energy levels, in 
fact, in order to have just two-level degeneracy u>i{p) — LOi + \{p) three parameters must 
be restricted in the Hamiltonian H(p), but there is only one p. However the wave packet 
velocity 

_ duj^ 

dp 
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P=PS (2.150) 



is non zero at the Fermi energy 



uii(pf) = 



(2.151) 



Low energy excitations are such that fermion states with p close to pf are excited and for 
particles relevant at low energies one has 



Uiip) = 






0((p- Pf ) 2 ) 



(2.152) 



The particle at the crossing point with u>i = line is called a right(left) mover when 

dui 



^| P=P/ >0 (<0) 



(2.153) 



Antiparticles and particles have the same velocity. Assuming locality u>i(p) is analytic in 
p except for degeneracy points and a smooth curve u>i (p) is defined in tOi — p space. The 
curve is closed since on a lattice p is defined modulo 27r/a and thus runs on a circle Si. 
Such curves must cross equally many times from Ui(p) < to u>i(p) > as they cross from 
u>i(p) > to u>i(p) < 0. One can orient the curve such that it is along the increasing p 
direction. The curve going up through u% = represents a right-mover and that going down 
a left-mover. From a topological point of view there must be the same number of upgoings 
and downgoings for the closed curve — see fig. (|2.8|) . Therefore, there appear equally many 




Figure 2.8: Dispersion relations for 1 + 1 dimensional lattice Weyl fermions. Each curve is 
closed since end points are identified. Each crossing point at u> — is a Fermi "surface" 
describing one species of Weyl fermions. 

right and left movers in the low energy regime. Let us introduce the shifted momentum 

P = p-P f (2.154) 

The dispersion relation reads 



diOi 
dp 



\p=0 P 



(2.155) 



This is a relativistically invariant dispersion relation for a massless mover along the right or 
left direction. This is true only when ^-\ p — Pf is considered equal to the velocity of light. 
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Let us now turn our attention to the 3 + 1 dimensional case. The eigenvalue equation 
for the Hamiltonian Ea. (|2.148J) reads 

Hffl1>(fl=Wifflil>(fi i = l,...,N (2.156) 

with N ordered eigenvalues u)\ > 101 > . . . > ujn- When two energy levels, the i-th and 
(i+l)-th coincide, i.e. are degenerate, at the momentum Pdeg, 

OJi{pdeg) = ^i+l(Pdeg) = (2.157) 

one can expand the N x N matrix H(p) near the degeneracy point pdeg- The most general 
Hamiltonian H^ 2 \p) for the two component spinor U(p) describing the i-th and (i+l)-th 
levels, is of the form 

H^ip) = UdegiPdeg) + {p ~ Pdeg)b + (p - Pdeg) fcV£ <T° + 0({p - p d egf) (2.158) 

with b and V constants. By introducing the shifted variables 

p=p-p de g UJ = U)-U>deg (2.159) 

H^(p) reads 

H^\p) = ~p k V*o a + ~p b (2.160) 

let us define further a new momentum 

Pa = u-{p-p deg )b (2.161) 

P a = ±~p k V* (2.162) 

The sign ± depends on the sign of detV which determines the relative handness of the p 
and p coordinates. Take conventionally p to be the right handed one. H^ 2 ' (p) becomes 

H {2) (p) = a ■ P (2.163) 

and one obtains the right and left handed Weyl equation 

PaU(p) = ±P U{p) (2.164) 

corresponding to the sign of detV. Wheter one has right- or left- handed Weyl fermions 
depends on the degeneracy point since the constant tensor V depends on pd eg ■ 

One may call the degeneracy point a right or left handed one according to Ea. f2.164j> . 
Each degeneracy point represents one species of Weyl fermions. If the Fermi energy surface 
lies at the degeneracy points of the i-th and (i+l)-th, each degeneracy point represents one 
species of Weyl particles. 

Let us consider curves in 4-dimensional to — p space or on the 3-dimensional dispersion 
relation surface defined by 

{&<*(&) | <o|w i (j)>=0} (2.165) 

where 

< a\u>i(p) >= aiV'f + a 2 ip ( 2 ) + ... + a N ip { ^ (2.166) 

The vector \a >= (01,02, ... ,ajv) is an arbitrarily chosen constant in the complex N- 
dimensional vector space, which may be chosen to be the basic vector corresponding to 
the field number one. Eq. (|2.166|l specifies a 1-dimcnsional curve in the generic case, in fact 
Ea. H2.166f) fixes two variables since < a\uji(p > is a complex number and it is a continuous 
and analytic complex-valued function of p. The set of curves always passes through all 
the degeneracy points, as can be shown explicitly using the continuum 2-componcnt Weyl 
equation near the degeneracy point pdeg- Then one can construct a new state 

piip) >= a\ui(p) > +0\wi+i(fl > (2.167) 
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near p — pdeg- One can always choose a and (5 such that 

< api(f) > = (2.168) 

which is the condition of Eg. ((2.166|) for the curve. Due to the existence of the Brillouin zone 
and the fact that the curves continuously pass through all the degeneracy points between 
the energy levels uii and Wj.fi, such a curve must be closed or consists of a number of closed 
curves. If the curves are oriented they must pass equally many times upward and downward 
between the energy levels u>i and u>i+i, through all the degeneracy points. Let us now orient 
the curves and show that, when the curve crosses the level to = between the levels uji 
and oJi+i upward (or downward), the crossing point, i.e. the degeneracy point is a right- (or 
left-) handed degeneracy point. The eigenstate determined by Eg. (|2. 156(1 is unique modulo 
a phase factor if one imposes the normalization condition < u)i(p)\u)i(p) >= 1. One can 
choose the phase of \uii{j>) > to be analytic in a simply connected region which should not 
include the degeneracy points. An orientation of the curves is assigned by means of the phase 
rotation of < a\uji(p) > on a small circle Si around the curve (|2.166|l . The convention is that 
an increase of phase on 5*1 should form a right-handed screw together with the oriented curve 
when one takes a right-handed coordinate convention for p. Let us consider for example the 
case the curve is along p z > on the upper cone pq > near the right handed degeneracy 
points. The two component field U satisfies the equation 

PaU{p) = P a U{p) (2.169) 

One draws a circle S\ with radius R around a point on the curve at a distance d from the 
degeneracy point and let us take R <§; d. The point Q on S\ is parametrized in the following 
way 

p x = Rcos9 (2.170) 

p y = R sin 6 (2.171) 

Pz = d (2.172) 

The normalized eigenvector U of Eq. (|2. 169(1 is solved in the limit |cl 

u = ( iz ie ) ■ (2,173) 



Since one obtains 

on the curve R = 0, the vector < a\ must be < a\ = (0, 1) so that on Si 



U=( I ) (2.174) 



< a | Wl >= A e *e (2.175) 

2p z 

and according to the orientation convention the curve is oriented in the positive p z direction 
because it forms a right handed screw. Computing the case of po > and p z < 0, po < 
and p z < 0, po < and p z > one obtains fig. I |2.9h .) for the right handed degeneracy point. 
The case of left handed degeneracy point is depicted in fig. I |2.9b ). As one can see in fig. 
((2.9(1 . the curve continues through the degeneracy points conserving the orientation. The 
orientation determined in this way becomes the same all along the curve. The important 
fact is that the curve must be oriented away from the degeneracy point on the level sheet 
with right-handed particles (i.e. with positive helicity), while the curve must be inward 
oriented on the level with left-handed particles. Finally the oriented curve is closed and 
thus must go equally many times upward and downward through all the degeneracy points 
^i(Pdeg) = ^>i+i{Pdeg) = between the energy levels 0Ji(p) and u>i+i(p). When the Dirac sea 
(wj+i level) is filled and the Fermi surface lies at the degeneracy point u> = 0, there appear 
equally many species of right- and left-handed Weyl fermions. This concludes the proof of 
this no-go theorem of Nielsen and Ninomiya. 
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Figure 2.9: Around a degeneracy point all the curves of the type < a\oJi(p) > are oriented 
away on one level and inward on the other one. 



The strong interaction lattice models for QCD invented by Wilson, Susskind and Banks 
and Casher break chiral invariance in order to avoid the fermion doubling. Drell, Wcinstein 
and Yankiclowicz were able to construct a chiral invariant model for lattice QCD but it 
breaks locality of the interaction. There is some proposal in the literature to avoid the 
no-go theorem, giving up some of its hypothesis, and to put the Standard Model on the 
lattice; but up to now no generally accepted scheme does exist. 

Nielsen and Ninomiya pointed also out a similarity between the fermion system of lattice 
gauge theories and the electron system of crystals. More precisely they pointed out that 
there exists an effect analogous to the mechanism of the Adler-Bell-Jackiw (ABJ) axial 
anomaly in solid state physics. The main similarity between lattice fermions and "real" 
electrons in crystals is that in both systems there is only the lattice translational invariance 
which implies the Brillouin zone existence. The crystal electrons are described by a one 
component non-relativistic Schrodingcr equation, but the energy eigenvalues form bands. It 
is possible to write the electron theory as a multicomponcnt lattice fermion theory with a 
matrix Hamiltonian. At points where two energy bands of the electrons make a contact in 
the energy-momentum dispersion law space, the two band wave functions of the electrons 
can be approximated by the Weyl equation, and so, according to the no-go theorem, such 
degeneracy points occur in pairs of the right-handed and left-handed type. There is no 
net production of electrons like the abscence of the net production of Weyl fermions in 
lattice fermions theories when parallel electric, E, and magnetic, H, fields are put on. This 
leads to the conservation of axial charge. The ABJ anomaly manifests itself by transferring 
electrons from the neighborhood of the left-handed degeneracy point to the right-handed 
one in energy-momentum space. In a parity non-invariant zero-gap semiconductor the effect 
analogous to the ABJ anomaly gives rise to a peculiar behavior of the conductivity of the 
electric current in the presence of the magnetic field. 

In chapter 5 I shall compute the VEV of the operator F — ip L tp L ip R ip R , describing an 
umklapp process in which two right-movers are annihilated and two left-movers are created. 
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I shall show that VEV of the umklapp operator is non-zero due to the coupling of left- and 
right-movers by the gauge field. Right-handed fcrmions are annihilated and left-handed are 
created at a degeneracy point so that there is no net production of electrons. 
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Chapter 3 



Antiferromagnetic spin chains 



The title of this chapter deserves an entire book to be properly developed, and in fact 
some books, reviews and long articles about spin chains do already exist |J2|. Spin chains 
interested mathematical physicists for their exact solvability, field theorists for the possibility 
to test their methods by comparation with exact results and condensed matter physicists for 
the possibility to understand the spectra of real physical systems. It is my purpose to report 
a self-contained analysis of antiferromagnetic spin chains, useful to understand chapter 5 and 
6 where I shall show that the multiflavor lattice Schwinger models in the strong coupling 
limit are effectively described by the spin-1/2 SU(J\f) antiferromagnetic Heisenberg chains. 

A common tool to study spin systems is spin-wave theory. Spin-wave theory was devel- 
oped by Anderson, Kubo and others [35] in the 1950s and it is a very powerful approach 
to quantum magnetic systems in dimensions greater than one, predicting long-range order 
and gapless Goldstone bosons. The situation remained clouded for magnetic chains be- 
cause it was known that long-range order could not occur in one dimension, but the Bcthe 
ansatz predicts massless excitations. There is "quasi-long-range-order" corresponding to a 
power law decay of spin-spin correlators and there are gapless excitations which are not true 
Goldstone bosons. 

Quantum spin chains have been extensively studied in the literature, starting from the 
seminal paper by Hans Bethe |3fi| in 1931 for the spin-1/2 case, where he introduced an 
ansatz to write down the eigenfunctions of the Heisenberg Hamiltonian, describing a chain 
with periodic boundary conditions. Seven years later Hulthen |37| was able to compute the 
ground state energy of the antiferromagnetic Heisenberg chain. Wc had to wait until 1984 to 
know exactly the spectrum of this model, when Faddeev and Takhtadzhyan |3()| analyzed the 
model using the algebraic Bethe ansatz and showed that the only one-particle excitation is 
a doublet of spin-1/2 quantum excitations with gapless dispersion relation. This excitation 
is a kink rather than an ordinary particle and is called spinon. All the eigenstates of the 
antiferromagnetic Heisenberg Hamiltonian contain an even number of kinks, nevertheless 
the kinks are localizable objects and one can consider their scattering. There are no bound 
states of kinks in this model. 

The one dimensional spin systems are not only unusual because they are exactly solvable 
but also for their incompatibility with long range order. Actually solvability and absence of 
long range order are deeply related concepts; systems with spontaneously broken symmetries 
are more difficult to describe and resist analytical solutions. All isotropic antiferromagnetic 
quantum spin chains with short range interactions exhibit quantum spin-liquid ground states 
— i.e. states with short range antiferromagnetic correlations and no order. Moreover these 
chains have strange quantum elementary excitations above these ground states, that are not 
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ordinary spin waves but are usually called spinous — neutral spin-1/2 kinks. 

As R. B. Laughlin 38 points out:"spinons are not "like" anything familiar to most of us, 
but are instead an important and beautiful instance of fractional quantization, the physical 
phenomenon in which particles carrying pieces of a fundamental quantum number, such as 
charge or spin, are created as a collective motion of many conventional particles obeying 
quantum mechanical laws. The fractional quantum number of the spinon is its spin. It is 
fractional because the particles out of which the magnetic states are constructed are spin 
flips, which carry spin 1." 

One significant result of my research [23 0Uj has been to show explicitly how spinons 
appear in lattice gauge theories. More precisely I showed that the massless excitations of 
the two-flavor Schwinger model coincide with the spinons of the spin-1/2 antiferromagnetic 
Heisenberg model. 



Section (|3.1|) shows the failure of spin wave theory in one dimension, due to the infrared 
divergencies. Spin-spin correlation functions are computed in the spin wave theory and are 
compared with exact expressions. 

In section H3.2[l I review the Haldane-Shastry model @T], a spin model discovered in- 
dependently by Haldane and Shastry and in the context of one dimensional Bose gas by 
Sutherland. The Haldane-Shastry model is a realization of the spin-1/2 Heisenberg chain in 
which the quantum disordered ground state and spinon excitations are particularly easy to 
understand. 

Section i|3.3[l is devoted to review the Bethe ansatz solution of the antiferromagnetic 
Heisenberg chain. Moreover I shall compare the exact solution given in jHO] with a study of 
finite size chains of 4,6, and 8 sites. An original result presented is the thermodynamic limit 
coefficient of the states with N — 2 domain walls composing the ground state. The spin-spin 
correlation function J2 X < 9- s -\S x ■ S x +2\g-s. > originally computed by M. Takahashi 42 , 
is derived in a slightly different way. The vacuum expectation value of the square of the 
vector V = ^ x S x A S x +i is computed in an approximate way. 

In section H3.4(l I introduce the spin-1/2 SU(J\f) antiferromagnetic Heisenberg chains. To 
provide an intuitive picture I study the SU(3) two sites chain and I find the ground state. 
A short review of the literature on SU(Af) antiferromagnetic Heisenberg chains is provided. 



3.1 Failure of spin wave theory in D=l 



In the limit of large spin S — > oo the quantum nature of a magnetic system can be forgotten, 
and one can consider the spin S as a classical vector in the three dimensional space. In fact 
in this limit the spin commutator is much smaller than the square of the spin variables 

[S a , S b ] = ie abc S c = O(S) « 0(S 2 ) (3.1) 

The classical ground states of the Heisenberg model are states with spins parallel for the 
ferromagnetic case 
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or neighbouring spins antiparallcl for the antiferromagnetic case 



numm... 

The classical ground state of the antiferromagnet is called Neel state. By introducing the 
raising and lowering spin operators 

X ~~ X ^ X V / 

the antiferromagnetic ( J > 0) Hcisenberg Hamiltonian reads 

N 

(3.3) 



f 1 

Hj = J 2_^ { i>x^x+l + 9 \ * "x+l + ^x ^x+1, 
x = \ ^ 



The classical ferromagnetic ground state coincides with the quantum one, but the Neel state 
is not an eigenstate of Hj due to the raising and lowering spin operators action. 

It is an interesting question to investigate whether or not long-range antiferromagnetic 
order occurs 

< S*S$ >= ±m (3.4) 

where the sign depends whether x belongs to the even or odd sublattice. The spin wave 
approach to this problem is to perturb away from the limit in which the Neel state is exact. 
This happens of course for a spin-S Heisenberg antiferromagnet in the limit S — > oo. In 
this limit one expects that the fluctuations of S 3 around +S on one sublattice and —S on the 
other sublattice are very small. It is customary to study the effects of these fluctuations using 
the Holstein-Primakov transformation, which consists in parametrizing the spin operators 
in terms of boson operators a and a^ on the sublattice of even sites A and b and W on the 
sublattice of odd sites B. On the sublattice A one writes 

S 3 = S-a^a (3.5) 

S~ = (2S)*ot(l-^)* (3.6) 

Eas. H3.5l3.6(l for spin operators correctly reproduce the spin algebra. The state with maximal 
S 3 = S is a state with zero bosons. Applying S~ on this state one creates a boson and so 
lowers S 3 by one. The state with S 3 = —S has 2S bosons and is annihilated by S~~ . On 
sublattice B one writes 

S 3 = -S + tfb (3.7) 

S~ = (2S)Hl- b ^)h (3.8) 

Using the Holstein-Primakov representation given by Eas. l|3.5l3.()|) and Eas. l|3.7l3.8|) and 
studying the small fluctuations around the Neel state so that a) a and tfb are O(l) <§C S in 
the large S limit, one may develop an expansion in i by expanding the square roots in the 
definition of S~ in Eas. H3.6l3.8(l . To quadratic order the Hamiltonian H3.3(l reads 

JV 

H=JJ2{S 2 + S(ala x + bl +1 b x+1 + a x b x+1 + 4&l+i)} (3-9) 

x=l 

where higher orders in -4 have been dropped. In momentum space the Hamiltonian Ij3.9|l 
reads 

Hj = 2JSJ2 {lk{a k b- k + a{bl k ) + (a{a k + b{b k )} (3.10) 



k 
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where 

Ik = y = cos/c = 7_ fc . (3-11) 

The Hamiltonian (|3.10|1 describes free bosons. Let us take the Bogoliubov transformation 

C k = U k a k - V k bl k D k = U k b k - V k al k (3.12) 

where, since a and b are bosonic operators, one has 

\U\ 2 -\V\ 2 = l . (3.13) 

To put to zero the off diagonal terms inverting Eas, 13.12(l one has 

7 (C/ 2 + V 2 ) + 2UV = (3.14) 

Solving Ea. l|3.13|l and Eq. (|3.14|l one gets 

1 1 



U k = 
V k = 



i i 



V 1 - Ik + Vl+7fc 



y/1 - 7fe - y/T+ 



Ik 



2 (l-7, 2 ) 
Using Eas. l|3.12|l the Hamiltonian Ea. l|H.l()|l reads 

H.j = 2JS£(1 - j 2 )HctC k + D\D k ) 



(3.15) 
(3.16) 

(3.17) 



The free excitations created by C and D are known as spin waves and correspond to the 
infinitesimal deviations of the spins away from the Neel state. The dispersion relation for 
spin waves is 

E{k)=2JS{l- 1 l) 1 * (3.18) 

that vanishes linearly at k = 

E{k)^v\k\ (3.19) 

where the effective "light velocity" is 



2JS 



(3.20) 



In a field theoretical language spin waves are gapless due to the Goldstone theorem. In fact 
there are two Goldstone modes C and D corresponding to the breaking of 5*0(3) down to 
the SO(2) rotations about the z axis. Low energy long-wavelength spin waves correspond to 
configurations where all regions are locally in some Neel ground state but the direction of the 
sublattice magnetization vector makes long- wavelength rotations. The two different modes 
C and D have spin S 3 = ±1 and correspond to raising S 3 on one sublattice or lowering it 
on the other. 

Now one can compute the reduction in the sublattice magnetization due to quantum 
fluctuations 



<s . w>=s ./g 



dkl 

2 



1 



1 



r . (3.21) 

V 1 -7| 

Quantum spins like quantum harmonic oscillators have zero point motion. In one dimension 
one finds that Ea. l|3.21|l is divergent for small wavelengths 



A. .s"(,i - = -*'•,,•, .-S = - /gl = -oo 



(3.22) 



Ea. (|3.22|l indicates that the Neel state is destabilized by quantum fluctuations in one di- 
mension even in the large S limit. The divergence of Eq. l|3.22|l is in agreement with the 
Coleman and Mermin- Wagner theorems |43j preventing the spontaneous breaking of con- 
tinuous symmetries in (l+l)-dimensions due to infrared divergences connected with the 
Goldstone bosons. 
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3.1.1 Spin-spin correlators in spin wave theory 



We have seen that the Neel state is destabilized by quantum fluctuations in D=l even in 
the classical limit S — > oo. The spin— spin correlation function to next order in i is 



< S Q ■ S r >« ±S 2 



l-^log- + 0(i) 

irb a S z 



(3.23) 



Eq. (|3.23(l indicates that at large S the system is ordered at small distances while the disorder 
occurs at exponentially large scales 

£ ~ e nS (3.24) 

More precisely, one should distinguish the r odd and even cases. To sec this, it is instructive 
to compute the finite distance correlators from spin wave theory. If r belongs to the odd 
sublattice, the spin-spin correlator reads 



< OO ' Or ^ — — & 

= -S 2 
while if r belongs to the even sublattice one has 



1 1 (1 l^l_cosrkc^ ) + 0{s _ 2 , 



S y N 
1 



| sinfcj 



l + -(l-/ (r))+0(5- 2 ) 



(3.25) 



< <Sq • &r 5> — — J 



1 



S { N ^ 



1 v-^ 1 — cos rk 



I sin k I 



) + o{s- 2 ) 



-S 2 



1 + -(1-I e (r)) + 0(S- 2 ) 



In Ens.J3.25l3.2fil) 



Io{r) 



Ur) = 



2Wo 



1 — cos rk cos k 



2tt 



2;r 



| sinfc| 
1 — cos rk 
I sinfcl 



-dk 



dk 



(3.26) 

(3.27) 
(3.28) 



and up to the distance r — 10 one gets the values given in table (|3.1|l The spin-spin correlator 



Table 3.1: I (r) and J e (r) 



Io{r) 


h(r) 


'o(l) = 1 


W) = 1 


^»(3) = g 


m - 1 


'.(5) = & 


'.(«) = If 


w = m 


7 «( 8 ) = nib 


m - m 


^e(lO) = ff 



G(r) =< 5*0 ■ j r > are given in table i|3.2|) . It is interesting to compare table i|3.2|) with the 
numerical correlators given in table (|3.7|l . While it is quite impressive the agreement of G(r) 
for r = 1,2 with the exact results, increasing the distance r the agreement becomes very 
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Table 3.2: Spin- wave theory spin-spin correlation functions 



r 


G(r) 


1 


-0.4317 


2 


0.1817 


3 


-0.0073 


4 


0.0756 


5 


0.1625 


6 


0.0119 


7 


0.2716 


8 


-0.0335 


9 


0.3525 


10 


-0.0689 



poor already for r = 3, and starts to have no meaning at r = 5, where one gets even a wrong 
sign, opposite to the right answer. I think that this is and interesting lesson that one can 
learn from the failure of spin-wave theory. While, in studying short distance correlations 
(r = 1,2) one gets very good results in the spin- wave theory, quantum disorder prevails 
for larger distances and so a Neel state approximation cannot pick up the main features of 
the system. The expected scale for the appearence of quantum disorder for S = % 1S > from 
Ea.l|3C2iJl. f; = ef w 5. 



3.2 The Haldane-Shastry spin model 



The Haldane-Shastry spin model [H] is a generalized spin-1/2 Heisenberg Hamiltonian with 
long range inverse-square exchange interaction between spins. The inverse-square interaction 
makes the model more easy to solve than the original Heisenberg model and the description 
of the excitation spectrum in terms of spin-1/2 spinons is particularly easy to understand. 
The model can be considered as a discretized version of the Calogero-Sutherland model |4*4*j 
describing a one-dimensional Bose gas with inverse square repulsions. 

This paragraph shortly reviews the approach to the Haldane-Shastry model given by 
Laughlin et al. in [J2|. My purpose is to stress some feature of one-dimensional antifer- 
romagnetic chains that are more intuitive to understand in this model than in the Heisen- 
berg model. The Haldane-Shastry Hamiltonian describes N spin-1/2 particles distributed 
equidistantly on a circle, coupled with an exchange interaction that decays proportional to 
the inverse square of the chord between the spins 



TV *--' \Z a — Zg\ 



(3.29) 



where each site is parametrized by the roots of the unity 



^-1 = 



(3.30) 



and S a is the spin operator on site a. The Hamiltonian Eg. (|3. 29(1 is invariant, due to periodic 
boundary conditions, under translations generated by T 



T 



'J a 



(3.31) 
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and under global rotations in the spin space generated by 

N 

S = Y,S X ■ (3.32) 

x=l 

Of course, 



S,T 



(3.33) 



Moreover the Hamiltonian Eq. 13.29(1 is invariant under the special internal Yangian sym- 
metry |35] 

H HS ,K\=Q (3.34) 

where 

A = |E(? ± ?)(^ A ^) (3-35) 

The ground state of this model has the same functional form as the fractional quantum Hall 
ground state | 46II47| : it is a spin singlet in the case of N even and exhibits quantum disorder: 
i.e. it has the same relation to the antiferromagnetically ordered Neel state that a quantum 
liquid has to a conventional crystal. 

In the hard-core bosons picture, the spin f on a site corresponds to the presence of 
a boson, the spin J, to the abscence of a boson. The ground state wavefunction is the 
complex coefficient assigned to each configuration with N/2 bosons. The ground state of 
the Hamiltonian 13.29[) reads 

JV JV 

2 2 

ipg. s . (zi, . . . , zn ) = Y[ {z a - zp) 2 Y\_ z °< (3.36) 



a</3 a=l 



and its energy is 



E 9-s.=-J(^)(N+^) (3.37) 

where z\, . . . ,z» denote the localization of the N/2 bosons. It is easy to verify that the 
ground state wave function Eq. H3.36[) is real, and it is translationally invariant with momen- 
tum or 7r for N/2 even or odd. In fact the translated expression of one site of tpg.s. is 
obtained by multiplying each of its arguments by z — e % ~^~ and due to the fact tpg. s . is an 
homogeneous polynomial of degree (N/2) 2 one gets 

■ NTT 

i/) g . s .(ziz,.. -,znz) = e I ^~ip g . s .(z ll .. -,zn) (3.38) 

One can immediately verify that ipg.s. is a spin singlet. Due to the fact that Zj is constrained 
to take N/2 possible values, i.e. that there are N/2 spins up (and N/2 down), S z ip g . s . = 0. 
Moreover it is the lowest weight state 

S-^ g . s . = (3.39) 

It is worth stressing that it is very peculiar of this model the possibility to write down 
in a closed form, so simple to read, the exact quantum ground state (|3.36|l . In general, 
it is very difficult to write down in an easy compact form the ground state of a quantum 
antifcrromagnetic chain. In the case of the Hciscnbcrg antiferromagnet this is not possible 
for a general N. 

ipg.s. exhibits quantum disorder, i.e. the spin-spin correlation function 

< b a - bp >= — 7- (3.40) 

< Wg.s.Wg.s. > 
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goes to zero as | z a — zp | — > oo so that it is a spin-liquid state. However the fall-off of Eq. (|3.4()j) 
is slow and this is typical of half-integer spin systems that exhibit gapless excitations 48 . 
Different is the case of integer spin systems that are strongly-disordered spin liquids with 
exponential decay of correlations on a length scale £ and have an energy gap for excitations 
A = ?y, where v is the spin- wave velocity of a nearby ordered state. 

The ground state Eq. i|3.36[) is not degenerate, but the state 

JV 
2 

1 \{zl (3.41) 

ct</3=l a 

is almost degenerate with ipg.s. 



^' gs ( Zl ,...,ZN ) = Y[ { Z a ~Zfif 
a<0=l 



H H sW g . s .>=-J{^){N - -)W g . 8 > 



(3.42) 



and it has momentum 7r greater than ip g . s .- The state ?// s can be thought as ip g . s . plus a 
pair of spinons excited out in a singlet state of momentum ir. When the number of sites N 
is odd, the ground state cannot be a spin singlet, but it is a spin-1/2 state, i.e. it describes 
a spinon of momentum ±7r/2. The wave function 



ip a (zi,---,ZN-i) = JJ (z a - z 7 ) Y[ ( z 7 ~ z s) 2 Yi 



(3.43) 



7=1 



7<<5=1 



7=1 



describes a J, spinon at site a and it can be interpreted as a J. spin on site a surrounded by 
an otherwise featureless singlet sea. The linear combination of the states given by 



N 



4>m{Zl, 



,Z M 



) = 52(z* a ) m Mzi,...,ZM) 



(3.44) 



with < m < „ 1 , is an exact eigenstate of Hhs with eigenvalue 



1, J,2ti\ 2 ,JV-1 



- m) \ \tp m > 



lluslr,,, -=<j-J(_)(n--) + -(-)^m( 2 



The state \ip m > represents a J. spinon travelling with crystal momentum 

■k 2ir 1 . 

q = -N-—(m + -) [mod 2tt) 



due to the definition 
Rewriting Eg. 13.45(1 as 



ip m (ziz, . . . , zmz) = e iq ip(zi, . . . , zm) 



U »«\ r " -\~ J( ^) {N+ ^^ ) + Eq ^ r "' ' 



one obtains the dispersion relation 



E «=2 



'2 



(o) 2 -? 2 



(mod 7r) 



(3.45) 

(3.46) 

(3.47) 

(3.48) 
(3.49) 



The allowed momenta for the spinon are only half of the Brillouin zone (the inner — ? < q < 
5 or the outer part — ir < q < — ? and ^ < q < ir, depending on wheter N — 1 is divisible 
by 4). The spinon dispersion relation is linear in q close to ±5 with velocity 



7T J2?r 
2 ft AT 



(3.50) 
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Spinons are "relativistic" excitations travelling with velocity v sp i non and they exist because 
quantum disorder exists whenever the total spin in the unit cell is half integral. A spin 
liquid singlet ground state can exist only if the total number of spins is even. If N is odd the 
total spin cannot be less than 1/2. Since in the thermodynamic N — > oo there can be no 
difference between N even and odd, the spin system must have already had a neutral spin-i 
spinon excitation at the beginning, i.e. Ea. lj3.47p describes a spinon. In the case of N odd, 
the ground state is 4-times degenerate and is given by Ea. (|3.47p for m — and ^p- and 
their | counterparts. Even if more than one spinon is present in the chain, they keep their 
identity. In the thermodynamic limit spinons are non-interacting excitations, some kind of 
quantum kink-excitations. Spinons are neither bosons nor fermions but they are semions, 
i.e. particles obeying 1/2 statistics. The two-spinon wave function is given by 

i>AB = \\i z j - z a)(z 3 - Z B ) \\_{ z j ~ z k) 2 (3.51) 

3 3<k 

where za and zb are not necessarly lattice sites, and it has the property 

1 N 

ip AB {zi,. .. ,zn_ 1 ) = (z A z B ) ~~iPab( z i,-- ■,zjv_ 1 ) . (3.52) 



The phase for adiabatic motion of spinon A in the presence of spinon B, called Berry 
phase, is 



1 



i ^ . d , . , 

< i>AB\ z A-R WAB > + < Z A^\ VABWAB > 

dz A OZA 



1 N 

jCl-y) (3-53) 



2 < iPab\^ab > 
By moving spinon A around the loop to exchange it with B, one gets the phase 

A$s=/i(l-v)— = ± J* ( mod 2 *) ( 3 - 54 ) 



One would get in the case of bosons or fermions 

A$ B = 

A$ F = n (3.55) 

Fractional statistics has nothing to do with the symmetry properties of ipAB under the 
exchange of A with B, but it has to do with state-counting. 

Last I want to say that the the Yangian symmetry operator Eq. (|3.35p represents a scaled 
spin current. Acting on the propagating spinon Ea. H3.44p one gets 

A Z |V-m >= \^P- - m \ IV'm > (3.56) 

So the eigenvalue of A z is proportional to the spinon velocity 

dE q _2*JIN-l ^ 



dq N { 4 

The action of A z on a multispinon state is the sum of all the spinon velocities. 

I analysed the model at half filling, i.e. when every site is occupied by one particle. In 
the presence of doping the possible configurations on each site are spin up, spin down and 
empty site. When the filling is less than half, new quantum excitations called holons are 
created. Holons are charged quantum excitations of spin 0. New phenomena appear in 
presence of holons and it is necessary to add to Hhs & term describing the holon dynamics. 
Since my thesis is concerned with the relationships between spin chains and field theories, 
the half filling constraint is required. The reader interested in knowing more about holons 
can look up Ref. |49| . 
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3.3 The Bethe Ansatz solution of the antiferromagnetic 
Heisenberg chain 



In this paragraph the Bethe ansatz solution of the antiferromagnetic Heisenberg chain |3U| 
is discussed in detail. Moreover I shall compare the exact solution given by L.D. Faddeev 
and L.A. Takhtadzhyan in [201 with a study of finite size chains of 4, 6 and 8 sites. I show 
that already these very small finite systems exhibit spectra that match very well with the 
thermodynamic limit solution. I suggest to the reader interested in the subject the references 

m 

The Bethe ansatz is a method of solution of a number of quantum field theory and sta- 
tistical mechanics models in two space-time dimensions. This method was first suggested 
by Bethe jSE] in 1931 from which takes its name. Historically one can call this formulation 
the Coordinate Bethe ansatz to distinguish it from the modern formulation known as Al- 
gebraic Bethe ansatz. The eigenf unctions of some (l+l)-dimensional Hamiltonians can be 
constructed imposing periodic boundary conditions which lead to a system of equations for 
the permitted values of momenta. These are known as the Bethe equations which are also 
useful in the thermodynamic limit. The energy of the ground state may be calculated in 
this limit and its excitations can be investigated. 

I review here the method applyied to the study of the antiferromagnetic Heisenberg 
chain |3()j . In particular it is shown that there is only one excitation with spin-1/2 which 
is a kink: physical states have only an even number of kinks, therefore they always have 
an integer spin. The one dimensional isotropic Heisenberg model describes a system of TV 
interacting spin-^ particles. The Hamiltonian of the model is 

N 1 

Hj = jJ^iSx ■ S x+1 - -) . (3.58) 

x = l 

where J > ( J < would describe a ferromagnet) and the spin operators have the following 
form 

S x = li <g) 1 2 «> . . . <g> y®---® 1 w ■ ( 3 - 59 ) 

They act nontrivially only on the Hilbert space of the x th site. Periodic boundary conditions 
are assumed. 

The Hamiltonian (|3.58|) is invariant under global rotations in the spin space, generated 

by 

N 

S = Y,S X ■ (3.60) 

x=l 

Due to the periodic boundary conditions, under translations generated by the operator T, 

f :S X -^ 5,_! (3.61) 

the Hamiltonian is invariant and [S^T] = 0. 

In order to diagonalize Hj it is convenient to use an eigenfunction basis of operators 
commuting with Hj, so obviously S 2 , S z and also T. Let us sketch the Coordinate Bethe 
ansatz technique. One has to introduce the "false vacuum" 

|fi>=]7®IT>* (3.62) 
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s+\ T>* 


= 


S*\ T>x 


- slt>» 


*a; = 


z Oj, zt O™ 


^> = 


Tl«> 



with 

(3.63) 
(3.64) 

where 

(3.65) 

(3.66) 

5?2| " > = f ( f + m> (3 ' 67) 

f\Q > = \Cl > . (3.68) 

All the other basis vectors have S 3 < ^ and one can get them by properly acting on |fi > 
with the lowering operators S~ . Let us start from the generic state with S 3 = % — 1 , with 
N — 1 spins up and M = 1 spin down 

N 

\M = 1 >=^2<f) x \x> with \x>=S~\n> (3.69) 

x=l 

where the coefficients <j> x must be such that \M — 1 > is a translationally and rotationally 
invariant state 

N N N 

f\M = l>=^<t> x \x-l>= y £ j <l> x+ i\x>=n y £,<t>x\x> (3-70) 

x—1 x—1 x—1 

and from Eci. (|3.70J) one gets 

(j) x + l = /10a; (3-71) 

&,+! = m x «^i x^N (3.72) 

0! = M^iv - n N <h. — ■♦ M W = 1 ■ (3.73) 



Setting 0i = l one has <f> x — fi x l . There are N possible values for (j> x , corresponding to the 
N roots of the unity. One of these roots corresponds to a state with S = ^ 

N 

^2\x>=S~\n> — ► /i=l (3.74) 

x = l 

while the other N — 1 roots have S = ^ — 1 . 

The generic case with M spins down is more complicated. Let us consider the case M = 2 
to understand what happens 

N 

\M = 2>= Y^ <f>{xi,x 2 )\xi,X2 > with \xi,x 2 >=S- i S- 2 \n> . (3.75) 

x±<x 2 = l 

By requiring the translational invariance of the state \M = 2 > one has 

f|M = 2 >= T\M = 2 > (3.76) 

and for x 2 < N one has 

<j>(xi + l,x 2 + l)=T < f>(x 1 ,x 2 ) (3.77) 

that would easily give 

4>{x 1 ,x 2 ) = ^'"V^ 2 ^ 1 ) T = ^i^ 2 (3.78) 
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but due to the periodic boundary conditions one has to find coefficients <j){x\, X 2 ) that satisfy 
not only Ea. (|3.77j) but also 

<f>(l,x+l)=T<f>(x,N) . (3.79) 

Eg. (|3. 79(1 is no more satisfied by i|3.78[l . Bcthe proposed the following ansatz for the coeffi- 
cients <t>{x\, X 2 ) 

4>{xi,X2) = A l>2 ^- 1 i%- x + A 2 , 1 ^- 1 iiT- 1 . (3.80) 

Eq. (|3.77I3.79|I are satisfied if the following equations hold 

A 12 = A 21 ^ (3.81) 

Aii = A 12 ^ . (3.82) 

By imposing to \M = 2 > to be an highest weight state 

S + \M = 2>=0 (3.83) 

taking into account Eq. H3.81I3.82J) and introducing the following change of variables 

Ma = ^f (3.84) 

A Q 2 

one gets the so called Bethe ansatz equations 

A« + f ^ A Q - A/3 + « 



In the general case of M spins flipped one gets 



\ - i M \ _ \ 

(^f=-n ; 7 * . (3.86) 



The energy and the momentum of a given state with M spins down can be expressed in 
terms of the parameters X a 



E 



M 



M T M 



2 ^-^ A 2 

a=l q=1 a ' 4 



M M \ _ i 

P M = ilnT=Y,P a =iJ2 ln T r • (3 - 88) 

a=l a=l ■! 

Energy and momentum are thus additive as if there were M independent particles and the 
A Q must satisfy the Bethe ansatz equations l|3.86|) in order for Em and Pm to be eigenvalues 
of the Hamiltonian and momentum operators. 

The solution of the antiferromagnetic Heisenbcrg chain is reduced to the solution of the 
system of the M algebraic equations (|3.86(l . This, in general, is not an easy task. It can be 
shown 30 , however, that, in the thermodynamic limit N —* 00, the complex parameters A 
have the form 

Xi=\j,L+H , l = -L,-L + l,...,L-l,L; (3.89) 

where L is a non-negative integer or half-integer, A^l is the real part of the solution of 
(|3.86|) and I shall define shortly the set of allowed values for the integer index j. The A's 
that, for a given \j,l, are obtained varying I between [— L, L] by integer steps, form a string 
of length 2L + 1, see fig. (|3.1(l . This arrangement of A's in the complex plane is called the 
"string hypothesis" (SOj. In the following I shall verify that, even on finite size systems, the 
"string hypothesis" is very well fulfilled. 
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Figure 3.1: Strings for L = 0, §, 1, 



In a generic Bethe state with M spins down, there are M solutions to (|3.86|) . which can 
be grouped according to the length of their strings. Let us denote by vl the number of 
strings of length 2L + 1, L = 0, |, . . .; strings of the same length are obtained by changing 
the real parts, Aj,l, of the A's in l|3.89[) : as a consequence j — 1, . . . , vl- If one denotes 
the total number of strings by q one has 



9 = $>z , M=^(2L+1H 



(3.90) 



The set of integers (M, q, {vl}) constrained by 1|3.90[) . characterizes Bethe states up 
to the fixing of the q numbers Xj.l] this set is called the "configuration". Varying M, q 
and vl, one is able to construct all the 2^ eigenstates of an Heisenberg antiferromagnetic 
chain of N sites [HO]. The energy and momentum of the Bethe's state, corresponding to 
a given configuration — within exponential accuracy as N — > oo — consist of q summands 
representing the energy and momentum of separate strings. For the parameters Aj,£ of 
the given configuration, taking the logarithm of (|3.86() the following system of equations is 
obtained in the thermodynamic limit 



where 



2JVarctan j ' L \ = 2ttQ jm + Y] Y] $l x l 2 (A J>£l - \ k ,L 2 ) 

Li + ^ ' 



L 2 fc=l 



Ll+L2 A » 

< &l 1 l 2 (A) = 2 yj (arctan— + arctan- -) 

Li Lj -\- 1 

Z,=|L 1 -L 2 |^0 



(3.91) 



(3.92) 



Integer and half integer numbers Qj,l parametrize the branches of the arcotangents and, 
consequently, the possible solutions of the system of Eas. l3.91|l . In ref.(H0j it was shown that 
the Qj,l are limited as 



with Q™ ax given by 



}T X < Qi,l < Q2,l <■■■< Q VL ,L < QT 



Q 



max 
L 



(3.93) 
(3.94) 
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and 

2min(i 1 ,£ 2 ) + l if P\ ^ L 2 



J ^ L ^ = {2 Ll li ' if l\ = U . ^ 

The admissible values for the numbers Qj.l are called the "vacancies" and their number for 
every L is denoted by Pl 

Pl = 2Q r [ iax + 1 . (3.96) 

The main hypothesis formulated in J3D] is that to every admissible collection of Qj,l 
there corresponds a unique solution of the system of equations l|3.91[) . The solution always 
provides, in a multiplet, the state with the highest value of the third spin component S 3 . 

Let us now consider some simple example. The simplest configuration has only strings 
of length 1, i.e. all the A's are real. The singlet associated to this configuration 

N 
M = q = v Q = — , v L =0 , L > , (3.97) 

is the ground state. The vacancies of the strings of length 1, i.e. the admissible values of 
Qj t 0, due to eas. (|3. 9313. 9413. 95j) . belong to the segment 



N 1 „ TV 1 
h - < Q, ■ o < 



- + t;<Qj,o<--t; ■ (3-98) 



Therefore they are N/2. All these vacancies must then be used to find the N/2 strings of 
length 1. As a consequence this state is uniquely specified and no degeneracy is possible. 

Next I consider the configuration that provides a singlet with 1 string of length 2 and 
all the others of length 1: 

N N N 1 

M = - , q = - - 1 , u Q = --2 , v,_ = 1 , v L = , L > - . (3.99) 

For the strings of length 1 the number of vacancies is again N/2; for the string of length 2 
there is one vacancy and the only admissible Qj,i equals 0. Thus, since the number of strings 
of length 1 is fo = -f- — 2, there are two vacancies for which Eas. H3.91(l have no solution; they 

are called "holes" and are denoted Q\ ' and Q 2 ■ This configuration is determined by two 
parameters: the positions of two "holes" which vary independently in the interval l|3.98() . 

There is another state with only 2 holes: the triplet corresponding to the configuration 

N 
M = q = v = —-1 , v L = , L > (3.100) 

The number of vacancies for the strings of length 1 equals ^ + 1 , while vq = ^ — 1. 

The excitations determined by the configurations H3.99I3.100JI belong to the configuration 
class called in 30 Maf- The class M.af is characterized as follows: the number of strings 
of length 1 in each configuration belonging to this class differs by a finite quantity from N/2, 
z/q = y ~ ko where fco is a positive finite constant, so that the number of strings of length 
greater than 1 is finite. From l|3.96|l one then has 

(3.101) 
__ L > (3.102) 

L'>0 

so that 

V 

(3.103) 



Po 


= 


— + fc - 2 ^ v L 

L>0 






Pl 


= 


2k - 


-2J2 J(L,L' 

L'>0 


)VL> , 






Po 


N 
> — 
~ 2 


, Pl < 2fc , 
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From (|3. 101|1 follows that the number of holes for the strings of length 1 is always even and 
equals 2 only for the singlet and the triplet excitations discussed above. One can imagine the 
class M-af as a "sea" of strings of length 1 with a finite number of strings of length greater 
than 1 immersed into it. It was proven in |30| that, in the thermodynamic limit, the states 
belonging to M-af have finite energy and momentum with respect to the antiferromagnetic 
vacuum, whereas each of the states which corresponds to a configuration not included in the 
class M-at has an infinite energy relative to the antiferromagnetic vacuum. 

Let us now sketch the computation of the thermodynamic limit ground state energy. 
Eas. 1(3.91(1 for the ground state have the form 

ttQ 1 N/2 
arctan2Aj = — — - + — >J arctan(Aj — Xk) . (3.104) 

fc=i 
Taking the thermodynamic limit N — > oo, one has 

Qr^X , -\<x<\ , \j^X(x) , (3.105) 

and Eas. 1(3. 104(1 can be rewritten in the form 

arctan2A(x) = irx + / arctan(A(a;) — Xvyyjdy . (3.106) 

Upon introducing the density of the numbers X(x) in the interval dX 



PW = ^ ( 3 - 107 ) 

dx \x=x(X) 

and differentiating Eas. ((3. 106(1 . one gets 

i r 00 - ~lll i 

PW = 7^ — ^e- jA|?l d£=T 'z-^r • (3.108) 

Hy ' 2tt J_ 00 1 + e-IH 2cosh7rA v ; 

The density p(A) introduced in this way is normalized to 1/2. It is now easy to compute the 
energy and the momentum of the ground state 

N_ 

2 nOO TAT ^.OO 1 

E g .,. = J2 e a = N / e(X)p{X)dX = -— / dX — = -JiVln2 (3.109) 

^ J-oo 4 J-oo (A 2 + |)coshvrA 

K 

2 /.OO IU roe AT 

P 9 .s.^ Pa = N j_j{X )p{ X)dX^-- j_JX^^- x = -, (mod 2.) . 

(3.110) 

According to Ea. l|3. 110(1 . P g . s . — (mod27r) for y even, and P g . s . = n (mod 2tr) for y odd. 
The ground state, as expected, is a singlet, in fact the spin S is given by 

N ^ N f°° 

<?= £l= JV/ P (X)dX = . (3.111) 

Let us analyze the triplet described by Eg. 1(3. 100(1 : Eas. ((3.91(1 take the form 

-l 

N ' N 



7T0- 1 

arctan2Aj = J + — 2_. arctan(Aj — Xk) (3.112) 



fc=i 



40 



where now the numbers Qj lie in the segment [— f, -j] and have two holes, Q{ and Q 2 
with Q^ < Q 2 ■ Taking the thermodynamic limit one gets 

o {h) o {h) O 1 

^-*i , ^_-z 2 , ^r^x + -(9(x-x 1 ) + 6(x-x 2 )) (3.113) 

where 9(x) is the Heaviside function. Eas. (|3.112[l become 

arctan2A(x) = irx + — (9(x — x%) + 8{x — x 2 )) + / arctan(A(x) — X(y))dy . (3.114) 



Ea. (|3.114|l gives, for this triplet, the density of A, p(X) = 



dx 



Pt {\) = p(X) + l(a(A - Ai) - a(X - A 2 )) (3.115) 

where p(X) is given in Ea. 1(3. 108(1 and 

1 f°° 1 

°w = ->i;Lu^® e * ■ (3 - 116) 

Ai and A2 are the parameters of the holes, A; = X(xi), i — 1,2. The energy and the 
momentum of this state measured from the ground state are now easily computed 



e T (Ai,A 2 )=iV / e(A)(p t (A)-p(A))dA = e(A 1 ) + e(A 2 ) (3.117) 

{• OO 

Pt(Ai,A 2 ) = 7V / p(\)(pt(\)-p(\))d\=p(\i)+p(\2) (mod 2tt) (3.118) 



where 



e(A)=/ e( MM A-^M=J^^ X (3.H9) 



/•°° 7T 

P(^) = / p(/i)cr(A — /i)d/i = arctansinh7rA , — n < p(X) < . (3.120) 

J-00 2 

From Eas. (|3.11i)l3.120|) one gets 

e= sinp . (3.121) 

The momentum j»t(Ai, A 2 ) varies over the interval [0, 27r), when Ai and A 2 run independently 
over the whole real axis. The spin of this state can be computed by the formula 



S = - {a(X - Ai) + ct(A - X 2 ))dX = 1 . (3.122) 



Let us finally analize the singlet excitation characterized by the configuration 
Denoting by As the only number among the A^i/2 which characterizes the string of length 
2 and by A^ the numbers A J:0 for the strings of length 1, Eas. ((3.91(1 read 



arctan2A 3 - = -^- + — $(A, - As) + — J2 arctan(Aj - X k ) (3.123) 

fc=i 

f-2 

arctanAs = ■= ^ $(X S - \j) (3.124) 



with 



2 
$(A) = arctan2A + arctan-A . (3.125) 
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The y — 2 numbers Qj vary in the segment [— ^ + h,^f — h]] among them there are the 

two holes Q\ and Q 2 ■ Taking the thermodynamic limit one finds the density of A's for 
the singlet 

p(X s ) = p(A) + i(a(A - Ai) + a(X - A 2 ) + c^(A - \ s )) (3.126) 

where p and a were given in Eas. (|3. 1081 IH.llHjl and where 

1 r°° 9 

•w—kL'-^t-i&hsrj ■ m2T > 

In |30| it was demonstrated that the string parameter As is uniquely determined by the A's 
parametrizing the two holes 

Ah) .(h) 

Xs= Al ± 2 ■ (3.128) 

In [23] it was also proved the remarkable fact that the string of length 2 does not contribute 
to the energy and momentum of the excitation, so that the singlet and the triplet have the 
same dispersion relations 

£5(Ai,A 2 ) = e T (Ai,A 2 ) = e(A 1 ) + e(A 2 ) (3.129) 

Ps(Ai,A 2 ) = p T (Ai,A 2 )=p(A 1 )+p(A 2 ) (mod 2tt) . (3.130) 

The spin of this excitation is, of course, zero 

(2cr(A)+w(A))dA = (3.131) 

The only difference between the state whose configuration is given in Eq. (|3.1UU|l and the 
state of Eg. (|3. 99(1 is the spin. 

To summarize, the finite energy excitations of the antifcrromagnetic Heisenberg chain 
are only those belonging to the class Mat an d are described by scattering states of an 
even number of quasiparticles or kinks. The momentum p of these kinks runs over half the 
Brillouin zone —n < p < 0, the dispersion relation is e(p) — 4^ sinp, Ea. ((3.121(1 . and the spin 
of a kink is 1/2. The singlet and the triplet excitations described above are the only states 
composed of two kinks, the spins of the kinks being parallel for the triplet and antiparallel 
for the singlet. For vanishing total momentum all the states belonging to Mat have the 
same energy of the ground state so that they are gapless excitations. Since the eigenstates 
of Hj always contain an even number of kinks, the dispersion relation is determined by a 
set of two-parameters: the momenta of the even number of kinks whose scattering provides 
the excitation. There arc no bound states of kinks. 



3.3.1 Finite size antiferromagnetic Heisenberg chains 

Let us now turn to the computation of the spectrum of finite size quantum antiferromagnetic 
chains by exact diagonalization. We shall see that already for very small chains, the spectrum 
is well described by the Bethe ansatz solution in the thermodynamic limit. Furthermore, 
an intuitive picture of the ground state and of the lowest lying excitations of the strongly 
coupled two-flavor lattice Schwinger model emerges, due to the mapping of the gauge model 
onto the spin chain — see chapter 5. 

The states of an antiferromagnetic chain are classified according to the quantum numbers 
of spin, third spin component, energy and momentum \S, S 3 , E,p >. For a 4 site chain the 
momenta allowed for the states are: 0, ^, 4p mod 2tt. The ground state is 

\g.s. >= 10,0,-37,0= -±=(2| UU> +2| ITIT> -| TTU> -| TUT> -I UTT> -| 4TTI>) 

V 12 

(3.132) 
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This state is P-parity even. In fact, by the definition of P-parity given in En. (|5.29|l . the 
P-parity inverted state ((3.132(1 is obtained by reverting the order of the spins in each vector 



> appearing in 1(3. 132(1. e.g. | UTT>- 



TTU>- 



The A's associated to the ground state (solution of the Bethe ansatz equations 1(3.86(1 ) 
are Ai = —7-7= and A2 = ttt^- There is also an excited singlet 

2yo 2yo 



|o,o,-J,tt >= -^=(| ||TT> -1 ITU> -| TUT> +1 TTU>) 



(3.133) 



It is P-even, so that it is a S p = + excitation, with the same quantum numbers (the 
isospin is replaced by the spin) of the lowest lying singlet excitation of the strongly coupled 
Schwinger model discussed by Coleman [HOI ■ The state 1(3.133(1 also coincides with the excited 
singlet described by the configuration H3.99fl . It has only two complex A's which arrange 



themselves in a string approximately of length 2, Ai = — A2 = i\ 8 l 17 and there are two 



holes with Q\ 



(/,) 



\ and Q 



(h) 



There are also three excited triplets, whose highest weight states are 



1 



ll,l,-J,- > = 


= -A 


|1,1,-2J,7T> = 


1 

V4 


I,W.£> ■ 


1 

7t 



(I ITTT> +*| TITT> -I TUT> -*| TTU>) (3.134) 

(I ITTT> -I UTT> +1 TTIT> -| TTT4>) (3.135) 



(IITTT>-*ITITT>-|TTIT>+«ITTTI>) • (3.136) 



Among these, only the non-degenerate state with the lowest energy has a well defined P- 
parity l|3. 135(1 . It is a S p = 1 _ like the lowest lying triplet of the two- flavor strongly coupled 
Schwinger model. The degenerate states can be always combined to form a P-odd state. 

We thus see that within the states in a given configuration there is always a representative 
state with well defined parity, the others are degenerate and can be used to construct state of 
well defined energy and parity. Moreover the parity of the representative states (with respect 
to the parity of the ground state) is the same of the one of the lowest-lying Schwinger model 
excitations in strong coupling. 



All the triplets in (|3.13(i|l have one real A and two holes; they can be associated with 
the family of triplets (|3. 100(1 . In table ((3.3(1 I summarize the triplet A's and Q^'s. The 



Table 3.3: Triplet internal quantum numbers 



TRIPLET 


A 


Qi'° 


Q^ 


|l,l,-J,f > 


1 
2 


-1 





|1,1,-2J,7T> 





-1 


1 


|l,l,-J,f > 


1 
2 





1 



spectrum exhibits also a quintet, whose highest weight state is 

|2,2,o,o>HTTTT> 



(3.137) 
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E/J 



O- 



{3- 



tt/2 



-e- 



■€* 



■€) 



3tt/2 



-e- 



O QUINTET 
Q TRIPLETS 
□ SINGLETS 



2- 



-*- 



■H 



P 



Figure 3.2: Four sites chain spectrum 



I report in fig. (|3.2|) the spectrum of the 4 sites chain. 

Let us analize the spectrum of the 6 site antiferromagnetic chain. The momenta allowed 
for the states are now 0, f, 3p,7T, 4f, ^ mod 2ir. The ground state is 



\g.s.> = |0,0,- 



■-(5 + Vl3),7r>= 

2 V26 - 6-/L3 



{imut>-itmu> 



+ 



i- yi3 

6 

I UUTT> 
4-a/13 



(I TUTU> -I UTUT> +1 iTUTT> -| TUTU> +1 UTUT> -| ITTITI> 



ITITTI> -I UTTU> +1 iTTUT> -| TTUU> +1 TUUT>) 
(I TTUU> -I TTIIIT> +1 TIUTT> -| UiTTT> +1 IITTTI> -| iTTTU>)} (3.138) 



This state is odd under P-parity. The spectrum of the six sites chain is reported in fig. (j3.3(l . 
There are 9 triplets in the spectrum. In |51| it was already pointed out that the number of 
lowest lying triplets for a finite system with N sites is N(N + 2)/8, so for N = 6 there are 6 
lowest lying triplet states. In order to identify these 6 states among the 9 that are exhibited 
by the spectrum of fig. l|3.3|l . it is necessary to compute their A's and their Q's. In this way 
in fact, I can find out which are the triplets characterized by two holes and thus belonging 
to the triplet of type (|3.10()|l . In table (|3.4|l I report the internal quantum numbers of the 
lowest lying triplets. The Q( h ''s vary in the segment [— |, |]. The highest weight state of 
the triplet of zero momentum and energy —(J/2) (5 + v5) reads 



|0,0,-^(5 + \/5),0> 



-A- 



-3 + V5 



^45 - 15\/5 2 

|TTUTT>+|TUTTT>) 

(I mm> +i mm> +i mm> +i tttiti> 

(i - V5)(l ITTITT> +1 TUm> +1 TITTIT>)} 



(1 44TTTT> +1 mm> +1 TTTTI4> +1 TTTiit> 



tnnt> 



nntt>) 

(3.139) 



One can get the triplet of energy —(J/2) (5 — \/5) from (|3.139|l by changing y/Z — > —y/Z, As 
can be explicitly checked from l|3.139l) . the two non-degenerate triplets of zero momentum 
are then P-parity even, namely they have opposite parity with respect to that of the ground 
state, as it happens for the lowest lying triplet excitations of the two-flavor Schwinger model. 
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For what concerns the degenerate triplets of momenta 7r/3 and 57r/3 (or 2tt/3 and 4tt/3) 
they do not have definite P-parity, but it is always possible to take a linear combination of 
them with parity opposite to the ground state. 



Table 3.4: Triplet internal quantum numbers 



TRIPLET Ai A 2 Q^ Q { 2 h) 




M 1 5+\/5 7 n - J5-2V5 / 5-2^5 1 1 
1^-' 2 «,w^ y 20 V20 2 2 




11 1 5-\/5 7 a x /5+2V5 J5+2VE 3 3 
I- 1 ' ' 2 J, u ^ V20 V20 2 2 


Mi 5rir x n/S+v/tt Vl-v^r 3 1 
I 1 ' 1 ' 2 > 3 ^ 8 8 2 2 




|i i 7+vT7 r 2ir x -2^-^-2+2^17 -2-\/3+V-2+2\/l7 3 1 
1 ' ' 4 J ' 3 " 2+2VT7 2+2 v / 17 2 2 




M i 7+VT7 7 4ir x 2V3-V-2+2VT7 2V3+V-2+2VT7 1 3 
1 ' ' 4 J > 3 " 2+2VT7 2+2 v / 17 2 2 


|1 1 5 7 5tt ^ v^-^ v^+V^ 1 3 
|l, 1, 2 J i 3 •-* 8 8 2 2 



The remaining three triplets in fig. (|3.3l) have no real A's and are characterized by a string 
of length 2 and four holes for Q = — I, — |, |, |, i.e. do not belong to the type (|3.1UU|) . More 
precisely, two triplets have a string approximately of length 2, due to the finite size of the 
system, while the triplet of momentum 7r has a string exactly of length 2. In table (|3.5|) I 
summarize the quantum numbers of these triplets. 



Table 3.5: Four holes triplet internal quantum numbers 



TRIPLET 



Ai 



1.1, 



7— y/l7 7 2i. 
4 J > 3 ^ 



2V3-iV '2+2x717 2y / 3+i- s /2+2xAT7 



-2+2V17 



-2+2V17 



11,1,— J,7T > 



11,1,-2=^2.7, £ > 



2y / 3+t\/2+2v / T7 2\/3-t\/2+2\/T7 



2-2V17 



2-2V17 



In fig. (|3.3l) it is shown that the spectrum exhibits five singlet states. The lowest lying 
state at momentum -k is the ground state. Then there are three excited singlets characterized 
by the configuration with two holes 1)3. 99|) . i.e. they have one real A and a string of length 
almost 2. In table l|3. fij) I summarize their quantum numbers. Among these singlets, those 
which are not degenerate, have P-parity equal to that of the ground state (odd) as it happens 
in the two-flavor Schwinger model. The non- singlet in fact reads 



|0,0,-3J,0> = 



{I mm> +1 utut> +1 initt> +i mm> +1 umt> +i itnn> 



- i unn> -i mm> -i nmi> -i mut> -i mm> -i ninr>} .(3.140) 
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The degenerate singlets are again not eigenstates of the P-parity, but it is always possible 
to take a linear combination of them with the a P-parity that coincides with that of the 
representative state <|3.14U[1 of the configuration. 



Table 3.6: Singlet internal quantum numbers 



SINGLET 


A 


As 


QP 


Q 2 ' l) 


|0,0,-3J,0> 








-1 


l 


|0,0,-2J,f > 


14 


-2+3^ 
V3(4+v^) 





l 


|0,0,-2J,f > 


14 


2-3\/2 

V3(4+v^) 


-l 






The remaining singlet 1 0, 0, — 5 ^ 13 J,tt > it is not of the type 
acterized by a string approximately of length 3 with A^i 

•^3,1 = - 



5+2y / 13 
12 ' 



It is char- 
A2 1 = and 



5+2VT3 



12 



Even in finite systems very small like the 4 and 6 sites chains, the "string hypothesis" 
is a very good approximation and it allows us to classify and distinguish among states with 
the same spin. 



combination of all the 



states with S spins up and ^ spins down. These states 



The ground state of the antiferromagnetic Heisenberg chain with N sites is a linear 

' N ' 

N 
2 

group themselves into sets with the same coefficient in the linear combination according to 
the fact that the ground state is translationally invariant (with momentum (it) for — even 
(odd)), it is an eigenstate of P-parity and it is invariant under the exchange of up with down 
spins. The states belonging to the same set have the same number of domain walls, which 
ranges from N, for the two Neel states, to 2 for the states with ^ adjacent spins up and ^ 
adjacent spins down. 



The ground state of the 8 sites chain is 



hs. >= 



where 



(1^8 > +a|4 1} > +/?l4 2) > +l\4 1} > +<5|V>i 2) > +e|^i 3) > +CIV>2 >) 



(3.141) 



l^8>=innnn>+iiTinriT> 




(3.142) 


\4 1} >= 1 TUUTU> +1 UUUTT> +translated 


states 


(3.143) 


IV4 2) >= 1 TTmm> +1 UUTUT> +translated 


states 


(3.144) 


1^4 >= TT44TT44> ^translated states 




(3.145) 


IV>f >= 1 TTTUm> +1 1UTTUT> translated 


states 


(3.146) 


\4 3) >= 1 TTIIITU> +1 UTTTUT> +translated 


states 


(3.147) 


\ip2 >= TTTTUU> +translated states . 




(3.148) 
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E/Jf 



A 



tt/3 



2tt/3 



7T 



4tt/3 



5tt/3 



2^ 

7±r 



-1/2 



(5- yggVg 



-(> 



-(7- V17)/2' 



e 



-(5- \/5)/2 



^5 



■G 



O 



O 



p 



O 



-3/2 



-5/2 



■e 



-(Vl7 + 7)/4 



■^> 



■e 



& 



o 



■e 



€> 



e 



-(VS + 5)/2 



€> 



-(Vl3 + 5)/2 
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Figure 3.3: Six sites chain spectrum 



By direct diagonalization one gets 



a - 


= -0.412773 


P -- 


= 0.344301 


1 = 


= 0.226109 


5 = 


= -0.087227 


e = 


= 0.136945 


c = 


= 0.018754 


N -- 


= 2 + 16a 2 + 



8/3 2 + 4 7 2 + 16(5 2 + 16e 2 + 8C 2 = 6.30356 



(3.149) 
(3.150) 
(3.151) 
(3.152) 
(3.153) 
(3.154) 
(3.155) 



The energy of the ground state is 



E g , s , = -5.65109J 



(3.156) 



;in2 



Eq. (|3.156|l differs only by 1.8% from the thermodynamic limit expression E g . s . = 
—5.54518. Moreover also the correlation function of distance 2 Ea. 1)3. 158(1 computed for the 
8 sites chain is G(2) = 0.1957AT, value which is 7% higher than the exact answer Ea. 1)3. 159(1 . 

In the analysis of finite size systems I were able to find the coefficient (3 of the first set 
of states containing N — 2 domain walls in the ground state. These states are obtained 
interchanging two adjacent spins in the Neel states. The (3 is for a generic chain of TV-sites 







N + 2E g 
N 



l-21n2 



(3.157) 
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3.3.2 Spin-spin correlators 

The explicit computation of spin-spin correlations is far from being trivial since the correlator 
G(r) =< g.s.\So ■ S r \g.s. > is not known for arbitrary lattice separations r. For r — 2 it was 
computed by M. Takahashi 05| m his perturbative analysis of the half filled Hubbard model 
in one dimension. For r > 2 no exact numerical values of G(r) are known. In ^ 5*2") were given 
two representations of G(r), while in |53l 154) the exact asymptotic (r — ► oo) expression of 
G(r) was derived. 

In order to explicitly compute the second order energies Eq. (|5.62|) and Eq. (|5.63[l one has 
to evaluate the correlation function 



1 N 
G{ ^ =wE< 9-s.\S x ■ S x+2 \g.s > (3.158) 



x=l 



which has been exactly computed in 42 j and is given by 

G(2) = ~(l-161n2 + 9C(3))= 0.1820 . (3.159) 

In the following I shall show how the knowledge of this correlator allows one to compute 
explicitly the first three "emptiness formation probabilities" , used in Ref. [S2j in the study 
of the Heisenberg chain correlators, G(r). The isotropy of the Heisenberg model implies 
that 

N N 

J2<9-s.\S x -S x+2 \g.s. >=3j2<9-s.\Sl-S 3 x+2 \g.s. > . (3.160) 

X— 1 X — 1 

Let us introduce the probability P3 for finding three adjacent spins in a given position in the 
Heisenberg antiferromagnetic vacuum. Taking advantage of the isotropy of the Heisenberg 
model ground state and of its translational invariance, it is easy to see that the correlator 
(I3.160[) can be written in terms of the Pa's as 

N 1 

£ < g.a.\Sl ■ S 3 x+2 \g.s >=N-2( P 3 (TTT) + P^UT) - ^(TTI) - P^UTT) ) • (3.161) 

x = l 

The factor 2 appears in (|3.161|l due again to the isotropy of the Heisenberg model: the 
probability of a configuration and of the configuration rotated by 7r around the chain axis, 
are the same. 

In |S2] the so called "emptiness formation probability" P{x) was introduced. 

X 

P(x)=<g.s.\l[P j \g.s.> , (3.162) 



where 



Pj = -(crf + l) (3-163) 



and <jj is the Pauli matrix. P(x) determines the probability of finding x adjacent spins up 
in the antiferromagnetic vacuum. One gets 

P'(TTT) = P(3) (3.164) 

PflTT) = P(1)-2P(2) + P(3) (3.165) 

P(ITT) = P(TTI) = P(2) - P(3) (3.166) 



so that Eq. (|3.158|l reads 



G{2) = 2P(3) - 2P(2) + \p{1) ■ (3.167) 
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Using the exact value the correlator G(2) computed in 02] from ((3.167(1 and from the known 
values of P(2) and P(l) given in [S] 

P(l) = \ (3.168) 

P(2) - i(l~ln2) (3.169) 

(3.170) 
one gets 

P(3) = i(l_71n2) + |c(3) ■ (3.171) 

For the general emptiness formation probability P(x) of the antiferromagnetic Heisenberg 
chain, an integral representation was given in |52] . but, to my knowledge, the exact value of 
P(3) l|3. 171(1 was not known. 

I illustrate now the computation of a spin-spin correlator which appears in the mass 
spectrum of the two-flavor lattice Schwinger model 

N 



< g-s.\V ■ V\g.s. >= ^2 (< 9-S-\(S x ■ S y )(S x+1 ■ S y+1 )\g.s. > 

x,y=l 

N 

- < g.s.\(S x ■ S y+1 )(S x+1 ■ S y )\g.s. >) - ^ < 9- s -\S x ■ S x+1 \g.s. > .(3.172) 



x=l 



It is possible to extract a numerical value from Ea. l(3.172|l only within the random phase 
approximation [421 155| . For this purpose it is first convenient to rewrite the unconstrained 
sum over the sites x and y as a sum where all the four spins involved in the VEV's lie on 
different sites, 

< g.s.\V ■ V\g.s. > = ^2 (< g.s.\(S x ■ S y )(S x+1 ■ S y+1 )\g.s. > 

y^x±l 

- 3 

- < g.s.\(S x ■ S y+1 )(S x+1 ■ S y )\g.s. >) + -N 

N N 

- ^ 5Z < 9- s -\S* ' S x+ i\g.s. > -^2 < g-s.\S x ■ S x+2 \g.s. >(3.173) 

X— 1 X— 1 

and then factorize the four spin operators in Eq. (|3.173|1 as 
< g.s.\V ■ V\g.s. > = N^(< g.s.\S Q ■ S r \g.s. > 2 - < g.s.\S Q ■ S r+1 \g.s. >< g.s.\Si ■ S r \g.s. >) 

r=2 

3 1 N N 

+ ^N--J2<9-s-\S x -S x+1 \g.s.>-J2<9-s-\S x -S x+2 \9.s.> (3.174) 

x— 1 a; — 1 

Of course, Eg. (|3. 174(1 provides an answer larger than the exact result; terms such as < 
(. ..)(...) > yield negative contributions which are eliminated once one factorizes them in 
the form < (...) >< (...) >. This is easily checked also by direct computation on finite size 
systems. 

The spin-spin correlation functions G(r) are exactly known for r = 1,2. For the spin-spin 
correlation functions G(r) up to a distance of r = 30 the results are reported in table 1(3.7(1 

E31ESI 

For r > 30, one may write |53( 
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Table 3.7: Spin-spin correlation functions 



r 


G(r) 


r 


G(r) 


1 


-0.4431 


16 


0.0305 


2 


0.1821 


17 


-0.0296 


3 


-0.1510 


18 


0.0274 


4 


0.1038 


19 


-0.0267 


5 


-0.0925 


20 


0.0249 


6 


0.0731 


21 


-0.0242 


7 


-0.0671 


22 


0.0228 


8 


0.0567 


23 


-0.0223 


9 


-0.0532 


24 


0.0211 


10 


0.0465 


25 


-0.0206 


11 


-0.0442 


26 


0.0196 


12 


0.0395 


27 


-0.0193 


13 


-0.0379 


28 


0.0183 


14 


0.0344 


29 


-0.0181 


15 


-0.0332 


30 


0.0172 



™ <> 3 /~2~ 1 m 3 f Mi 156C(3)-73 ,, 3 «, , , 4 , 

GW = 4 V ^^JW) [ " ~* 9{) + 384 *(r? + 0(g(rf) 



^((-ir + l + 0(.g(r)) + 0(l)] 
2r r* 



with g{r) satisfying 



and 



9(r)= CW (1+ 2 3(r)ln(9(r))) 



C(r) = ln(2V2^e 7+1 r) . 
Ea. (|3.176|l may be solved by iteration. To the lowest order in ij one finds 

1 1 



9 ^'~ C(r) C(r) 2 
Inserting l|3.178|l in Eq. (|3.175|) leads to 



lnC(r) 



G(r) « ^tt^v^WII + ^y InC(r)] + 0(-^) 



Inserting Ea. l|3.179|) in l|3.174|) . one finally gets 

<g.s.|y-f|t/.s. >=0.3816iV 



(3.175) 

(3.176) 

(3.177) 

(3.178) 

(3.179) 
(3.180) 



Last, I report the following exact three spin correlators that have been used in the 
determination of the mass spectrum of the two-flavor lattice Schwinger model 



N 

<g.s.\J2 S * S *+A2\9-s-> = 

x=l 

N 

<g.s.\J2 S xS^ +1 Sl +2 \g.s.> = 

x=l 



(3.181) 
(3.182) 
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A' 



<g.s.\J2S-S+ +1 S 3 x+2 \g. S .> = (3.183) 

x=l 

N 

<g.s.\J2 S x S l + i S x+2\g-s-> = (3.184) 

x=l 
N 

<g.s.\Y,SZSl +1 St +2 \g.s.> = (3.185) 



So, on a spin singlet, not only the VEV of J2 x =i ^x * s zer0 ; but also every VEV with an 
odd number of S 3 . 



3.4 SU(J\f) quantum antiferromagnetic chains 



It is my purpose to introduce spin-1/2 antiferromagnetic Heisenberg chains where "spins" 
are generators of the SU(Af) group. In the limit J\f = 2 one has the usual antiferromagnetic 
Heisenberg chain discussed in the previous section. An U(Af) spin-1/2 quantum antiferro- 
magnetic chain is described by the Hamiltonian 

JV 

H U j {N) = jY / S ab (x)S ba (x + 1) (3.186) 

x=l 

where S ab (x) a, b = 1, . . . , J\f are the generators of U{M) satisfying the Lie algebra 

[S ab (x), S cd (y)} = (S a d(x)6bc - S cb {x)6 ad )S xy (3.187) 

and they can be conventionally represented by fermion bilinear operators 

S ab (x)=i;i x ^ bx - 5 -f . (3.188) 

The representation of the algebra on each site is fixed by specifying the fermion number 
occupation 

P {x) = Y, s ^ x ) ■ ( 3 - 189 ) 

a=l 

By fulfilling the global neutrality condition X^=i P( x ) ~ 0> one ma y choose 

&!.*-{/„ I'™ (3.190, 

or viceversa the opposite choice for x even-odd. Eq. I|3. 190(1 restricts on each site the Hilbert 
space to a representation with Young tableau of m rows for x even and J\f — m rows for x 
odd. For each site x p(x) is the generator of the U(l) subgroup of U(Af). 

Let us use the basis T a = (T a )* , a = 1, . . . , M 2 - 1, of the Lie algebra of SU{N) in the 
fundamental representation such that tr(T a T f3 ) = S afj /2 and [T a ,T^] = if a ^T^, where 
fa/3~f are ^g s t ru cture constants. By means of 

T: b T? d = \5 ad 5 bc - ^SabScd (3.191) 

and redefining the group generators 

S: = VLWto (3.192) 

GO 



one can rewrite the Hamiltonian 1(3.186(1 as 

N 

hT ] =JY.P^ + ^+ hS J U(N) (3-193) 



x=l 

where 

N 

Hj UW = JE S "? + i (3.194) 



is the Hamiltonian of an SU(Af) quantum antiferromagnet. From Ea. 1(3. 193(1 it is clear that 
by fixing p(x) Hj is reduced to Hj 

The representations of the SU(AT) spins relevant for the relationship of the spin Hamil- 
tonians ((3.194(1 with the Schwinger models are different when M is even or odd. When TV" 
is even I shall consider the representation on each site with Young tableau of one column 
and A/"/2 rows so that p{x) — 0. When M is odd I shall consider on one sublattice the 
representation with Young tableau of one column and (A/" + l)/2 rows and of one column 
and (N — l)/2 rows on the other sublattice. 

As an explicit example let us consider the case M — 3. The generators T a of the SU(3) 
group in Eg. 1(3. 192(1 are the Gell Mann matrices. By denoting u,d and s the three flavors, 
the eight spin operators 5" read 

S x = ^Lfe+V-LVw (3.195) 

Si = -ii>i x ipdx + iw dx i>ux (3.196) 

Si = ^LV-W-^Lfe (3.197) 

St = <ALV^+CVw (3.198) 

S x = -#L^*+#LV^ (3-199) 

S 6 X = ^LV^+Cfe (3.200) 

Si = -#L^, + #Lfe (3-201) 

S 8 X = i(C^+^Lfe-2VLV' s:c ) (3.202) 

Let us find the ground state of the Hamiltonian ((3.194(1 for a chain of two sites. Taking the 
representations with one particle on site 1 and two particles on site 2, the ground state with 
energy E g , s , = — 16J/3 reads 

\g.s.>=-L(\u d s >-\d U s >+\s U d >) (3.203) 

The state (|3.203(l is a singlet of SU(3), i.e. it is annihilated by the Casimir S 2 = 0i + S 2 ) 2 = 
12 a =i [(Si) 2 + (Si) 2 + 2S1S2] ■ If one chooses the representation with two particles on site 
1 and one particle on site 2, the ground state degenerate with ((3.203(1 is 

\g.s. >'= -J=(| * U>-\ U s d>+\ U d s>) (3.204) 

Diagonalizing the translation operator T with the Hamiltonian 13.194(1 , one gets two degen- 
erate ground states 

\ G .S. >± = ' g - 3 - > ^ S - >' (3.205) 

In Eq. l(3.205|) the linear combination with the + has momentum zero and the one with 
the — has momentum 7r. By studying this very simple example, one can infer that the 
thermodynamic limit N — > 00 analysis for a generic SU(Af) chain is very involved. 



61 



Unfortunately, no analysis with a level of completeness such as the one given in [23 f° r 
the SU(2) case does exist for a generic symmetry group SU(Af) with spins in the represen- 
tation in which I am interested. In [S7| SU(Af) antiferromagnetic models were solved for a 
particular spin representation such that the Hamiltonian (|3.194|l becomes 

N 

h suW)^jJ2p xx+1 (3.206) 

x=l 

with P X x+i the operator that permutes whatever objects occupy sites x and x + 1. The 
spectrum of Eq. (|3.206|l was shown to exhibit massless excitations. 

A large Af expansion approach has been performed in |58j for an SU(Af) antiferromag- 
netic chain characterized by spins living in a representation with Young tableaux of one row 
on one sublattice and N — 1 rows on the other sublattice. In the case of Young tableaux 
of one column it was found that the ground state is twofold degenerate and breaks transla- 
tional and parity symmetry. Moreover the elementary excitations are massive non relativistic 
solitons in the large Af limit with a mass of 0(Af). 

In (SI] the Lieb-Shultz-Mattis theorem jHO] was generalized to SU(Af) spin chains. The 
theorem proves that a half-integer-S spin chain with essentially any reasonably local Hamil- 
tonian respecting translational and rotational symmetry either has zero gap or else has de- 
generate ground states spontaneously breaking translational and parity invariance. In [SU| 
it was proved the existence of a unique ground state of the SU(Af) chains where Af is even 
for spins in the antisymmetric Af tensor representation. Under the assumption of a unique 
ground state which must be a SU(Af) singlet, an infinitesimal energy gap was found for all 
the representations of SU(Af) whose Young tableaux contain a number of boxes not divisible 
by Af. Of course the degeneracy of the ground state trivially implyes a zero gap. 

Using the Lieb-Shultz-Mattis I can state that when Af is even the ground state of the 
Hamiltonian (|3.194|l is unique and there are gapless excitations for the spin representation 
with a Young tableau of one column and Af/2 rows. When Af is odd the ground state is 
twofold degenerate, as I illustrated with the simple example of the two site SU(3) chain. It is 
doubtful if there exist gapless excitations in this case and at present time I am investigating 
this problem. 
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Chapter 4 

The one-flavor lattice Schwinger 
model 



In this chapter I study the one-flavor lattice Schwinger model in the strong coupling limit 
using staggered fermions and the hamiltonian approach to lattice gauge theories [H] ■ Even 
though the continuous axial symmetry is broken explicitly by the staggered fermions, the 
discrete axial symmetry remains and appears in the lattice theory as a translation by one 
site. In the continuum gauge model the chiral symmetry is broken by the anomaly. An 
interesting issue which arises in the lattice regularization of gauge theories is how the lattice 
theory produces the effects of the axial anomaly. Lattice field theory are manifestly gauge 
invariant by construction and therefore, to produce continuum theories, they must find a 
way to violate axial current conservation. Normally, in lattice gauge theories axial anomalies 
are either cancelled by fermion doubling or else the lattice regularization breaks the axial 
symmetry explicitly. The lattice Schwinger model in the hamiltonian formalism with one 
species of staggered fermions represents the unique example where neither of these occur. 
The effects of the anomaly are not cancelled by doubling since the continuum limit of (1 + 1)- 
dimcnsional staggered fermions produce exactly a Dirac fermion which is the matter content 
of the model. Furthermore, even if the continuum axial symmetry is explicitly broken, the 
remaining discrete axial symmetry must be spontaneously broken. The mass operator ipijj 
is odd under discrete axial transformations, if this symmetry is unbroken then < 'ijjift >= 0. 
However, in the continuum theory it is known that 62 , 

<V^(x))=-£-^= . (4.1) 

where 7 = 0.577... is the Euler constant. The expectation value in Eq. (|4.1|l is the chiral 
condensate. LJ I shall show that the chiral condensate i|4.1|) is non-zero also on the lattice 
and the discrete axial symmetry is indeed spontaneously broken. In addition to the analysis 
of the chiral symmetry breaking pattern on the lattice, by performing a strong coupling 
expansion around a gauge invariant ground state I shall derive results converging more 
rapidly to the continuum than the results provided in |63j . The one- flavor Schwinger model 
has also been used as an example of the fact that quantum link models may reproduce the 
physiscs of conventional hamiltonian lattice gauge theories 65 . 

J It is worth to stress that the sign of the chiral condensate in Ea. 14.11 is just a convention. The chiral 
condensate has always the opposite sign of the fermion mass m that appears in the Lagrangian of the theory: 
it is the order parameter for a spontaneous symmetry breaking of the chiral symmetry. The chiral condensate 
is the analogous of the magnetization for a spin model and the mass m plays the role of the magnetic field. 
So even if one is studying a massless model, in order to decide the sign of the chiral condensate, one has to 
consider a small mass term m that then one sends to zero. The minus sign which appears in Eg. 14. II is due 
to the fact that one considers the m > 0+ limit in the Schwinger model Lagrangian. 
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4.1 Definition of the model 



The continuum Schwinger model |28l IM] is the prototypical example of a solvable model 
where the anomaly occurs. It is 1+1-dimensional electrodynamics with a single charged 
massless Dirac spinor field. The action is 

S = J d 2 a#( V + 7 ^)</> - 4^2 V") (4-2) 

It is invariant under gauge transformations, 

iP(x) — ► e ie{x) iP{x), ^(x) — > ip^{x)e' ie{ - x) (4.3) 

A^ — » A„ + c^(x) (4.4) 

and, formally, under the axial phase rotation, 

ip(x) — ► e i " lsa tp(x),^(x) — ► ^(x)e-' nsa (4.5) 

(where 7 5 = ij ^ 1 ). 

At the classical level the above symmetries lead to conservation laws for the vector and 
axial currents, 

f{x)=^{x) 1 ^{x), j|(x)^(x) 7 V^(x), (4.6) 

respectively. At the quantum level both currents cannot be simultaneously conserved. If 
the regularization is gauge invariant, so that the vector current is conserved, then the axial 
current obeys the anomaly equation 66 

dMx) = $^F^x). (4.7) 

Z7T 

Moreover, the correlation functions of the model do not exhibit axial symmetry. 

There are only neutral particles in the spectrum of the Schwinger model; notably, the 
bound state boson created by the axial current operator, 

3^(x) = V^J^x) (4.8) 

appears in the spectrum as a free pseudoscalar with mass -y=. 

In terms of the boson field, the action is 

S= />z(^<OT + |U 2 ) (4.9) 

and the correlation functions for the vector and axial currents are obtained from the corre- 
lation functions of derivatives of the bose fields. The coupling constant e c has the dimension 
of a mass and the model is super-renormalizable. 

Since the Schwinger model is solvable it has often been used as a field theory where 
methods of lattice gauge theory, particularly the strong coupling expansion can be tested 
and compared with exact results of the continuum model |63| . 

In this chapter, I shall revisit the lattice quantization of the Schwinger model. My main 
purpose is to demonstrate two results. One is to show how the effects of the anomaly appear 
through spontaneous symmetry breaking. The other is to demonstrate that a strong coupling 
expansion around a gauge, parity and charge conjugation invariant lattice ground state can 
produce results in good agreement with the continuum. My results turn out to be quite 
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different from those of previous strong coupling computations 1631 . The difference arises 
from my careful treatment of the normal ordering of the charge density operator, which 
must be defined so as to be compatible with all of the discrete symmetries of the theory. 
The necessity of careful normal ordering of the charge operator was pointed out in |24| . My 
vacuum state, vacuum energy and energies of elementary excitations are different from those 
in [631 167) and, as I shall show, they give a more accurate extrapolation to the continuum 
limit. 

The Hamiltonian and gauge constraint of the continuum Schwinger model are 

H = J ' dx(^-E 2 (x)+^(x)a(id x + eA{x))iP(xy) (4.10) 

d x E(x) + ^ {x)i/)(x) -0 (4.11) 

A lattice Hamiltonian and constraint which reduce to these in the continuum limit are: 

Hs = ^ E E l ~ £ E(V4 + ie^, - ^e-^^+i) (4.12) 

X X 

E x -E x ^ + i>ti> x --~Q, (4.13) 

where the fermion fields are defined on the sites, x = — y, — | + 1,...,|, gauge and the 
electric fields, A x and E x , on the links [x; x + 1], N is an even integer and, when N is finite, 
I use periodic boundary conditions. When N is finite, the continuum limit is the Schwinger 
model on a circle (HSJ E2 • The coefficient t of the hopping term in H4.12[) plays the role of 
the lattice light speed. In the naive continuum limit, ex = e c and t = 1. However, I shall 
keep ex and t as parameters which can be adjusted to fit the lattice values of quantities such 
as the mass gap and the chiral condensate to those which are known in the continuum. 

The non vanishing (anti-)commutators for the lattice variables are 

[A X ,E V ] = iS xy , {ip x , iph = S xy (4.14) 



The Hamiltonian and gauge constraint exhibit the discrete symmetries 

• Parity P: 

A x —+ -A_ x _ u E x — - -£U_i, 4>x — » (-W-*. 4 — ► {~l) x rl>U (4-15) 

• Discrete axial symmetry T: 

A x — ► A x+1 , E x — ► E x+1 , t/j x — > i> x +i, i>t — ► i> x+ i (4-16) 

• Charge conjugation C: 

A x — > -A x+ i, E x — ► -E x+ i, ip x — > ipt +1 , ipt — * ipx+i (4-17) 

When A x — 0, the spectrum of the hopping Hamiltonian is e(pa) = tsmpa(p £ [0, — ]). 
In the low energy limit it resembles a massless relativistic spectrum for excitations near 
the Fermi level. The two intersections of the energy band with the Fermi level provide the 
continuum right and left moving massless fermions. The electron ground state is invariant 
under charge conjugation only when the Fermi level is at e(pp) = (i.e. when exactly half 
of the fermion states are filled and ^2 X < p(x) >— 0). In the remainder of this paper I 
consider only this case of half-filling. 
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One has two possile ways to put the spinors ip = 



1p2 



on the lattice in the staggered 



fermion formalism. One can put upper components on even sites and lower components on 
odd sites or viceversa. The mass operator, which reduces to ipi]) in the continuum limit is 
the staggered charge density 



M(x) = t-±)l (^ _ 1/2 ) 



(4.18) 



for the first choice of staggered fermions or 



M (x) = ( ^- (^ - 1/2) 



(4.19) 



in the other case. 



These operators l|4. 1811 1|4. 1911 are scalars under parity and are even under charge conju- 
gation but change sign under a discrete axial transformation. 

The lattice Schwinger model is equivalent to a one dimensional quantum Coulomb gas 
on the lattice. To see this one can fix the gauge, A x = A (Coulomb gauge). Eliminating 
the non-constant electric field and using the gauge constraint, one obtains the effective 
Hamiltonian 



Hs — H u + H p = 



lh_ E *+ e -^Yp{ X )V{x-y)p{y) 



2N 



2 ^ 



x.y 



it 
2a 



E^i+i 



where the charge density is 



e lA ^ x 



V4e-^Vx+i) 



p{x) = ipiip x - 



i 



and the potential 



1 



JV-l 

V y ' N ^ 4sin 2 ^ 

n— 1 IS 



(4.20) 



(4.21) 



(4.22) 



is the Fourier transform of the inverse laplacian on the lattice for non zero momentum. The 
constant electric field is normalized so that [A, E] — i . The constant modes of the gauge 
field decouple in the thermodynamic limit N — ► oo. The gauge fixed Hamiltonian l|4.20|l 
can also be written as a quantum spin model. Consider the Jordan- Wigner transformation 

i/> x = l[[2iS 3 (y)]S-(x) , V4 = H[-2iS s (y)]S+(x) , i>U> x = S 3 (x) + 1/2 (4.23) 



y<x 



y<x 



Then, the Hamiltonian is 



H = 



t 



^r 



2a z ^ 



(e iA S + (x + l)S-(x) + e- lA S-(x + l)S+{x))+- 



\S 3 (x)V(x-y)S 3 (y)+e 2 L E 2 /2N 



x x.y 

(4.24) 
For the moment, let us ignore the constant modes E and A. The first term in the Hamiltonian 
is the quantum x-y model which has a disordered ground state (of course, exact solution of 
this model is obtained by the inverse Jordan- Wigner transformation |6()l 169] ) . The second 
is a long-ranged Ising interaction. I will prove in the next section that, in spite of the low 
dimensionality of this system, the latter term has a Neel ordered ground state. This is a 
result of the infinite range of the Coulomb interaction. When both terms are present in the 
Hamiltonian, they compete, the first one favoring disorder and the second one favoring order. 
Later, when I extrapolate my strong coupling results to the continuum limit, I shall assume 
that the ordered ground state persists for all positive values of the dimensionless constant 
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e\a 2 /t. The strong coupling limit is where this constant is large and the ground state has 
Neel order. The continuum limit of the Schwingcr model is found where this constant is 
small. I shall assume that the Neel order persists in that limit. |3 



4.2 The low-lying excitation spectrum 



I shall examine the Schwinger model in the strong coupling limit. To solve the lattice 
Schwinger model in the strong coupling expansion, it is necessary to find states which are 
annihilated by the generator of gauge transformations 

g x = E x - E x _ x + ^J> x - - (4.25) 

and, at the same time, are eigenstates of the unperturbed Hamiltonian 

X 

I shall quantize the gauge fields in the functional Schrodinger picture where the wave- 
functions are functions of the variables A x and the electric field operators are defined as 

i dA x 

In [21] it was proved that there were two gauge invariant ground states of the system. In 
fact, when the number, N, of lattice sites is even the two states 

\ip> = ]J -0i|O > g- 1 ^- A - 1 ^ n - (4.28) 

x=even 

\X> = LI ^ x \0>e^* A *^ 1 (4.29) 



odd 



are degenerate gauge invariant ground states of H u . They are gauge invariant since are 
annihilated by the operator Q x . To see that they are indeed ground states, note that when 
H u acts on gauge invariant states, it is effectively the sum of two operators, the zero mode 
energy e 2 L E 2 /2N where E = J2 X E x /y/N and the Coulomb energy 



Ha-° LU ^ 



- ^ p{x)V{x - y)p{y) = ^§J2 P(- k )T^TTP( k ) ■ ( 4 - 30 ) 

2 L — ' 2iV f— -( 4 sin % 



where I have used the Fourier transform 

^x-^ = ^E e ^W^=^»"e{-^,...,f-l} (4.31) 

p 

Furthermore, gauge invariant states have vanishing total charge, 

X>£^x - \) = (4.32) 

X 

The Coulomb energy satisfies the bound 

= £ l^^ ( ^^-^ )2 = eiaiv/32 (4 - 33) 



2 If it did not persist, but there was a second order phase transition at some critical value of the constant 
to a disordered state, the continuum limit would be found at that phase transition. 

G7 



(since (ip^ip x — \) 2 = 4)- The states \if) > and \x > are eigenstates of the charge density 
with eigenvalues p(k) = ±^8k,-K. These charge densities saturate the lower bound for the 
Coulomb energy. Furthermore, E\tp >= = E\\ > so \ip > and \\ > minimize the energy 
of the zero mode. Therefore, they are degenerate ground states of H u . 

The ground states are invariant under parity and charge conjugation but they both 
break chiral symmetry, i.e. the symmetry under translations by one site. The ivariance 
under parity is easy to see 

P\ip>= JJ (-l) :c '0i ;c e- i2: -~ j4 — - 1 ^^|0 >= |V > (4.34) 

x=even 

and analogously 

P\x>=~\x> (4.35) 

For what concerns charge conjugation one has 

x=even 

(4.36) 



a;— even 



where the last equality is true since, being 



one has 



]T<0|p(:r)|0>=-^ (4.37) 

X 

£ < 0\C- l p(x)C\0 >= ^ (4.38) 



n ^i°>= n ^+icio> (4.39) 



x— even x—even 



Analogously one finds 

C\ X >= |x > (4.40) 

Under discrete chiral symmetry one gets 

r|^ >= Ix > ; r| x >=|V> (4.41) 

The order parameter of the chiral symmetry, the mass operator M = -^- ~^2 x Al(x), is 
diagonal in these states and has eigenvalue \. |f| In order to have a chiral condensate of 
negative sign as in the continuum, one has to consider the mass operator M'(x) — —M(x). 
It is interesting that the ground states are degenerate and are not invariant under chiral 
symmetry, even in finite volume and in one spatial dimension; it is possible due to the 
long-ranged nature of the Coulomb interaction 03] . The chiral symmetry is spontaneously 
broken in the infinite volume limit. It is essential that the number of sites of the lattice, N, 
is even. 

Of course once perturbations from the hopping term in the fermion Hamiltonian are 
taken into account, when the volume is finite, there are high order terms (of order N) in 
the strong coupling perturbation theory which mix the ground states that are constructed 
by perturbing |^> > and \x >■ The mixing amplitude is 

1 N 2 \ f \N aN-1 N \ 

Since the amplitude 14.42JI is negative and real, whatever the diagonal elements are, the 
state with lower energy is the state 

^(I^>+|X>) (4-43) 

3 One has to consider M(x) given by Eg. 14. 181 when applying it to \ip > and M(x) given by Eg. 14. 191 
when applying it to \\ >. 
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This is an eigenstate of the discrete axial transformation, with eigenvalue 1. Thus the true 
ground state is chirally symmetric. 

In order to avoid this restoration of the symmetry, it is necessary to take the infinite N 
limit so that the two ground states are not mixed to any finite order of perturbation theory. 
Thus, in the perturbative expansion one only has to consider diagonal matrix elements 
and, consequently, perturbation theory for non-degenerate-states. I shall then choose as 
unperturbed ground state \ip >. The results concerning the energies would have been the 
same had I chosen \\ >. 

Excitations are created by the operators 

2 

and have energy -k-\i — j\ greater than the ground state energy. A single electron, in order 
to have a gauge invariant state, must be attached to a line of electric flux which goes from 
the electron to infinity. This is not allowed at all in a finite volume with periodic boundary 
conditions and on the infinite lattice, the string of electric flux would have infinite energy. 

2 

Thus electric charges are confined and the string tension is given by the electric charge -^-. 

Let us now start the perturbative computation of the excitations mass gap. I shall work 
in the Coulomb gauge VA = 0, so that A is ^-independent and the ground states are 

(4.45) 
(4.46) 



\i>> 


= n ^i°> 




x=even 


\x> 


= n ^io>. 




ir— odd 


The Coulomb potential reads 




fc^0 4Mn 2 


= ^N {x " v)2 -\ 



~ yl + ihv {N2 ~ l) - (447) 

The Schwinger Hamiltonian, rescaled by the factor e? a, may then be written as 



H = H + eH h = 4~( H C + H p ) (4.48) 

e L a 



with 



and 



^ = E[^r ! - I ^]^V(y) (4.49) 

x>y 

H h = -* E^+i 6 " 4 ^ - ^ e ~"V*+i) = -KR - L ) (4-50) 

X 

e =4^- (451) 



To the fourth order in e the perturbative expansion for the energy is given bj|j 



E* = ST + M 2) + ^f (4-52) 



4 Let us observe that for the ground state perturbative expansion the relative minus sign in the hopping 
Hamiltonian is irrelevant whereas it will play an important role in the calculation of the excited state energy. 
This is easily seen if one notices that at the second order in t one gets 

E^ = <^|(£-fl)_5* (R-L)\i,>= 



< ^l L T7nT^ R W > + < ^\ R -(o) L ^ > 

V 



? {0) -H ' r r ' eW-Ho 
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where 

4 0) =< V>IW > (4-53) 

4 2) =< mt-J^—H^ > (4-54) 

■% --Ho 

^ -< W H hTM H Hh F W H Hh F (0) H Hh W > + 
^ - -«0 ^ - -«0 iity - -HO 

" < ml ^ ~ W >< *\rf h J+ H h \1, > (4.55) 

h i, ~ H {% ~ -Ho)" 

and 1 — II^i is a projection operator which projects a generic state \s > on \tp > 

(l-n ip )\s>=\i>> . (4.56) 



The phases e tA in the hopping Hamiltonian are irrelevant until one reaches the N th order 
in the perturbative expansion. As a matter of fact, only at this order the transverse electric 
field becomes important. By direct evaluation one gets 

E^ } = -4N (4.58) 

E^ ] = 1927V (4.59) 

Ej, = —- 4Ne 2 + 192iVe 4 . (4.60) 

Next I compute the energies of the low-lying excitations to the same order in the per- 
turbative expansion. The lowest-lying excitations are degenerate and are created by the 
operators 



L(x) = ipte iA ip x+1 (4.61) 



R(x) = ^ x+1 e- lA ^ x (4.62) 



which move a particle one lattice spacing to the left and right respectively. When they act 
on the vacuum \ip > they create the two states 

\x,R> = R{x)\i> >= ^ x+1 e- %A ip x \iP > (4.63) 

\x,L> = Ux)\il> >= Vie^Vx+ilV' > (4-64) 

which are degenerate in energy (see Eg. (14.92(1 1. In order to resolve the degeneracy of these 
two states, one must find the first nontrivial order of the perturbative expansion in which 
the two states are mixed; this turns out to be the second order. It is easy to show that the 
second order of perturbative expansion is diagonal in the two states which are created by 
the parity odd and even operators 

#_(ar) = R(x)+L(x) (4.65) 

H+{x) = R(x) -L{x) . (4.66) 

In fact 

PH±(x) = ±H ± (-x) . (4.67) 

The parity odd operator creates the state with lower energy. Thus, it is this the state 
which will have finite energy gap in the continuum limit and which will correspond to the 
pseudoscalar boson of the continuum Schwinger model. In the continuum one gets this 
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massive excitation by applying the vector flux j — j® to the vacuum. The excitation 
created by H+ is the scalar excitation which in the continuum limit decouples from the 
spectrum. Since one is interested in the masses of these excitations, one has to compute 
their energies at zero momentum. 

On the lattice the construction of the pseudoscalar excitation at zero momentum is 
provided by 

N N N 

\ 0> = ^El^)>=7f ^H^xM>=- / =J2^ie iA ^ x+1 +4 +1 e- iA i; x M>= 

x—1 x—1 x—1 

N N 

= -7^Y,j 1 ^>=-7^Y,^ x ^>- ( 4 - 68 ) 

v x—1 v x—1 

One could also consider the scalar massive excitation 
1 N 1 N 

\<p> = = —= y^\ip(x) >= —= y^ h + (x)\iJj >= 

v X—1 v X—1 

N N 



The states \6 > and \(p > are characterized by a dimer-antidimer configuration; that is by a 
pair of particles followed or preceded by a pair of antiparticles. 

The energy of the state \9 >, at the second order in e, is given by 

E f ] = E Jj < %)K X ~ R )MR ~ L)\6{x) > (4.70) 

x,y 

where 

Ae = ^V (471) 

with Eg = E!l + A, A = | and 1 — lie a projection operator onto \6 >. A state obtained 
by applying the hopping Hamiltonian to the state \9 > has a dimer-antidimer pair; as a 
consequence its energy will be raised of 2A with respect to E\ ' . |j The computation of Eg 
(see section (|4.2.2[) ) leads to 



E ( g z> =-— N+ — = -AN + 8 (4.72) 

d 

N 2 
< O\HlAlH h \0 >= -^ - -^ = 16N - 32 . (4.73) 

Next one should compute 

E <? ] = ^E < 0(y)\HlA e HlA e H h AoH h \e(x) > + 
~4s E < 0{y)\HlA e H h \6{x) > ]T < 6(y)\HlA 2 e H h \0{x) > . (4.74) 



N 2 

x,y x,y 



5 One must be careful in applying one more time the hopping Hamiltonian. Since the Coulomb Hamilto- 
nian and the hopping Hamiltonian do not commute (see section 14.2.21 . there is a state created by appling 
two times the hopping Hamiltonian to the state \8 > which has not a gap of 3 A with respect to the ground 
state, but it has a gap of 5A-characterized by three particles followed or preceded by three antiparticles, 
that is, it is created by applying L 3 or B? to the ground state . 
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Using the matrix elements given in paragraph 1(4. 2. 2|) . one gets 



?W 



3 „ 9 



Ef> = —N ~ — = 192iV - 576 . (4.75) 

The pseudoscalar excitation energy, to the fourth order in e, is then given by 

E e = (§} + \) + (- m + 8)e 2 + (192JV - 576)e 4 . (4.76) 



Using again the matrix element given in appendix C, one easily obtains the energy E v 
of the state \ip >, i.e. 

^ = (32+7) + (~ 4iV + 24 ) g2 + ( 192N ~ 1600 ) e4 • ( 4 - 77 ) 



Taking the difference between Eg {E^) and £?</,, one gets the masses of the pseudoscalar 
and scalar excitations, 

^f- = \ + 8e 2 - 576e 4 , (4.78) 

^- = - + 24e 2 - 1600e 4 . (4.79) 

e\a 4 

Eas. (|4. 7814. 79(1 provide the behavior of m p and m s for small values of z = e 2 = , \ 4 . 



As expected from the continuum theory, the pseudoscalar particle mass is smaller than 
the scalar particle mass. Furthermore, since both quantities are intensive, there is no N 
dependence in Eas. l|4.78l4.79|l . 

Now let us compute the dispersion relation of the pseudoscalar and scalar excitations to 
the fourth order in perturbation theory. 

At the zero-th order one has no momentum dependence in the lattice pseudoscalar exci- 
tation since 

E e 0) (P) = JJ E < 8(v)\Ho\0(x) > e*(*-ri = ^ + A = ^ + \- (^) 

x-.y 

The momentum dependence starts only at the second order in e, where one gets 

E 2) (P) = Jf E < 0(y)\HlA 6 H h \e(x) > e***~* . (4.81) 

x,y 
Using the matrix elements given in paragraph 14.2.2(1 . one has 

1 9 

Ef\p) =- — (JV-4)- -r cos pa = -4N + 16 ~ 8 cos pa . (4.82) 

The dispersion relation, at the fourth order in perturbation theory, reads 

Ef\ P ) = i Y, < O(y)WtAeHlA 9 H h A H h \e(x) > e^"^ + 

x-.y 

~^E < 0{y)\HiK 9 H h \6{x) > e*<»-*> £ < 9{y)\rf h A$H h \0{x) > e*<*-«0 (4.83) 
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and leads to 

r-ili)/ \ 3 1 4 9 25 1 

e b (p) = -£E N - ^3 cos 2 P a + ^3 ( 2 C0S ^ a - cos pa ) ~ y A3 = 

= 192iV-32cos2pa + 256(2cospa-cos 2 pa)-800 . (4.84) 

Using Eas. H4.8()(l . (|4.81(l and (|4.84(l . the dispersion relation of the lattice pseudoscalar exci- 
tation is given by 

Ee(p) = Ef\p)+e 2 Ef\ P ) + e i Ef\p) + ...= 

= E e (0) + -(1 - cospa)e 2 + [—(1 - cos2pa) - —(1 - cospa) 2 ]e 4 = 
= Eg(0) + 8(1- cos pa)e 2 + [32(1 -cos 2pa) - 256(1 -cos pa) 2 ]e 4 . (4.85) 

Performing a small momentum expansion, one gets 

£ e (p) = E e (0) + (4e 2 + 64e 4 )(pa) 2 + . . . (4.86) 

The same procedure yields 



E v (p) = E^{p) + e 2 E^(p) + e 4 E^(p) 



= E e (0) - -(1 - cospa)e 2 + [^(l " cos2pa) + — - — (1 + cospa) 2 ]e 4 = 
= E 9 (0) - 8(1 - cospa)e 2 + [32(1 - cos2pa) + 1024 - 256(1 + cospa) 2 ]e 4 . (4.87) 

which, for small momentum gives 

E v (p) = E v (0) + (-4e 2 + 72e i )(pa) 2 + ... (4.88) 

The p 2 terms are consistent with the bosonic nature of the excitations. 

Ea. H4.86|l and Ea. l|4.88|l give us the masses of the pseudoscalar and scalar excitations 
from the curvature of the energy momentum relation 

P2 3 =(4e 2 + 64e 4 ) (4.89) 

2m% e L a A 

= (-4e 2 + 72e 4 ) (4.90) 



1 



2m^e 2 a 3 



4.2.1 Unperturbed energies of the dimer states 

In this section the energies of the excitations generated by applying the hopping Hamiltonian 
Hh to the ground state \ip >, are evaluated. The action of Hh on the state \tp > generates 
a dimer-antidimer pair. All the states with one dimer-antidimer pair, and with two dimer- 
antidimer pairs (obtained applying Hh twice), have the same energy gap with respect to the 
ground state. However, not all the states obtained applying (Hh) n , with n > 3, to \ij) > 
have the same energy. There are, for example, two possible energies, corresponding to a gap 
of 3A = 3e|a/4 and 5A, for the states with three dimer-antidimer pairs. 

The states with a dimer-antidimer pair are given by 

\Sl L >=H R<L {y)\i>>. (4.91) 
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Applying on these states the unperturbed Hamiltonian one gets 
H u \S 1 RtL > = ^.,l(2/)4^E(^ ± ^) 2 |^> = 



2 

X 

= K + ^fJ2^ i -J 1 ^y + s ^\ s kL> (4-92) 

X 

where y must be even if one applies R and must be odd if one applies L. One thus sees that 
all the states \S R L > have a mass gap of A. 



S ?R,l)(r,l) >= H R , L (z)H RiL (yM > (4-93) 



The states with two dimer-antidimer pairs are given by 

>(iJ,L)(ii,L) 
2 

Hu\Sf RtL){RtL) >= ^H R . L {z)H R . L {y) Y^{E X ± Sxz ± 8 xy f\^ > (4.94) 

X 

It is easy to show that all these states have an energy gap of 2A. 

Based on Ea. H4.92J) and Ea. l|4.94JI one might expect that the states with three dimer- 
antidimer pairs have an energy gap of 3A. This is actually true for all the states of this 
type except for the 2 states with three particles followed or preceeded by three antiparticles. 
These states have an energy gap of 5 A. If one considers the state 

\S RRR >=R(w)R(z)R(yM>, (4.95) 

one easily obtains that 

2 

H U \S 3 RRR > = ^R(w)R(z)R(y)J2(E x + 5 XW + S xz + 5 xy ) 2 \^ > = 

X 

~r ^OxwOxz ~r ^0 X w^xy \ ^0 X z^xy "r 

+ S xw + S xz + S xy )]\S 3 RRR > . (4.96) 

To generate a state with three antiparticles followed by three particles, one must have 

z = y - 2 (4.97) 

w = y - 1 (4.98) 

The energy gap for this state is 5A. The same result is obtained for the state \S\ LL >= 
L(w)L(z)L(y)\'>p > with z = y + 3 and w = y + 2, which contains three particles followed by 
three antiparticles. 

The reason for these different behaviours, lies on the fact that the commutation relation 
between Hq and Hh is 

[H ,H h ] = ^HH + ^fc-Y^WUi^k-tU-^JPi (4-99) 

If the second term in the r.h.s. was absent, all the states with n dimer-antidimer pairs would 
have an energy gap of nA. 



4.2.2 Perturbative matrix elements 



In paragraph I provide the matrix elements needed in the strong coupling expansion expres- 
sion of the mass gap. The six matrix elements that arise in the computation of the second 
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order bosonic excitation energies, are 



Mi = < ip\R{y)Rk e LL{x)\il> >= 

IN 1 

= "^(y- 1)^-^(1-^) (4.100) 

M 2 = < rl>\L(y)LA g RR(x)\il> >= 

= -^(Y-l)^-i(l-M (4.ioi) 

M 3 = < i>\R(y)LA RL{x)\ip >= 

1 JV 1 

= "^(y- 3)^ + ^(1-^) (4.102) 

M 4 = < ip\L(y)RA LR(x)\ij >= 

IN 1 

= "^(y- 3)^ + ^(1-*^) (4.103) 

M 5 = - < ^|i?(y)iA eJ Li?(x)|^ >= 

= —r($y,x-i + 8y, x +i) (4.104) 

M 6 = - < il)\L(y)RAeRL(x)\il) >= 

= -^(^,x-i + ^,x+i) (4.105) 

where Mi and M3 are defined for x and y odd, M2 and M4 for x and y even, M5 for x even 
and y odd, and A^ for x odd and y even. 

To compute the fourth order energy gap one needs also the matrix elements M[, i = 
1, ..., 6 obtained by replacing Kg with Ag in M{ 



M[ = M^ = M^ = M' 4 = —(--l)5 xv + —(l-5 xy ) (4.106) 

M' 5 = M' & = --^(d w ,x-i + ^,x+i)--T2(l-V»-i-Wi)(4-107) 



with the same constraints on x and y as before. 

One must then evaluate the twenty matrix elements arising when three Hr and three 
Hl are combined in all possible ways in the expression in the first line of Eq. 1)4. 74|) 

< 6(y)\(L - R)A (L - R)A e (R - L)A g {R - L)\8{x) > (4.108) 

Since hermiticity of the Hamiltonian requires that the elements obtained changing R with 
L are equal one has only to compute the ten matrix elements, 



M " = jjJ2 M 'i( x ^ = ^Y< < ^ R ^ RAf>RAeLAeLL ( x W > ( 4 - 109 ) 

x.y x.y 

M 'i = jfJ2 M %( x >y} = ^J2 < ^ R (y s > LA < )LA( > RA < >RL ( x M > ( 4 - n °) 

x,y x,y 

M 'z = ^J2 M 3^y) = ^J2 < ^ R ^ RAeLAeRAeLL ^ x ^ > ( 4 - m ) 

x.y x.y 

M " = ^^M^(x,y) = l^<^|i?( 2 /) J RA e LA e LA e i?L(x)|^> (4.112) 

x,y x,y 

< = jjj2 M s( x ^ = ^J2 < ^ R ^ LA ° RAeRA <> LL ( x W > ( 4 - 113 ) 

x.y x.y 

M e' = ^J2 M ^ x ^ = ^J2 < ^ R ^ LAeRA6LAeRL ( x ^ > ( 4 - 114 ) 



N ^ D v a ' N 

x,y x,y 
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M " = ^J2 M "( x >y) = -jjY. < ^ R ^ RAsLAsLAsLR ( x ^ > ( 4 - 115 ) 

x.y x,y 

x,y x,y 

M $ = ^E M 9'( a; 'y) = -^E < ^l^) iA ^A e i?A ( ,Li?(x)|V> (4.117) 

M "o = ^Y, M ^ x ^y^ = -^T, < ^ R ^ LAsLAsLAsRR ^^ > ■ ( 4 - 118 ) 



7V^ ^~'<" N 

x,y x,y 

where M"(x, y), ..., Mg(x, y) are defined for x and y odd, M"(x, y), ...,M" Q (x, y) are defined 
for x even and y odd. The results are 

M['(x,y) = --Ls x , y -^(N-2)(N-4)6 x , y + 

(5 y>x+ 2 + 6 yiX _ 2 )-^(N-4){l-5 x ,y) (4.119) 



4A 3 

1 

4A3 

1 

8A 3 

1 ,„ „w, „ , 1 



Mg(x,y) = - T ^(N-4)(N-6)S x , y (4.120) 

Mg(x,y) = MUx > y) = M' 5 '(x,y) = M^(x ) y) = -^(N-4)(N-6)d x ,y + 



and 



M?(x, y) = M£(x, J/) - M£(a:, V) = Mio(*, 2/) = 



< = -±--^L(N-2){N-4) (4.123) 



1 

8A 3 ' 
1 



Mg = M£ = M£ = Mg = Afg' = - ^-^ {N - 4) (N - 6) (4. 124) 



M? = M£ = M£ = M{' = +^3(iV-4)(JV-6) (4.125) 



4.3 The chiral condensate 



In the continuum Schwinger model, the phenomenon of dynamical symmetry breaking of 
the chiral symmetry, is due to the anomaly. The order parameter is the mass operator 
M(x) = ip(x)ip(x) which acquires a non zero vacuum expectation value, giving rise to the 
chiral condensate (H5] 

Xc =< i>(x)i>(x) >= - e —m c = - J-^=. (4.126) 

Z7T Z7T \/TT 



In this section I compute the lattice chiral condensate to the fourth order in pertur- 
bation theory. In the staggered fermion formalism (having put particles on even sites and 
antiparticles on odd sites), one has 

$(x)1>(x) — { -^(^ x - ~) (4.127) 

a I 

The lattice chiral condensate may be obtained by considering the mass operator 

JV 

M 



^-E(-l)^^ (4.128) 



Na 

x=l 
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where the extra minus sign is there to give the same sign to the lattice and continuum chiral 
condensates (see footnote 1 of this chapter) , and evaluating the expectation value of (|4.128|l 
on the perturbed states \p^, > generated by applying Hh to \ip >. One has 

|p >= \if> > +e\p\ > +e 2 |p 2 > (4.129) 

where 

\p\>=~H h \^> (4.130) 

14 >= ^ H hHh\i> > ■ (4.131) 

The lattice chiral condensate is then given by 

= < Pv ,|M| Pv ,> = < il>\M\j, > +e 2 < p\\M\p\ > +e i < p\\M\p\ > 
XL < P^ K > < V# > + e2 < V\ \p\ > +e 4 < V\ \V% > 

One gets the following expressions for the wave functions 

< ip\ip > = 1 (4.133) 

<pl\p\> = ^ (4.134) 

< P >1> ^ ^^ (4.135) 

and for the M operator 

< ij)\M\il) > = -^- (4.136) 

2a 

<p\\M\p\> = -^-(^-2) (4.137) 

<p%\M\p%> = -^L-(?L-4)-2-(N-3). (4.138) 

To the fourth order in e, the lattice chiral condensate is given by 

XL =--(£- ^~ 2 e 2 + ^) = - l -{\ - 32 £ 2 + 15366 4 ). (4.139) 

a A I\ z A* a 2 



4.4 Lattice versus continuum 



In this section I want to extract some physics from the lattice results I obtained; to do this, 
one should compare the answer of the strong coupling analysis of the lattice theory with the 
exact results of the continuum model. For this purpose one should extrapolate the strong 

.2 

coupling expansion derived under the assumption that the parameter z — \ 4 <C 1 to the 

region in which z» 1; this region corresponds to the continuum theory since e L a — ► 0, 
z — > oo. To make this extrapolation possible, it is customary |29j to make use of Pade 
approximants, which allow to extrapolate a series expansion beyond the convergence radius. 
Strong coupling perturbation theory improved with Pade approximants should be compared 
with the continuum theory. 

As we shall see, the gauge invariant strong coupling expansion here proposed, provides 
a very accurate estimate of the observables of the continuum theory, already at the second 
order in powers of z. This strongly suggests that, expanding around the gauge, parity and 
charge conjugation invariant ground state \ip >, leads to a perturbative series converging to 
the continuum theory faster then the one used in |(S3I Ifi7| . 
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Let us first compute the ratio between the continuum value of the meson mass m c = 
e c / 's/tt and the lattice coupling constant e^, by equating the lattice chiral condensate 1(4.139(1 . 
to its continuum counterpart ((4.126|) 

1/1 \ e 7 

- - - 32z + 1536z 2 =— m c . (4.140) 

a \ 2 J 2tt 

Of course, Ea. 14. 140(1 is true only when Pade approximants are used, since - as it stands - 
the l.h.s. holds only for z <C 1, while the r.h.s. provides the value of the chiral condensate 
to be obtained when z ~ oo. Using the relation 

a=(^)i- (4.141) 

4z e L 

one gets from (14.140(1 

m r ,40. i 2-7T , 1 9s , . , . 

^ = ^ 1 ^ ( -2~ 32Z + 1536Z) (4 ' 142) 

As in Ref. JH3J, due to the factor z*, the fourth power of 1(4.142(1 should be considered in 
order to construct a non-diagonal Pade approximant. Since the strong coupling expansion 
has been carried out up to second order in z, one is allowed to construct only the [0, 1] Pade 
approximant for the polynomial written in 1(4.142(1 . One gets 

mc 4 = 1 2tt 4 z_ 

K e L J 4t 2y e-r J l + 256z V ' 

Taking the continuum limit z — > oo one has 

(!^£)4 = ^^I (4144) 

Next let us compute the same ratio by equating the pseudoscalar mass gap given in 1(4.78(1 
to its continuum counterpart m c 

e\a f- + 8z-576z 2 J = m c . (4.145) 

Again Ea. 1(4. 145(1 is true only when Pade approximants are used. 

Dividing both sides of Eq. 14.145(1 by e^ and taking into account that 

e L a=(t.)i, (4.146) 

one gets 

— = (r) i (l+8^-576z 2 ), (4.147) 

eL 4z 4 

Taking the fourth power, as I did for the chiral condensate equation, and constructing the 
[1,0] Pade approximant for the r.h.s. of Ea. 14.147(1 one gets 

One may now take the limit z — > oo, obtaining 

(— ) 4 = J ■ (4-149) 

eL o 

Equating Ea. 1(4. 144(1 and Ea. 14. 149(1 one gets an equation for t which gives 

t = -?— = 1.049 (4.150) 

84e 7 
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which lies 4.9% above the exact value. It is conforting to see that the lattice theory gives 
a light velocity very close to (but greater than ) 1. Putting this value of t in Ea. H4.144|l or 
Eq. (|4.149|l one has 

771 

— = 0.609 (4.151) 

which lies 7.9% above the exact value -h=. 

It is worth to stress that I can reproduce very well the continuum results, even if I use 
just first order (in z) results of the strong coupling perturbation theory. 

A second way jSSj to extrapolate the polinomial in (|4.142[l is to equate it to the ratio 

( 1 + bz) x 
—-— — '— = 1 - 64z + 3072z 2 (4.152) 

(1 + az) x+ i 

so that one can take the limit z — > oo in (|4.142[1 . Setting t = 1 one gets 

^ = ^^ (4.153) 

&l eT a x+ i 

It is worth to say that this extrapolation works for a large range of values of the parameter 
x between and 6 giving an answer that lies between 6.9% below and 10% above the exact 
answer. In particular for x = 1 one gets — = 0.596 which lies 5.6% above the exact value. 
For x = | one gets — = 0.593 wich differs from the continuum value of the 5.1%. For 

x = To one nas on ly 0.7% of error. In the same way one can extrapolate the polinomial in 

(I4-T47I) 

(1 + bz) x+ i 

±-— — , —— = 1 + 32z - 2304z 2 (4.154) 

(1 + az) x x ' 

so that one can then take the limit z — > oo in (|4.147|l . Setting t — 1 one gets 

^ = -Wr + ^ (4-155) 

eL 44 a x 

for x — 1 one gets — = 0.617 which lies 9% above the exact value, against the 30% error 
of ref. 63 . For x = 0.5 one has — = 0.595 which lies only 5.4% above the exact answer. 

In the continuum Schwinger model the ratio between the masses of the scalar and the 
pseudoscalar particles, m s /m p , equals 2, since m s is the lowest eigenvalue of the continuum 
mass spectrum, which starts at 2m p |63| . I shall now evaluate this ratio using the lattice 
expression of m p and m s I|4.78I4.79|) . m s /m p is given by 













m s \ + 24z - 1600z 2 




h.s. 


of 






m p \ + 8z - 576z 2 


Expanding the r. 


(4. 


156| 


in power series of z 












—= = 1 + 64z - 6144z 2 


one may perforn: 


i the 


[1: 


■1] 


Pade approximant 












m s 1 + 160z 












m p 1 + 96z 


For z — ► oo 










m s 

— = 1-67 



(4.156) 



(4.157) 



(4.158) 



(4.159) 

which lies 16% below the continuum value. My result coincides with the one obtained in 
|63| at the same order in the perturbative expansion. 
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Taking the [1, 1] Pade approximants of Ea. H4.89|l and Ea. H4.90|l one can compute in an 
independent way the ratio between the masses of the scalar and pseudoscalar particles 

4 = -t^ («*> 

m% 1 — 16z 

For z — > oo 
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m 



p 



= 1.125 (4.161) 



which lies 43% below the exact value. This means that one should go to the next order 
in perturbation theory in Ea. l|4.89|l and Ea. 14.901) to reproduce better the continuum limit, 
since in Eq. (|4.89J) and Eq. (|4.90[) one has no zeroth order term. Another test (also suggested 
in [HSj) to check the validity of the lattice computations, may be performed by computing 
the quantity 

z 3 / 4 d / , i, p \ 1 z dm p 

This should equal 1/4 if the lattice theory has to reproduce the continuum result. 



D = Z \/4 + 8z ~ 32Z(1 " 1?6Z) ■ (4 ' 163) 



Using (|4.78|> one gets 



Constructing the [0, 1] Pade approximant and taking the z — ► oo limit, one has 

D = z 32 z =^= — = 0.182 (4.164) 

l + 176z 11 v ; 

which lies 27% below the desired 0.25. This coincides with the result obtained in (H3|, using 
a strong coupling expansion up to the order z 4 . The agreement of my results with the 
continuum theory is very encouraging, since only terms up to the order z 2 have been used 
in the gauge invariant strong coupling expansion proposed in this paper. 



4.5 Summary 



The investigation, that has been exposed, is aimed at constructing an improved strong cou- 
pling expansion for the lattice Schwinger model. I showed that, gauge invariance together 
with the discrete symmetries of the model, force the ground state to be a Neel state. More- 
over gauge invariance requires the ground state to have electric fields related to the charge 
density operator by the Gauss's law. Thus, for large e|, the ground state energy is of order 
e\ rather than 4r 63 ; as a consequence, the strong coupling expansion here constructed 
provides a very accurate extrapolation to the continuum theory already at the second order 
in z = . 4 4 . Furthermore, the Coulomb energy of elementary excitations is also affected, 
since — besides the Coulomb self energy — there is also an interaction energy with the charge 
of the ground state. 

The strongly coupled lattice Schwinger model is equivalent to a one-dimensional antifcr- 
romagnetic Ising spin chain with a Coulomb long range interaction, which admits the Neel 
state as an exact ground state. Due to this, in my approach, chiral symmetry is broken 
spontaneously in the lattice model already at zero-th order in z. I expect that chiral sym- 
metry breaking should persist for all coupling, since the critical coupling for D = 1 should 
be at e 2 L — 0. I conjecture that, for small e\, this behavior is a manifestation of a Peierls 
instability — the tendency of a one dimensional Fermi gas to form a gap at the Fermi surface. 
This happens with any infinitesimal interaction. 
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In the continuum theory the Schwinger model exhibits the Nambu-Goldstone phenomenon: 
global chiral symmetry is spontaneously broken, but no Goldstone boson appears in the spec- 
trum since the local current is anomalous |66j . On the lattice there is no anomaly due to 
the Nielsen- Ninomiya theorem |27| . 

The present analysis shows that in the strong coupling limit, chiral symmetry is sponta- 
neously broken at all orders, since the states \ip > and \x > mix only at a perturbative order 
comparable with the volume of the system. Being chiral symmetry replaced by the discrete 
symmetry representing invariance under tranlations by one lattice site, effects related to the 
breaking of the chiral symmetry on the lattice should come from the coupling between the 
two sublattices. This is manifest in my strong coupling calculations. 
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Chapter 5 

The two-flavor lattice Schwinger 
model 



In this chapter I study the 5C/(2)-flavor lattice Schwinger model in the hamiltonian formal- 
ism using staggered fermions. The existence of the continuum internal isospin symmetry 
makes the model much more interesting than the one-flavor case; the spectrum is extremely 
reacher, exhibiting also massless excitations and the chiral symmetry breaking pattern is 
completely different from the one-flavor case. I shall demonstrate [691 0U| that the strong 
coupling limit of the two-flavor lattice Schwinger model is mapped onto an interesting quan- 
tum spin model — the one-dimensional spin-1/2 quantum Heisenberg antiferromagnet. The 
ground state of the antiferromagnetic chain has been known since many years j^E] and its 
energy was computed in |M7| : the complete spectrum as been determined by Faddeev and 
Takhtadzhyan 30 using the algebraic Bethe ansatz. 

The two-flavor lattice Schwinger model with non-zero fermion mass m has been analysed 
in |70j in the limit of heavy fermions m 3> e 2 ; good agreement with the continuum theory 
has been found. 

There are by now many hints at a correspondence between quantized gauge theories 
and quantum spin models, aimed at analyzing new phases relevant for condensed matter 
systems ^3 El El EH E01 • Recently Laughlin has argued that there is an analogy between 
the spectral data of gauge theories and strongly correlated electron systems 25 . Moreover, 
certain spin ladders have been shown to be related to the two-flavor Schwinger model 71 . 

The correspondence between the SU(2) flavor Schwinger model and the quantum Heisen- 
berg antiferromagnetic chain provides a concrete computational scheme in which the issue 
of the correspondence between quantized gauge theories and quantum spin models may be 
investigated. Because of dimensionality of the coupling constant in (l+l)-dimensions the 
infrared behavior is governed by the strong coupling limit, and it is tempting to conjecture 
the existence of an exact correspondence between the infrared limits of the Heisenberg and 
two-flavor Schwinger models. I shall derive |39l I4()j results which support this conjecture. 
For example the gapless modes in the spectra have identical quantum numbers; within the 
accurancy of the strong coupling limit, the gapped mode of the two- flavor Schwinger model 
was also identified in the spectrum of the Heisenberg model. 

In this chapter I present a complete study |40| of the strong coupling limit of the two- 
flavor lattice Schwinger model. I firstly compute explicitly the masses of the excitations to 
the second order in the strong coupling expansion; this computation needs the knowledge of 
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the spin-spin correlators of the quantum Heisenberg antiferromagnetic chain. The continuum 
massless two-flavor Schwinger model does exhibit neither an isoscalar Oipift) nor an isovector 
{{[)a a il)) chiral condensate, since this is forbidden by the Coleman theorem 0S]- On the 
lattice I show that these fermion condensates are zero to all the orders in the strong coupling 

expansion. Moreover I find that the pertinent chiral condensate is < ip L ip L ip R tpft > an d 
I compute its lattice expression up to the second order in the strong coupling expansion. It 
should be noticed that, in abscence of gauge fields, the chiral condensate is zero, is different 
from zero only when the fermions are coupled to gauge fields. This can be viewed as the 
manifestation of the chiral anomaly in this model. I finally compare the lattice results with 
the ones of the continuum theory. 



5.1 Definition of the model 



The action of the 1 + 1-dimensional electrodynamics with two charged Dirac spinor fields is 



r 2 
S = J d 2 x Y^ati-lnV + iM^a - £a V" 



(5.1) 



The theory has an internal SUl(2) ® SUr(2)-Q.&vot isospin symmetry; the Dirac fields are 
an isodoublet whereas the electromagnetic field is an isosinglet. It is well known that in 
1 + 1 dimensions there is no spontaneous breakdown of continuous internal symmetries, 
unless there are anomalies or the Higgs phenomenon occurs. Neither mechanism is possible 
in the two-flavor Schwinger model for the SUl(2) ® SUn(2)-symmetry. isovector currents 
do not develop anomalies and there are no gauge fields coupled to the isospin currents. The 
particles belong then to isospin multiplets. For what concerns the U(l) gauge symmetry 
there is an Higgs phenomenon |flfi| . 

The action is invariant under the symmetry 

SU l (2)®SUr(2)®U v (1)®Ua(1) (5.2) 

The group generators act on the fermion isodoublet to give 



SU L (2) 
SU R {2) 

Mi) 

Ua(1) 



Ms) — (e l8a ^ PL U Mx) , Ta{x) — Mb (x)(e- ie ^ p ^) ha (5.3) 

Mx) — > (e i$ "^ PH U Mx) , ~a(x) — » Mx) (e-^^) b a (5.4) 

Mx) —> (e l9{x)1 U Mx) , i>Ux) — 4(x) {e-^^) ba (5.5) 

Mx) — (e^ 5l ) ab Mx) , Mx) — > i>\(x) (e- M75l ) bQ , (5.6) 



where a a are the Pauli matrices, 9 a , 9{x) and a are real coefficients and 

Pl = 5(1-75), Pb = \{1 + 1 5 ) ■ (5.7) 

At the classical level the symmetries (l5.3H5.6|l lead to conservation laws for the isovector, 
vector and axial currents 

£{x)r = Mx)rP R (YUMx) (5.8) 

f a {x)L = Mx)^Pl(yUMx) (5.9) 

Mx) = MxMUbMx) (5.10) 

%(X) = ^MMl^abMx) (5.11) 

It is well known that at the quantum level the vector and axial currents cannot be simul- 
taneously conserved, due to the anomaly phenomenon |2fi| . If the regularization is gauge 
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invariant, so that the vector current is conserved, then the axial current acquires the anomaly 
which breaks the L^(l)-symmctry 



«(*) = 2 



Z7T 



(5.12) 



The isoscalar and isovector chiral condensates are zero due to the Coleman theorem ^H]; in 
fact, they would break not only the 17a (1) symmetry of the action, but also the continuum 
internal symmetry SU L {2)®SU R {2) down to SU V {2). There is, however, a SU L {2)®SU R {2) 
invariant operator, which is non-invariant under the !7A(l)-symmetry; it can acquire a non- 
vanishing VEV without violating Coleman's theorem and consequently may be regarded as 
a good order parameter for the C/A(l)-breaking. Its expectation value is given by [721 I62| 



ih im nl\ * * /j/i * ' 



=>7 



< f >=< vi 'i>i 'r R 'r R ' >= {-r) 



' 47T 7T 



(5.13) 



It describes a process in which two right movers are anihilated and two left movers are 
created. Note that F, being quadrilinear in the fields, is actually invariant under chiral 
rotations of 7r/2, namely under the discrete axial symmetry 



1p a (x) — > -y 5 1pa{x) 1pa(x) ~* -Ipaix^B 



(5.14) 



As a consequence, this part of the chiral symmetry group is not broken by the non- vanishing 

vev of f inrra . 

As we shall see in section (|5.4|) . the lattice theory faithfully reproduces the pattern of 
symmetry breaking of the continuum theory; this happens even if on the lattice the SU(2)- 
fiavor symmetry is not protected by the Coleman theorem. The isoscalar and isovector 
chiral condensates are zero also on the lattice, whereas the operator F acquires a non- 
vanishing VEV due to the coupling of left and right movers induced by the gauge field. The 
continuous axial symmetry is broken explicitly by the staggered fermion, but the discrete 
axial symmetry Ij5.14|l remains. 



The action l|5.1f) may be presented in the usual abelian bosonized form |73|. Setting 

1 



^Pal^a 



-.e^d^ a , a = 1,2 



(5.15) 



the electric charge density and the action read 

jo =: ip\ipi + V4^2 := —j=d x {^i 



d A x 



a^ifl"*! + ift^afl"* 



2vr 



$9 



($1 



By changing the variables to 



$+ 
$_ 



1 

V2 
1 

V2 



($ x + $ 2 ) 
($i - $ 2 ) 



one has 



<W 2 



s= fd 2 x Qa„*+0"fc + + ia M *_a"4_ - ^ 



(5.16) 
(5.17) 

(5.18) 
(5.19) 

(5.20) 



The theory describes two scalar fields, one massive and one massless. $ + is an isosinglet 
as evidenced from Ea. (|5.16jl : its mass ras = \r^ e c comes from the anomaly Ea. H5.12JI 66 . 
Local electric charge conservation is spontaneously broken, but no Goldstone boson appears 



<S4 



because the Goldstone mode may be gauged away. <&_ represents an isotriplet; it has rather 
involved nonlinear transformation properties under a general isospin transformation. All 
three isospin currents can be written in terms of $_ but only the third component has a 
simple representation in terms of $_; namely 

a - a 3 

3 M (x) =: t(^)7p(y)aA(i) := 

- \ ■ ?i(i)7^iW -Mx^Mx) := (27r)^"ft / 4_ . (5.21) 

The other two isospin currents jn(x) and jn(x) are nonlinear and nonlocal functions of 
$_ [73]; a more symmetrical treatment of the bosonized form of the isotriplet currents 
is available within the framework of non abelian bosonization |74j . For the multiflavor 
Schwinger model this approach has been carried out in |75|. providing results in agreement 
with |73j . 

The excitations are most conveniently classified in terms of the quantum numbers of 
P-parity and G-parity; G-parity is related to the charge conjugation G by 

G = e™^C . (5.22) 

$_ is a G-even pseudoscalar, while <&+ is a G-odd pseudoscalar 

$_ : I PG = 1-+ (5.23) 

$+ : I PG = 0— . (5.24) 

The massive meson $ + is stable by G conservation since the action (|5.20|) is invariant under 
$ + — > -$ + . 

In the massive SU(2) Schwinger model — when the mass of the fermion m is small 
compared to e 2 (strong coupling) — Coleman [73J showed that - in addition to the triplet 
$_ (I PG = 1 I") the low-energy spectrum exhibits a singlet I PG = ++ lying on top of the 
triplet $_. In this limit the gauge theory is mapped to a sine-Gordon model and the low 
lying excitations are soliton-antisoliton states. When m — > 0, these soliton-antisoliton states 
become massless |76| : in this limit, the analysis of the many body wave functions, carried out 
in ref.|75], hints to the existence of a whole class of massless states with positive G-parity; 
P-parity however cannot be determined with the procedure developed in [76| . These are not 
the only excitations of the model: way up in mass there is the isosinglet I PG = , ($+), 
already discussed in ref. |73| . The model exhibits also triplets, whose mass — of order ms 
or greater — stays finite [7^]; among the triplets there is a G-even state q 

The Hamiltonian, gauge constraint and non- vanishing (anti-)commutators of the contin- 
uum two-flavor Schwinger model are 



H = dx 



! -E 2 (x) + J2l=i 4>l(x)a (id x + eA(x)) ^ a (x)\ (5.25) 

d x E(x) +£Li^JO*OlMaO~0 (5-26) 

[A(x),E(y)]=i5(x-v),{il>a(x)>'l>tiv)}= s <* s (*-v) • ( 5 - 27 ) 

A lattice Hamiltonian, constraint and (anti-) commutators reducing to l|5. 2515. 2615. 271) in 
the naive continuum limit are 

2 N N 2 

x—l x—1 a—1 

E x - E x -! + ^LVi,. + V4,,^2,x - 1 - , (5.28) 

[A x , E y ] = id XlV , ypa,x,i>l iy \ =8ab8xy ■ 
1 K. Harada private communication. 

85 



The fermion fields are defined on the sites, x — 1, ...,7V, the gauge and electric fields, A x 
and E X) on the links [x;x + 1], N is an even integer and, when N is finite it is convenient 
to impose periodic boundary conditions. When TV is finite, the continuum limit is the two- 
flavor Schwinger model on a circle | 68| . The coefficient t of the hopping term in (|5.28|l plays 
the role of the lattice light speed. In the naive continuum limit, ex, = e c and t = 1. 

The Hamiltonian and gauge constraint exhibit the discrete symmetries 



Parity P: 

A x — ■* -A- a -i, E x — ► -E^-i, ip a>x — > (-lfVa,-x, ^L — ■* (-1)*^I,_ X 

(5.29) 



• Discrete axial symmetry T: 

Ax " Ac+i, E x > E x+1 , l/j a>x ► 1p a ,x+l, ll>t tX > i>l,x+l 

• Charge conjugation C: 

A x ► -A x + 1 , E x ► -E x+1 , i> a . x ► ll>l iX+1 , Vl,x > *Pa,x+l 

• G-parity: 

A x — ► —A x+ i, E, x — ► —E x +i 

1pl,x ► ^2, z +l, 4>l, x " V>2,x+1 

^2,x > --01.X+1) V4,x ► -V'l.x+l • 



(5.30) 
(5.31) 

(5.32) 



The lattice two-flavor Schwinger model is equivalent to a one dimensional quantum 
Coulomb gas on the lattice with two kinds of particles. To see this one can fix the gauge, 
A x — A (Coulomb gauge). Eliminating the non-constant electric field and using the gauge 
constraint, one obtains the effective Hamiltonian 



H* 



H u + H p 



it 



e 2 T .a 



^ E 2 + ^iJ2p(x)V(x-y)p(y) 



x.y 



-Tn EE^Ul^ - ^a, X e- iA ll)a,x+l) 



x a—1 



where the charge density is 



p(x) = i>\j)l,x + Iptx&.x - 1 



and the potential 



JV-l 



n*-v) = ^£ 



^i2-Kn(x — y)/N 



4 sin 2 ^ 



(5.33) 



(5.34) 



(5.35) 



is the Fourier transform of the inverse laplacian on the lattice for non zero momentum. The 
constant modes of the gauge field decouple in the thermodynamic limit N — ► oo. 



5.2 The strong coupling limit and the antiferromagnetic 
Heisenberg Hamiltonian 



In the thermodynamic limit the Schwinger Hamiltonian I J5.33L rescaled by the factor e 2 a/2, 
reads 

H = H a + eH h (5.36) 



N(> 



with 

'{x-yf 



Ho = E 



N 



(x-v) 



p(x)p(y) , (5.37) 



x>y 

H h = -i(R-L) (5.38) 

and e = t/e 2 L a 2 . In Ea. l|5.38|) the right R and left L hopping operators are defined (L = RJ) 
as 

N N 2 TV 2 

a;— 1 a: — 1 a— 1 x — 1 a— 1 

On a periodic chain the commutation relation 

[R,L] = (5.40) 

is satisfied. 

I shall consider the strong coupling perturbative expansion where the Coulomb Hamil- 
tonian (|5.37|) is the unperturbed Hamiltonian and the hopping Hamiltonian (|5.38(l the per- 
turbation. Due to Ea. H5.34|l every configuration with one particle per site has zero energy, 
so that the ground state of the Coulomb Hamiltonian i|5.37|l is 2 N times degenerate. The 
degeneracy of the ground state can be removed only at the second perturbative order since 
the first order is trivially zero. 

At the second order the lattice gauge theory is effectively described by the antiferromag- 
netic Heisenberg Hamiltonian. The vacuum energy — at order e 2 — reads 

4 2) =< Hi / H h > (5.41) 

where the expectation values are defined on the degenerate subspace and n is the operator 
projecting on a set orthogonal to the states with one particle per site. Due to the vanishing 
of the charge density on the ground states of H q , the commutator 

[H ,H h ]=H h (5.42) 

holds on any linear combination of the degenerate ground states. Consequently, from 
Ea. (|5.41() one finds 

E [ q ] = -2<RL> . (5.43) 

On the ground state the combination RL can be written in terms of the Heisenberg Hamil- 
tonian. By introducing the Schwinger spin operators 

& = <x^M (5-44) 

the Heisenberg Hamiltonian Hj reads 

JV 

I 

x=l 

N 



N / i i 

Y j y-^L x R x --p{x)p(x + l) ] j (5.45) 



x=l 

and, on the degenerate subspace, one has 

N , _, s \ IN 



= E^^-l = Eh LA ■ ( 5 - 46 ) 



< Hj > , /^ i 

\x=l 
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Taking into account that products of L x and R y at different points have vanishing expecta- 
tion values on the ground states, and using Ea. (|5.46|l . Ea. Q5.43jl reads 

E^ ] = 4 < Hj > . (5.47) 

The ground state of Hj singles out the correct vacuum, on which to perform the perturbative 
expansion. In one dimension Hj is exactly diagonalizable 36, 34]. In the spin model a flavor 
1 particle on a site can be represented by a spin up, a flavor 2 particle by a spin down. The 
spectrum of Hj exhibits 2 N eigenstates; among these, the spin singlet with lowest energy is 
the non degenerate ground state \g.s. >. 

I shall construct the strong coupling perturbation theory of the two-flavor Schwinger 
model using \g.s. > as the unperturbed ground state. \g.s. > is invariant under transla- 
tions by one lattice site, which amounts to invariance under discrete chiral transformations. 
As a consequence, At variance with the one-flavor model |fil| . chiral symmetry cannot be 
spontaneously broken even in the infinite coupling limit. 

\g.s. > has zero charge density on each site and zero electric flux on each link 

p{x)\g.s.>=0 , E x \g.s. >= (x=l,...,N) . (5.48) 

\g.s. > is a linear combination of all the possible states with -^ spins up and -^ spins down. 
The coefficients are not explicitly known for general N. In chapter 3 I exhibited \g.s. > 
explicitly for finite size systems of 4, 6 and 8 sites. The Heisenberg energy of \g.s. > is 
known exactly and, in the thermodynamic limit, is |37l I3()j 

Hj\g.s. >= (-N In 2)\g.s. > . (5.49) 

Ea. (|5.49|) provides the second order correction Eq. l|5.47[l to the vacuum energy, E g . 



g.s. 



-AN In 2. 



There exist two kinds of excitations created from \g.s. >; one kind involves only spin 
flipping and has lower energy since no electric flux is created, the other involves fermion 
transport besides spin flipping and thus has a higher energy. For the latter excitations 
the energy is proportional to the coupling times the length of the electric flux: the lowest 
energy is achieved when the fermion is transported by one lattice spacing. Of course only 
the excitations of the first kind can be mapped into states of the Heisenberg model. 

In jHII] the antiferromagnetic Heisenberg model excitations have been classified. There it 
was shown that any excitation may be regarded as the scattering state of quasiparticles of 
spin-1/2: every physical state contains an even number of quasiparticles and the spectrum 
exhibits only integer spin states. The two simplest excitations of lowest energy in the ther- 
modynamic limit are a triplet and a singlet |30j : they have a dispersion relation depending 
on the momenta of the two quasiparticles. For vanishing total momentum (relative to the 
ground state momentum P g . s . = for — even, P g . s . = 7T for -j odd) in the thermodynamic 
limit they are degenerate with the ground state. 

In chapter 3 1 showed that even for finite size systems, the excited states can be grouped in 
families corresponding to the classification given in |3Uj . I explicitly exhibited all the energy 
eigenstates for N = 4 and N = 6. The lowest lying are a triplet and a singlet, respectively; 
they have a well defined relative (to the ground state) P-parity and G-parity - 1 h for the 
triplet and ++ for the singlet. Since they share the same quantum numbers these states can 
be identified, in the limit of vanishing fermion mass, with the soliton-antisoliton excitations 
found by Coleman in his analysis of the two-flavor Schwinger model. A related analysis 
about the parity of the lowest lying states in finite size Heisenberg chains, has been given 
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Moreover in [3D] a whole class — M.af — of gapless excitations at zero momentum was 
singled out in the thermodynamic limit; these states are eigenstates of the total momentum 
and consequently have positive G-parity at zero momentum. The low lying states of the 
Schwinger model also contain |76| many massless excitations with positive G-parity; they 
are identified pM HO] with the excitations belonging to Maf- The mass of these states in 
the Schwinger model can be obtained from the differences between the excitation energies 
at zero momentum and the ground state energy. The energies of the states \ex. > belonging 
to the class Maf have the same perturbative expansion of the ground state. Consequently, 
the states \ex. > at zero momentum up to the second order in the strong coupling expansion 
have the same energy of the ground state (|5.41l) . EeJ — — 4iVln2. To this order the mass 
gap is zero. Higher order corrections may give a mass gap. 



5.2.1 The low lying spectrum of QCD 2 

In section (|2.3.2[1 I already explained that QCD2 with colour group SU(2) and 1 flavor is 
also mapped onto the spin-1/2 quantum Heisenberg antiferromagnetic chain. This result 
was originally found by Langmann and Semenoff in |24| . It is worth stressing this fact here, 
and speculating about the model spectrum. 

I claim that the massless excitation spectrum of this non-abelian theory in the strong 
coupling limit is the same of the massless excitation spectrum of the two-flavor massless 
Schwinger model. Both the models exhibit an exact correspondence with the infrared limit 
of the Heisenberg antiferromagnetic chain. 

The quark-like excitations of QCD2 are gapless spin-1/2 kinks or spinons. They always 
appear in even number so that the only physical excitations have integer spin and an even 
number of spinons. This phenomenon is not exactly what one would call confinement. These 
quark-spinons are localizable objects and one can consider their scattering. There are no 
bound states of spinons in the spin-1/2 quantum antiferromagnetic chain j.'iOj . Spinons exist 
only in even number and this constraint is the spin model counterpart of the abscence of 
coloured states in QCD2- One flavor QCD2 exhibits only mesons in its spectrum and it 
is not possible to generate baryons. I can say nothing, at this level of analysis, about the 
massive excitations of lattice QCD2 which are created by the non abelian currents of the 
model acting on \g.s. >. These currents involve charge transport on \g.s. > and create 
electric flux. The non zero electric flux generates the masses of the excitations, with a 
mechanism similar to that of the two-flavor Schwinger model. 



5.3 The meson masses 



In this section I determine the masses for the states obtained by fermion transport of one 
site on the Heisenberg model ground state. My analysis shows that besides the G-odd 
pseudoscalar isosinglet with mass ms = bl-sJ^/tt, there are also a G-even pseudoscalar 
isotriplet 1 ^ and a G-odd scalar isotriplet l" 1 with masses of the order of ms or greater. 
The quantum numbers are relative to those of the ground state I P f = ++ for N/2 even 

!g.?. = °~~ for N / 2 odd- 
Two states can be created using the spatial component of the vector j 1 (x) Ea. H5.10Jl and 
isovector jl,(x) Eas. (|5.8l5.9[) Schwinger model currents. They are the G-odd pseudoscalar 
isosinglet I PG = and the G-even pseudoscalar isotriplet I PG = 1 '". The lattice 
operators with the correct quantum numbers creating these states at zero momentum, when 



<s<) 



acting on \g.s. >, read 

N 

S = R + L = ^2j 1 (x) (5.50) 

x=l 

N 

T + = (T_)t=i?( 12 )+L( 12 )=^ J |(x) (5.51) 



x=l 

N 



T {l = -^(R^+L^-RW-LW) = J2j 1 3 (x) ■ (5.52) 



V2 K 



x=l 



R( ab ) and L^ ab ^ in l|5. 5115.52(1 are the right and left flavor changing hopping operators 
(i> fc ) = (i?( afc ))t) 

N 
x=l 



The states are given by 



\S> = |0~ >= S\g.s. > (5.53) 

\T ± > = |1 _+ ,±1 >=T ± \g.s. > (5.54) 

|T > = |l-+,0>=To|ff.s. > . (5.55) 



They are normalized as 



<S\S> = <g.s.\S*S\g.s. >= -4 < g.s.\Hj\g.s. >=4JVln2 (5.56) 

<T + |T+> = ~(N+<g.s.\Hj\g.s.>) = lN(l-\n2) (5.57) 

and 

<T |T >=<T_|T_ >=<T+\T+ > . (5.58) 

In Ens. 15.5615.5715.581) < g. S .\g.s. >= 1. 

The isosinglet energy, up to the second order in the strong coupling expansion, is Es = 



4°)+ c 24 2 ) with 



4 „, _ <|H!>_ 1 , ( , 59) 



<S\HJA s H h \S> 
<S\S> 



■,(2) _ <■ Jl-Hfa 



4 = dc . (5-60) 



-h-s = - (cTT^ — an d 1 — lis a projection operator onto |5 >. On \g.s. > 

E s —Hq 

[^o,(n s ^) n 5] = (n + l)(n s ^r, (n = 0,l,...), (5.61) 

holds; Ea. l(5.60|l may then be written in terms of spin correlators as 



E (2) = E (2) , 4 _ Ex=l <g-S-|^--S , x+2-||g.S.> 

S 9S - < <?.s.|Hj|<7.s. > 



(5.62) 



One immediately recognizes that the excitation spectrum is determined once < g.s.\S x ■ 
S x +2\g-s. > is known. Equations similar to Ea. H5.62|) may be established also at a generic 
order of the strong coupling expansion. 

At the zeroth perturbative order the pseudoscalar triplet is degenerate with the isosinglet 
Ej, — E s =1. Following the same procedure as before one may compute the energy of 
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the states ()5.54|) and (|5.55(l to the second order in the strong coupling expansion. To this 
order, the energy is given by 



EP = 41 - A DS {T) 

4 < g.s.\Hj\g.s. > +5J2 x =i < g.s.\S x ■ S x+2 - \\g.s. > 
N+<g.s.\Hj\g.s.> 



(5.63) 



where in terms of the vector operator V = J2x=i $s ^ &x+i, one can write Ajjs(T) as 

A (T x 19 < g.s.HV^lg.s. > + < ff. S .|(^ 2 ) 2 |g. S . ;■■ 

A ^ (T±) = 12 N + <g,.\H J \g,.> (5 ' 64) 

a. s( to) = u lz sW wY s ;l ■ ^ 

N+ < g.s.\Hj\g.s. > 

The VEV of each squared component of V on the rotationally invariant singlet \g.s. > give 
the same contribution i.e. Ajjs{T±) = Ads(Tq): the triplet states (as in the continuum 
theory) have a degenerate mass gap. This is easily verified by direct computation on finite 
size systems; when the size of the system is finite one may also show that Ajjs is of zeroth 
order in N. 

The excitation masses are given by ms — ^^{Es — E g . s .) and mj = '-^(Et — E gs ). 

Consequently, the (iV-dependent) ground state energy terms appearing in E s and E T can- 
cel and what is left are only N independent terms. This is a good check of our computation, 
being the mass an intensive quantity. 

In principle one should expect also excitations created acting on \g.s. > with the chiral 
currents, in analogy with the one flavor Schwinger model where, as shown in ref. |63j . the 
chiral current creates a two-meson bound state. The chiral currents operators for the two 
flavor Schwinger model are given by 

f( x ) = ^(x^iPix) (5.66) 

f a {x) = ? a (o;)7 5 (fW b (x) . (5.67) 

The corresponding lattice operators at zero momentum are 

S 5 =R-L = J2tj 5 (x) (5-68) 

Tl =(^) t = J R (12) -i (12) =Ef=iJtW (5-69) 

T o =^(fi (11) -^ (11) -fi (22) +i (22) ) = E, A L 1 J3 5 (^) • (5-70) 

(5.71) 

The states created by l|5. 6815.6915. 70[) when acting on \g.s. >, are 

\S 5 > = |0 ++ >= S 5 \g.s. > (5.72) 

|T|> = \1+~,±1 >=Tl\g.s. > (5.73) 

|T 5 > = |i+-,0>=T 5 | 3 . s . > . (5.74) 

They are normalized as 

<S 5 \S 5 > <g.s.\S^S 5 \g.s. >= -4 < g.s.\Hj\g.s. > = 47Vlog2 (5.75) 

~(N+ < g.s.\Hj\g.s. >) = -, 



<Tl\Tl> = ~(N+<g.s.\Hj\g.s.>) = ~N(l - !,,->) (5.76) 



and 



< T a |T a >=< Tl\Tl >=< T^_\Tl > . (5.77) 
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Following the computational scheme used to study |5 > and \T >, one finds that the 
state \S 5 > is not a particle- like excitation since its energy depends on the volume of the 
system and this would lead to an extensive mass; consequently, there is no ++ massive 
singlet in agreement with |73]. For the triplet \T 5 > one gets 

E$ = 1 (5.78) 

E (2) = E { 2 ) , T,x=i < 9-S-\S x ■ S x+2 - \\g.s. > -4 < g.s.\Hj\g.s. > 
T5 9S ' N+<g.s.\Hj\g.s. > 



Now I can compute the mass spectrum up to the second order in the strong coupling 
expansion. Using Eq. (|3.159f) . the isosinglet mass reads as 



^ = 1x1.™ •- 
e 2 a 2 



1.9509 e z . (5.80) 

For what concerns the isotriplet mass, since the double sum in Eq. (|5.63|l is given by 



a (t\ o <g-s.\V-V\g.s. > 

ADs(T) - 8 N+<g. S .\H J \g,.> ' ^ 



using Ea.q3.180f). one gets 



e? a 2 



= - + 0.0972 e 2 . (5.82) 



'V- 



The existence of massive isotriplets was already noticed in |76| . and their mass in the 
continuum theory was numerically computed for various values of the fermion mass. In 
particular there is a G-parity even isotriplet with mass approximately equal to the mass of 
the isosinglet 0~~. 

The mass of the \T§ > isotriplet is 

H^i = - + 4.4069e 2 . (5.83) 

e z a 2 

Equations l|5.80f) . (|5.82() and l|5.83fl provide the values of ms, mx and m T s for small 
values of z = e 2 = -4-^ up to the second order in the strong coupling expansion. Whereas 

(|5.82() is only approximate (|5.80|l and H5.83(l are exact at the second order in the e expansion. 
In section (15.51) I shall extrapolate these masses to the continuum limit using the standard 
technique of the Pade approximants. 



5.4 Chiral condensate 



In the following I shall first prove that also on the lattice either the isoscalar (ipi/)) or the 
isovector (^p(x a ip) chiral condensates are zero to every order of perturbation theory. This 
should be verified by explicit computation since on the lattice the symmetry SUl(2)®SUr(2) 
is already broken by introducing staggered fermions; thus, there is no symmetry to prevent 
the formation of such chiral condensates. In the continuum theory, instead, the breaking of 
the SUl{2) <g> SUr(2) down to SUv(2) is prevented by the Coleman theorem [4*3*|. 

In the staggered fermion formalism the isoscalar condensate is given by 
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it is obtained by considering the mass operator 

N 



M = 2^ J2(~ i y^L^* + 4,^2,x) (5.85) 

x=l 

and evaluating its expectation value on the perturbed states \p g . s . > generated by applying 
Hh to \g.s. >. To the second order in the strong coupling expansion, one has 

\p g . a . >= \g.s. > +e\p].s. > +e>L- > + • ■ ■ (5-86) 

where 

\Pg...> = -H h \g.s.> (5.87) 



2 



\pIs.> = -^H h H h \g.s.> . (5.88) 



To the fourth order, (|5.85|l is given by 

_ <P g .s.\M\ Pg . s . > < g.S.\M\g.s. > +e 2 <pJ. a .|M||>J. a . > +e* < P 2 gs \M\p 2 gs ^ > 



<Pg.s.\p g .s. > <g.s.\g.s. > +e 2 < p£. a >£. a . > +e 4 <P^. s .|p^. s . > +••■ 

(5.89) 
It is very easy to see that Xisos. is zero to all orders in the strong coupling expansion. Let 
us introduce the translation operator 

f = e lpa ; (5.90) 

using 

TMf- 1 = -M (5.91) 

THnf- 1 = H h (5.92) 

and 

f\g.s.>=±\g.s.> (5.93) 

one gets order by order in the strong coupling expansion in Eq. (|5.89l) 

Xisos. = -Xisos. ■ (5.94) 

In Eq. H5.93[l the + appears when N/2 is even and the — when N/2 is odd. 

The isovector chiral condensate is given by the expectation value of the operator 

N 

2Na 



i N 
^^Ef-^t/^M (5-95) 



x = l 

on the perturbed states \p g . s . >■ Taking into account that 

TT.f- 1 = -S (5.96) 

one gets 

Xisov. Xisov. ; ^O.y/J 

also the isovector chiral condensate is identically zero. 

In the continuum there is, as evidenced in section l|5.1|l . only a non- vanishing chiral 
condensate associated to the anomalous breaking of the Ua{1) symmetry J75J|H5]. Since 

fiMW=fW^fW , (5-98) 



<):; 



the pertinent operator is F = ip L ip L ip R tp R ; its expectation value has been computed in 
[72~] [H2] and is given by 

< ^MVgV? >= (|^ e 2 = (g) 2 ™ 2 . (5.99) 



On the lattice one has 

^(x)^(*) — > ~«^a., x - < a+1 V«,«+i + 4 Q) " 4 a) ) • (5.100) 

The factor l/2a is due to the doubling of the lattice spacing in the antiferromagnctic bipartite 
lattice. Upon introducing the occupation number operators n x a = ipax^Pa,x, the umklapp 
operator F is represented on the lattice by 

^ = - I ^E^{(«i 1) -4i)(«i 2) -0 + (4 1) -4 1) )(4 2) -4 2) )} • (s.ioi) 

x—l 

The strong coupling expansion carried up to the second order in e = / 2 , yields 

< ^ >= <p g . a .|F|p g . 8 .> = <g.s.\F\g.s. >+£ 2 <p g , 8 .|F|p g , a , > 1Q 

<Pg.s.b ff . s . > <. g-s- Is- s. > +e 2 < P g . s .\pl. s . > 

for the lattice chiral condensate. Since 

<g.s.\g.s.> = 1 (5.103) 

< pUKs. > = -4<9-s.\Hj\g.s.> , (5.104) 

and taking into account that 

<g.s.\F\g.s.> = -^ < g.s.\Hj\g.s. > (5.105) 

<P 1 g .s.\F\p 1 9 .s. > = -^(-2<g. S .\(H J ) 2 \g.s.>-^<g. S .\H J \g.s.>+^N 

2 N 



■J2<g. S .\S x -S x+2 -^\g.s.>) , (5.106) 

" x=l 

from Eas. (|5.49|l and l|3.167fl . one gets 



< F >= 4(0.0866 - 0.4043e 2 ) . (5.107) 

a z 

A nonvanishing value of the lattice chiral condensate is due to the coupling — induced by 
the lattice gauge field — between the right and left fermions. This is the relic in the lattice 
of the V a(1) anomaly in the continuum theory. 



5.5 Lattice versus continuum 



I now want to compare our lattice results with the exact results of the continuum model; to 
do this, one should extrapolate the strong-coupling expansion derived under the assumption 
that the parameter z = e 2 = -J-^ <C 1 to the region in which 2>1; this corresponds to take 

the continuum limit since e L a — > when z — ► oo. To make the extrapolation possible, 
it is customary to make use of Pade approximants, which allow to extrapolate a series 
expansion beyond the convergence radius. Strong-coupling perturbation theory improved 
by Pade approximants should then provide results consistent with the continuum theory. 
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As we shall see the strong-coupling expansion derived in this chapter provides accurate 
estimates of the meson masses, already at the first order in powers of z. 

Let us now evaluate ms and the lattice light velocity t. I first compute the ratio between 
the continuum value of the meson mass ms = \ — e c and the lattice coupling constant e^ by 
equating the lattice chiral condensate, Eq. l|5.107|l . to its continuum counterpart Eq. (|5.99|l 

4(0.0866 - 0.4043z) = (— ) 2 ml . (5.108) 

a 1 Att 

Ea. H5.1U8f) is true only when Pade approximants are used since, as it stands, the left hand side 
holds only for z«l, while the right-hand side provides the value of the chiral condensate 
to be obtained when z = oo. Using 

j_ 

a = — ^ , (5.109) 

e L zi 

one gets from Eq. (|5.108[) 

(^) 2 = (^) 2 -(0.0866-0.4043z) . (5.110) 

As in Refs. |63ll61| . due to the factor z? , the second power of Eo. (|5.110|l should be considered 
in order to construct a non diagonal Pade approximant. Since the strong coupling expansion 
has been carried out up to second order in z, one is allowed to construct only the [0, 1] Pade 
approximant for the polynomial written in Eq. lj5. 110(1 . One gets 

ir) 1 -^!^ • (5 ' m) 

e_L e 7 t z 1 + 9.3371z 

and, taking the continuum limit z — > oo, one finds 

^ 4. 0,0008 

K e L ' K ei ! t 2 v ' 



Next I compute the same mass ratio by equating the singlet mass gap given in Eq. ((5.80(1 
to its continuum counterpart ms 

e 2 L a{- + 1.9509z) = m s . (5.113) 

Again, Eo. 1(5. 113(1 is true only when Pade approximants are used. Dividing both sides of 
Eq. ((5. 113(1 by e^ and taking into account that 

i 

e L a= t -\ (5.114) 

z* 

one gets 

— = ^r(l + 1.9509z) . (5.115) 

Taking the fourth power and constructing the [1, 0] Pade approximant for the right hand 
side of Eq. 1(5.1151) one has 

( — ) 4 = -(^ + 0.9754z) ; (5.116) 

eL z 16 



when z — > oo, Eo. ((5. 116(1 gives 



,ross 4 = t 2 Q g754 _ (5n?) 

eL 
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The numerical value of the hopping parameter i, determined if one equates Ea. (|5.117J) 
and Ed. H5.112p . is 

4-7T 

i = —0.1692 = 1.1940 (5.118) 

and lies 19% above the exact value. Putting this value of t in Eq. (|5.112|l or Eg. (|5. 1171) one 
gets 

— = 1.0969 (5.119) 
eh 

which lies 37% above the exact value y ^. It is comforting to see that the lattice reproduces 

in a sensible way the continuum results even if I use just first order (in z) results of the 
strong coupling perturbation theory. 

Using the value of t given in Eq. (|5. 118(1 one gets for the isotriplet mass Eq. 1(5.124(1 

— = 0.5143 . (5.120) 
By direct computation on an 8 sites chain one gets 

— = 1.3524 . (5.121) 

The discrepancy between Eq. (|5.120|l and Eq. (|5.121|l is mainly due to the approximation 
involved in the computation of < g.s.\V 2 \g.s. >. However, it is safe to believe that my 
lattice computation implies the existence of a massive isotriplet 1 h with a mass between 
the lower bound ((5.120(1 and the upper bound 1(5.121(1 . This is in agreement with the results 
provided for the continuum theory in |76j . 

Using a similar procedure one may also compute the mass of the triplet l" 1 . From 
Ea. (|5~S3l) one gets 

^^ = 1.3347 . (5.122) 

This triplet, being G-odd, is a scattering state of a singlet with a 1 h triplet, which are 
the fundamental excitations of the system. The mass of this l" 1 triplet should be larger 
than the mass of the massive 1 h triplet, which should be a scattering state of massless 1 l " 
triplets. 

Putting t = c = 1, i.e. e^a — -\ the lattice mass spectrum gets closer to its continuum 

24 

counterpart; for the isosinglet mass, one gets 

— = 0.9938 (5.123) 
e 

while for the isotriplet one gets 

— = 0.4695 . (5.124) 

Eq. ((5. 123(1 provides a value of the isosinglet mass lying 24% above the exact answer. Again 
the triplet mass is reproduced with lesser accuracy due to the random phase approximation 
used in the computation of the pertinent correlator; a better answer is given however by a 
direct computation on the 8 sites chain yielding the value 1.2346 for mx / e-L- 



5.6 Summary 



In this chapter I used the correspondence between the two-flavor strongly coupled lattice 
Schwinger model and the antiferromagnetic Heisenberg Hamiltonian (391 140 j to investigate 
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the spectrum of the gauge model. Using the analysis of the excitations of the finite size 
chains given in the chapter 3, I showed the equality of the quantum numbers of the states 
of the Heisenberg model and the low lying excitations of the two-flavor Schwinger model. I 
provided also the spectrum of the massive excitations of the gauge model; in order to extract 
numerical values for the masses, I explicitly computed the pertinent spin-spin correlators 
of the Heisenberg chain in chapter 3. Although the spectrum is determined only up to the 
second order in the strong coupling expansion the agreement with the continuum theory is 
satisfactory. 

The massless and the massive excitations of the gauge model are created from the spin 
chain ground state with two very different mechanisms: massless excitations involve only spin 
flipping while massive excitations are created by fermion transport besides spin flipping and 
do not belong to the spin chain spectrum. As in the continuum theory, due to the Coleman 
theorem [73], the massless excitations are not Goldstone bosons, but may be regarded as 
the gapless quantum excitations of the spin-1/2 antiferromagnetic Heisenberg chain |48|. 

In computing the chiral condensate I showed that, also in the lattice theory, the ex- 
pectation value of the umklapp operator F is different from zero, while both < %j)ip > and 
< ipa°"ip > are zero to every order in the strong coupling expansion. This implies that 
both on the lattice and the continuum the SU(2) flavor symmetry is preserved whereas the 
£7,4(1) axial symmetry is broken. The umklapp operator F is the order parameter for this 
symmetry, but being quadri-linear in the fermi fields, is invariant, in the continuum, under 
chiral rotation of 7r/2 and on the lattice under the corresponding discrete axial symmetry 
(|5.30(l (translation by one lattice site). This shows that the discrete axial symmetry is not 
broken in both cases. Our lattice computation enhance this result since the ground state of 
the strongly coupled two- flavor Schwinger model is translationally invariant. 

The pattern of symmetry breaking of the continuum is exactly reproduced even if the 
Coleman theorem does not apply on the lattice and the anomalous symmetry breaking is 
impossible due to the Nielsen-Ninomiya |27| theorem. At variance with the strongly coupled 
one-flavor lattice Schwinger model, the anomaly is not realized in the lattice theory via the 
spontaneous breaking of a residual chiral symmetry jH3> but, rather, by explicit breaking 
of the chiral symmetry due to staggered fermions. The non-vanishing of < F > may be 
regarded as the only relic, in the strongly coupled lattice theory, of the anomaly of the 
continuum two-flavor Schwinger model. It is due to the coupling induced by the gauge field, 
between the right and left-movers on the lattice. 
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Chapter 6 

The multiflavor lattice 
Schwinger models 



In this chapter I study the TV-flavor lattice Schwinger models in the hamiltonian formalism 
using staggered fermions. I illustrated in chapter 5 how the presence of a nontrivial SU(2)- 
flavor symmetry makes the spectrum much richer than the one-flavor model spectrum and 
changes drastically the chiral symmetry breaking pattern. 

The TV-flavor Schwinger models have many features in in common with four dimensional 
QCD: at the classical level they have a symmetry group C/l(TV) <S> Ur{N) — SUl(J\T) <S> 
SU R (Af) ®U v (1)®Ua(1) that is broken down to SU L {M) SU R (Af) <g> U v {\) by the 
axial anomaly exactly like in QCD 78 . The massless TV-flavor Schwinger models describe 
no real interactions between their particles as one can infer by writing the model action 
in a bosonized form. The model exhibits one massive and TV 2 — 1 massless pseudoscalar 
"mesons" [75] . 

On the lattice I shall prove that — at the second order in the strong coupling expansion — 
the lattice Schwinger models are effectively described by SU{Af) quantum antiferromagnetic 
spin-1/2 Heisenberg Hamiltonians with spins in a particular fundamental representation of 
the SU{N) Lie algebra and with nearest neighbours couplings. The features of the model 
are very different depending on if TV is odd or even. When TV is odd, the ground state energy 
in the strong coupling limit is proportional to e 2 L , the square of the electromagnetic coupling 
constant. In contrast, when TV is even the ground state energy in the strong coupling limit 
is of order 1. This difference arises from the proper definition of the charge density 

~ TV 

P( X )=^2^i,xi>a,x- Y (6.1) 

a=l 

where the constant TV/2 has been subtracted from the charge density operator in order 
to make it odd under the charge conjugation transformation. As a consequence, when TV 
is even, p(x) admits zero eigenvalues and the ground state does not support any electric 
flux, while when TV is odd the ground state exhibits a staggered configuration of the charge 
density and electromagnetic fluxes. 

In the continuum the Coleman theorem |43| prevents the formation of either an isoscalar 
chiral condensate < i/ji/j > or an isovector chiral condensate < i\)T a i\) > — where T a is an 
SU(AT) generator — for every model with an internal S'C/(TV)-flavor symmetry. We shall see 
that this feature should be reproduced on the lattice also for this class of models. 
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6.1 The continuum AA- flavor Schwinger models 

The continuum 5C/(A/")-flavor Schwinger models arc defined by the action 

f * _ i 

S = / d 2 x(J2 ^a( V + l^)^a - j^ V") (6.2) 

J a=l c 

where the Af fermions have been introduced in a completely symmetric way. Although the 
theory described by (|f>.2|> strictly parallels what has been shown in chapter 5 for the SU(2) 
model, I shall now report a detailed analysis both for the sake of clarity and to show that 
some difference appears between Af even and odd. 

The Dirac fields are an A/"-plet, i.e. transform according to the fundamental represen- 
tation of the flavor group while the electromagnetic field is an SU(Af) singlet. The flavor 
symmetry of the theory cannot be spontaneously broken for the same reasons as in the SU(2) 
case. The particles of the theory belong to SU(Af) multiplets. The action is invariant under 
the symmetry 

SU L (Ar)®SU R (Af)®U v (l)®U A (l) (6-3) 

The symmetry generators act as follows 

SU L (JV) : rj, a {x) — » (e i9 " Tapi -) ab M*) , M*) — Vfc {z^^^ba (6.4) 

SU R {N) : i/> a (x) — » {e if>aT " PH U Mx) , ^(i)-»^(i) {e- W ^ PL ) ba (6.5) 

U v (l) : Mx) — > (e lB(x)1 U M*) > 1&{*) — 4W (e- l9{x)1 ) b a (6.6) 

U A {\) : ^ a (x) — > (e"^ 1 )^ ^(s) , r/4(*) — > Vj(s) (e~ M75l W (6.7) 

where T Q are the generators of the SU(Af) group, 9 a , 9(x) and a are real coefficients and 

Pl = 1(1-1b), Pr = \(1 + J 5 ) ■ (6.8) 

At the classical level the above symmetries lead to conservation laws for the isovector, vector 
and axial currents 

&{*)* = ^_ a {x)^P R {T a ) ab ^ b {x) , (6.9) 

ft(x) L = ^ a (x)^P L (T a ) ab Mx) > (6.10) 

f{x) = ^ a (x)^l ab ^ b (x) , (6.11) 

%(x) = V»7 M 7 5 Wb(x) . (6.12) 

At the quantum level the vector and axial currents cannot be simultaneously conserved. If 
the regularization is gauge invariant, so that the vector current is conserved, then the axial 
current acquires the anomaly which breaks the symmetry Ua(X) |78| 

d^(x)=Af^-e fw F^(x) . (6.13) 

The isoscalar < tpip > and isovector < ipT a ip > chiral condensates are zero due to the 
Coleman theorem 031 > m fact they would break not only the UaO-) symmetry of the action 
but also the continuum internal symmetry SUl(N) ® SU R (AT) down to SUv(Af). There is 
an order parameter just for the breaking of the U a (1) symmetry [H21IZ21j the operator 

(6.14) 



<v£°...^\e 


-<»>=(|)7^)- 


iral rotation 




IpL -> 75^-L = - 


-ll>L , i>R -> Ib^R = i>R 



(6.15) 
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the operator (|6.14(l of course transforms as 

<^...^^K..^>^(-ir<^...^^...^> (6.i6) 

The umklapp operator is even under (|6.15(l when M is even and this implies that notwith- 
standing the fact that the continuous chiral rotations Ua(1) are broken by the non-zero VEV 
(|6.14|l . the discrete chiral symmetry l|6.15|l is unbroken. When Af is odd also the discrete 
chiral symmetry l|6.15|l is broken by the non-zero VEV l|6.14|l . 

The usual abelian bosonization procedure may again be applied provided that Af Bose 
fields are introduced [ZH1 IH3 IHD] 

: ^ al ^ a := ±=e^d v ^ a , a = 1, .. .Af (6.17) 

V 71 " 

The electric charge density and the action read 

M -. M 

o=l V 71 " a _ 1 

S = /^X>^*« + f|(X>«) a ) (6.19) 

J a=\ a=l 

The mass matrix is determined by the last term in Eq. (|6.19|l and must be diagonalized. The 
field degrees of freedom span the vector space on which the mass matrix is defined. The 
action must be expressed in terms of an orthonormal basis of field vectors, in order to have a 
properly normalised kinetic energy term. The original <I> a in Ea. (|6.19J) are orthonormal basis 
vectors, but they are not eigenvectors of the mass matrix. The mass matrix has one non- 
zero eigenvalue jfe with associated eigenvector — i= J2 a =i ^ a anc ^ an the other eigenvalues 
are zero. The remaining eigenvectors can be made orthonormal by the following change of 
variables 

$ a = 0£$ b (6.20) 

where the orthogonal matrices O^ are [HOj 

0\ = -±=(1,1,..., 1) , (6.21) 

VAf 



Ol = (1,1,..., -AT +1) , (6.22) 

6 ^/Ar(AT-iy ' ' v ; 

Ol = (l,l,...,-jy+2,0) , (6.23) 



Of = -±(1,-1,0, — ,0) . (6.24) 



V2 
In terms of these new fields $ a the action (|6.19(l reads 

AT 



S= f d 2 x{ 1 - J2 dfiadtfr - l^ 2 (^) 2 ) (6-25) 

•* a=l 

2 

where /i 2 = Af—. The action l|6.25|) describes Af non interacting fields, one massive and 
Af — 1 massless. The multiflavor Schwinger model can also be studied in the framework of 
non abelian bosonization 73] , where the relationship between isovector currents and bosonic 
excitations appears in a more symmetrical form |75| . 
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The Hamiltonian, gauge constraint and non- vanishing (anti-)commutators of the contin- 
uum TV-flavor Schwinger models are 



AT 



H= fdx[^-E 2 (x)+Y, ipt(x)a(id x + eA(x))i> a (x)] 

J ' a=l 

d x E{x)+ Y%=ii>Ux)i>a(x) ~ o 

[A(x),E{y)}=i6{x-y) , {i> a (x)> ^t(y)} = S ab S(x - y) 



(6.26) 

(6.27) 
(6.28) 



6.2 The lattice A^-flavor Schwinger models 



On the lattice the Hamiltonian, constraint and (anti-) commutators reducing to l|6. 2616. 2716. 28|) 
in the naive continuum limit are 



n N 

eia 



it 



N M 



x—l x—1 a— 1 



.v 



E x — E x _i 



\/^-x ? -^yl — ^Ox 7 y 



TV 



a=l 

\lpa,x,1p bty j = S ab S xy 



(6.29) 



The fermion fields are defined on the sites, x = 1, . . . , TV, gauge and the electric fields, A x 
and E x , on the links [x; x + 1], N is an even integer and, when N is finite it is convenient to 
impose periodic boundary conditions. When N is finite, the continuum limit is the TV-flavor 
Schwinger model on a circle 68 . The coefficient t of the hopping term in l|6.29|l plays the 
role of the lattice light speed. In the naive continuum limit, ei, = e c and 4 = 1. 



The lattice TV- flavor Schwinger model is equivalent to a one dimensional quantum Coulomb 
gas on the lattice with TV kinds of particles. To see this one can fix the gauge, A x = A 
(Coulomb gauge). Eliminating the non-constant electric field and using the gauge constraint, 
one obtains the effective Hamiltonian 



Hs — H u + H p 



^ + e J^J2 P (x)V(x-y)p(y) 



x.y 



it 



M 



-T J2J2^l-+ie iA ^a, x - Vi*e- iA Va,*+i) 



2« 



(6.30) 



where p(x) is given in Ea. H6.f(l and the Coulomb potential V(x — y) is given in Ea. l|4.22|l . 
The constant electric field is normalized so that [A, E] — i . The constant modes of the 
gauge field decouple in the thermodynamic limit N — > oo. In the thermodynamic limit the 
Schwinger Hamiltonian (|6.3U|) . rescaled by the factor e|a/2, reads 



H = H + eH h 



(6.31) 



with 



u " = T,i^w^- {x ~ vMx)p{v) 



(6.32) 

X >y 

H h = -i(R-L) (6.33) 

and e = t/e 2 L a 2 . In Eq. i|6.33[) the right R and left L hopping operators are defined [L = R') 
as 

N N M N M 

«=E^=EE^ ) = EEiA^ • ( 6 - 34 ) 
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On a periodic chain the commutation relation 

[R, L] = 
is satisfied. 



(6.35) 



When AT is even the ground state of the Hamiltonian (j6.32|l is the state \g.s. > with 
p(x) = on every site, i.e. with every site half-filled 



M 



H ^ 



.8. >= y|S-s. > 



(6.36) 



It is easy to understand that p(x) = on every site in the ground state by observing that 
the Coulomb Hamiltonian (|6.32|l is a non-negative operator and that the states with zero 
charge density are zero eigenvalues of l|6.32[l . \g.s. > is an highly degenerate state; in fact 
at each site x the quantum configuration is 



IK 

a=l 



10 > 



(6.37) 



The state l|6.37[l is antisymmetric in the indices a — 1, ... , 4r ; i.e. it takes on any orientation 
of the vector in the representation of the flavor symmetry group SU(J\f) with Young tableau 
given in fig. (|6.1|l . The energy of \g.s. > is of order 1, since it is non zero only at the second 



2 



Figure 6.1: The representation of SU(J\f) at each site when AT is even 
order in the strong coupling expansion. 



When M is odd the ground states of the Hamiltonian (|6.32() are characterized by the 
staggered charge distribution 

p{x) = ±\(-lV (6.38) 

since (|6.38(l minimizes the Coulomb Hamiltonian H6.32(l : one can have p(x) — +1/2 on the 
even sublattice and p(x) — —1/2 on the odd sublattce or viceversa. The electric fields 
generated by the charge distribution l|6.38|l are 



E x = ±-(-iy 



(6.39) 



102 



Since now 



H \g.s. >= — Ifif-s. > 



(6.40) 



the ground state energy is of order e 2 L . The states \g.s. > are highly degenerate since they can 
take up any orientation in the vector space which carries the representation of the SU(J\f) 
group with the Young tableaux given in fig. (|6.2() . 
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Figure 6.2: The representation of SU(M) at each site of the even sublattice and odd sub- 
lattice when N is odd 



Either when J\f is even or when TV is odd the ground state degeneracy is resolved at 
the second order in the strong coupling expansion. First order perturbations to the vacuum 
energy vanish. The vacuum energy at order e 2 reads 



E 



(2) 



-<H, 



n 



t 

h p(°) u 



H h > 



(6.41) 



where the expectation values are defined on the degenerate subspace of ground states and 
IT is a projection operator projecting orthogonal to the states of the degenerate subspace. 
Due to the commutation relation 



N- 1 



[Ho, H h ] = — —H h - 2 J2 [V(x -y)- V(x -y- 1)] (L y + R y )p(x) 



(6.42) 



Ea. (|6.41|) can be rewritten as 



E, 



(2) 



-2<RL> 



(6.43) 



On the ground state the combination RL can be written in terms of the Heisenberg Hamilto- 
nian of a generalized SU(J\f) antiferromagnet. By introducing as in chapter 5 the Schwinger 
spin operators 

S x = VL Wta (6.44) 

where T a are now the generators of the SU(Af) group, the SU(Af) Heisenberg Hamiltonian 
reads 

A" _1_ 
2Af f 



N .r N 

Hj = J2( 4 • S x+1 -- + VK - r p(x)p(x + 1)) = --J2 L X R 



(6.45) 



When J\f is even, on the degenerate ground states one has 

N . r , N 



< Hj >=< ^2(S X ■ S x+1 - — ) >=< -- Y^ L xRx > 



(6.46) 
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while when Af is odd one has 

< Hj>=<Y / (S x -S x+1 -—f r )>=<--Y / L x R x > ■ (6.47) 



a; — 1 x— 1 



Taking into account that the products of L x and R y at different points have vanishing 
expectation values on the ground states and using Eq. (|6.46|l or Eq.l |6.47|l . Ea. l|6.43|l reads 

E^ =A<H.,> . (6.48) 

The problem of determining the correct ground state, on which to perform the perturbative 
expansion, is then reduced again to the diagonalization of the SU(Af) Heisenberg spin-1/2 
Hamiltonian (|6.45(l . As I already pointed out in chapter 3, generalized SU(Af) antiferro- 
magnetic chains have not been yet analysed in the literature in such a detailed way as the 
SU(2) chains. Consequcntcly, the study of the lattice SU(Af) flavor lattice Schwinger mod- 
els become extremely complicated for Af > 2. Nonetheless, the computational scheme that 
I developed for the 17(1)- and 5C/(2)-flavor models in chapters 4 and 5, should work for a 
generic SU (Af )-flavor model. 

The ground state of the gauge models is very different when Af is even or Af is odd. When 
Af is even, the ground state \G.S. > of the spin Hamiltonian (|6.45|l is non-degenerate and 
translationally invariant, and since it is the ground state of the gauge model in the infinite 
coupling limit, there is no spontaneous breaking of the chiral symmetry for any SU(2Af)- 
flavor lattice Schwinger model. In contrast, when Af is odd, the ground state \G.S. > of the 
spin Hamiltonian Ij6.45|l is degenerate of order two and is not translationally invariant and 
consequently any S , f7(2A/'+l)-flavor lattice Schwinger model exhibits spontaneous symmetry 
breaking of the discrete axial symmetry. By translating of one lattice spacing \G.S. > one 
gets the other one. The A/"-flavor lattice Schwinger models excitations are also generated from 
\G.S. > by two very different mechanisms, that I already described for the two-flavor model 
in chapter 5. There are excitations involving only flavor changes of the fermions without 
changing the charge density p(x) which corrispond to spin flips in the SU(2) invariant model. 
These excitations are massless. At variance massive excitations involve fermion transport 
besides flavor changes and are created by applying to \G.S. > the latticized currents of the 
Schwinger models which vary the on site value of p(x). 

Very different is the case of the massive multiflavor Schwinger models [HI] • When Af is 
odd, the presence of a non-zero fermionic mass m removes the degeneracy and selects one 
of the two \G.S. > as the non degenerate ground state. When Af is even the ground state 
remains translationally invariant in the strong coupling limit e\ 3> m. In the weak coupling 
limit m S> e 2 L the discrete chiral symmetry is broken for every Af. 
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Summary and conclusions 



In this thesis I showed explicitly that the strong coupling limit of the multiflavor lattice 
Schwinger models is effectively described by quantum spin-1/2 generalized SU(Af) antifer- 
romagnetic Heisenberg chains. Exploiting the mapping existing between these gauge and 
spin models, exact results about antiferromagnetic chains have been used to analyse the chi- 
ral symmetry breaking on the lattice. In particular an interesting question has been how the 
effects of the axial anomaly appear in the lattice regularization. Since |27| fermion theories 
on a lattice have an equal number of species of left- and right-handed Weyl particles in the 
continuum limit, there is no axial anomaly in a lattice theory. 

I studied in detail the C/(l)-flavor and 5E/(2)-flavor models, archetypes for all SU(2Af + 
1) and SU(2Af) models. The results are very different depending on if the number of 
flavors is even or odd. The S [/(27V)- flavor Schwinger models exhibit non-degenerate and 
translationally invariant ground states in the strong coupling limit and their energy is of 
order 1. At variance, the SU(2Af + 1) models have two degenerate ground states, break the 
symmetry of translation by one lattice spacing and their energy is of order e 2 L . The difference 
between Af even and odd arises, since, also on the lattice, the charge density operator must 
be odd under charge conjugation; therefore the constant Af/2 should be subtracted from 
the charge density operator (|6.1|l . As a consequence, when the Af is odd, the ground state 
supports electric fluxes while this becomes impossible when Af is even. 

In the one-flavor model I proved that the discrete axial symmetry is spontaneously bro- 
ken on the lattice, reproducing the effects due to the axial anomaly in the continuum. In 
the infinite coupling limit the ground state of the gauge model is the ground state of an an- 
tiferromagnetic Ising chain with long range spin-spin Coulomb interactions. The Schwinger 
model excitations are created by acting on the ground state with the latticized currents. 
The excitation masses are determined in the strong coupling expansion by computing the 
excitation energies up to the fourth order in the strong coupling expansion. Using suitable 
Pade approximants, I extrapolated to the continuum the lattice results and I found that the 
lattice answers are in very good agreement with the continuum theory. 

In the two-flavor Schwinger model, I demonstrated that in the strong coupling limit the 
gauge model is effectively described by the spin-1/2 antiferromagnetic Heisenberg Hamilto- 
nian and the ground state of the gauge model in the infinite coupling limit is the ground 
state of the spin chain. Since the antiferromagnetic Heisenberg model is exactly solvable 
in one-dimension, it is very interesting to exhibit a gauge theory with the same low lying 
excitation spectrum. 

The antiferromagnetic spin chain admits spin-1/2 quarklike spinon excitations but has 
physical states with integer spin and an even number of spinons. Spinons behave like U(l) 
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quarks and they can be identified in conventional materials and models by high-energy 
spectroscopy and inconsistences in sum rules 25 . The two-flavor Schwinger model strictly 
parallels this scenario. The fundamental particles are spin-1/2 fermions but the spectrum 
exhibits only integer spin bosonic excitations [28| . I proved that the spinons have the same 
quantum numbers of the gauge model excitations. While spinons are the gapless excitations 
created from the spin chain ground state by spin flipping, the massive excitations of the 
gauge model are created from the spin chain ground state by fermion transport besides spin 
flipping and do not belong to the antiferromagnetic Heisenberg spectrum. The lattice mass 
spectrum, properly extrapolated to the continuum by using Pade approximants, reproduces 
in a satisfactory way the continuum results. The pattern of chiral symmetry breaking of the 
continuum model is exactly reproduced on the lattice, even if the Coleman theorem does not 
apply and the anomalous symmetry breaking is impossible 27 . In fact, both the isoscalar 
< ipt/) > and the isovector < ipa a ip > chiral condensates are zero to every order in the strong 
coupling expansion, while the expectation value of the umklapp operator F is different from 



The multiflavor Schwinger models are the gauge theories associated to the quantum spin- 
1/2 antiferromagnetic Heisenberg chains. It would be interesting to extend my analysis to 
QCD2 with different flavor and color groups. The number of colors J\fc determines the 
representation of the antiferromagnetic chain, since the spin is S — Nc /2. Varying J\fc one 
changes the spin S. In particular, by considering the SU(2)-HavoT case, i.e. QCD2 with two 
kinds of quarks, one finds that in the strong coupling limit the gauge model is effectively 
described by the antiferromagnetic Heisenberg chain with spin-S" representation. Since half- 
integer spin chains are expected to exhibit gapless excitations, while integer spin chains 
have a gap |48| . it would be an interesting problem to investigate if also QCD2 exhibits a 
gapless or gapped spectrum depending on if A/c is odd or even. Two dimensional QCD 
affords an excellent opportunity to study various dynamical questions of gauge theories, 
since many of its qualitative features are also valid in four dimensions. QCD2 still resists 
analytic solutions for general SU(Mc) colour groups, except in the planar limit Mc — * 
00 [H2IHSJ- In fact, 't Hooft |H2] analysed the QCD2 spectrum in the continuum exploiting 
a perturbation expansion with respect to 1/Afc and showed that the J\fc — > 00 limit, keeping 
gMc constant (g is the strong interactions coupling constant), corresponds to taking only 
the planar Feynman diagrams with no fermion loops |83| . The J\fc — » 00 limit of two flavor 
QCD2 should be compared with the S — > 00 of the antiferromagnetic Heisenberg chain. The 
study of the spin chain in the large S limit can be performed using the spin-wave theory 
and it would be fascinating to single out "gluon" and "hadron" degrees of freedom in the 
antiferromagnetic chain. 
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Pade approximants 



For the sake of clarity and to make the thesis self-contained, I define here what are Pade 
approximants 29 84. Although the method of Pade approximant has a long history, dating 
back more than one century, it is only in the 1960s that the procedure has become a part 
of the "bag of tricks" of the working physicist J2HI • When a power series representation of 
a function diverges, it indicates the presence of singularities. The divergence of the series 
reflects the inability of a polynomial to approximate a function adeguately near a singularity. 
The basic idea of summation theory is to represent f(z), the function in question, by a 
convergent expression. In Euler summation this expression is the limit of a convergent series, 
while in Borel summation this expression is the limit of a convergent integral. The difficulty 
with Euler and Borel summation is that all of the terms of the divergent series must be 
known exactly before the "sum" can be found even approximately. In realistic perturbation 
problems only a few terms of a perturbative series can be calculated as we have seen in 
chapters 4 and 5. Therefore, it is needed a summation algorithm, which requires as input 
only a finite number of terms of the divergent series. The Pade summation is a method 
having such a property. 

The idea of Pade summation is to replace a power series such as X^^Lo a nZ n by a sequence 
of rational functions of the form 

P£(z) = t?M°f (6-49) 

where one can choose Bq — 1 without loss of generality As it stands from Eq. (|6.49l) Pj^(z) 
is the ratio of a polynomial of degree N and one of degree M. The M + N + 1 coefficients 
j4o , A\ , . . . , A n , B\ , i?2 j • • • j Bm are determined in such a way that the first M + N + 1 terms 
in the Taylor series expansion of Pj^(z) match the first M + N + 1 terms of the power series 
S^Lo a nZ n - The rational functions Pm(z) are called Pade approximants. Constructing Pade 
approximants P^(z) is very useful; in fact, if X^Lo a n z ™ i s a power series representation of 
a function f(z), then in many instances Pj^(z) — ► f(z) as M,N — » oo, even if X^Lo ™ 2 " 
is a divergent series. The sequence of Pade approximants (|6.49|l with M = N is called the 
diagonal sequence. Let us sketch the computation of Pf (z). One has to expand Pf (z) in a 
Taylor series 

P?(z) = T ^^=A -A B 1 z + O(z 2 ) (z^O) (6.50) 

1 + Biz 

and comparing the series (|6.5U|) with the first two terms in the power series representation 

of f{z) = X) n =o a « z ™ one § e ts two equations: clq — Aq, ai = —AqBi from which 

r? = y^Tz (6 - 51) 

As it clear from this very simple example, to construct a Pade approximant Pj^(z) one does 
not need the full power series representation of a function, but just the first M + N + 1 terms. 
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Since Pade approximant involve only algebraic operations, they are more convenient for 
computational purposes than Borel summation, which requires one to integrate on an infinite 
range the analytic continuation of a function defined by a power series. Pade approximants 
often work quite well, even beyond their proven range of applicability. For an analysis of 
the convergence theory of Pade approximants, see Ref. 84 . 
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